(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[         0,          0]
NotebookDataLength[    151162,       3390]
NotebookOptionsPosition[    141116,       3185]
NotebookOutlinePosition[    141668,       3205]
CellTagsIndexPosition[    141625,       3202]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["ACE Model: Linear Stability Analysis", "Title",
 CellChangeTimes->{{3.78147159264887*^9, 3.781471597697865*^9}, {
  3.834510876960617*^9, 
  3.834510879383965*^9}},ExpressionUUID->"a73c02a0-3198-4cee-a01e-\
9ddb9794a7ac"],

Cell["\<\
Modeling the Public Health Impact of E-Cigarettes on Adolescents and Adults
Sara M. Clifton, Lucia M. Wagner\
\>", "Author",
 CellChangeTimes->{{3.781471986215602*^9, 3.781471991720435*^9}, {
  3.834510811806079*^9, 
  3.834510867638866*^9}},ExpressionUUID->"955d9765-5966-49fe-9395-\
d4bef08eb8c0"],

Cell["\<\
Department of Mathematics, Statistics, and Computer Science
St. Olaf College\
\>", "Department",
 CellChangeTimes->{{3.781471995254315*^9, 3.781472002954248*^9}, {
  3.8345108939604473`*^9, 
  3.834510905719079*^9}},ExpressionUUID->"a0f81619-815d-4d6f-89c2-\
4637a266c928"],

Cell[CellGroupData[{

Cell["\<\
Plug in right side of differential equations and denote variables\
\>", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.781476201147489*^9, 3.7814762303612967`*^9}, {3.7814765934921217`*^9, 
  3.781476599052796*^9}},ExpressionUUID->"28badf05-170b-427f-b6c8-\
cc4a2b50f70a"],

Cell[CellGroupData[{

Cell["Special case: \[Alpha] = 1", "Item",
 CellChangeTimes->{{3.781476265050036*^9, 3.781476288207602*^9}, {
   3.7814763275108747`*^9, 3.781476331218816*^9}, {3.781476402397532*^9, 
   3.7814764323909683`*^9}, {3.7814766701297197`*^9, 3.781476681448688*^9}, {
   3.781477784547036*^9, 3.78147781265799*^9}, {3.7814778624128103`*^9, 
   3.7814778946608353`*^9}, 3.781514150949218*^9, {3.78151422679467*^9, 
   3.781514290002342*^9}, {3.7815151562346563`*^9, 3.781515205910882*^9}, {
   3.834579124092229*^9, 
   3.834579130126539*^9}},ExpressionUUID->"763b3740-75ea-4318-bd0d-\
d9abc7667534"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"P", "[", 
   RowBox[{"x_", ",", "u_"}], "]"}], " ", ":=", " ", 
  RowBox[{"x", " ", "u"}]}], "\[IndentingNewLine]", 
 RowBox[{"DEs", " ", "=", " ", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", 
      RowBox[{"P", "[", 
       RowBox[{"c", ",", "uC"}], "]"}]}], " ", "+", " ", 
     RowBox[{"e", " ", 
      RowBox[{"P", "[", 
       RowBox[{"c", ",", "uC"}], "]"}]}], " ", "-", " ", 
     RowBox[{"c", " ", 
      RowBox[{"P", "[", 
       RowBox[{
        RowBox[{"1", "-", "c", "-", "e"}], ",", 
        RowBox[{"1", "-", "uC", "-", "uE"}]}], "]"}]}], " ", "-", " ", 
     RowBox[{"c", " ", 
      RowBox[{"P", "[", 
       RowBox[{"e", ",", "uE"}], "]"}]}]}], " ", ",", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", 
      RowBox[{"P", "[", 
       RowBox[{"e", ",", "uE"}], "]"}]}], " ", "+", " ", 
     RowBox[{"c", " ", 
      RowBox[{"P", "[", 
       RowBox[{"e", ",", "uE"}], "]"}]}], " ", "-", " ", 
     RowBox[{"e", " ", 
      RowBox[{"P", "[", 
       RowBox[{
        RowBox[{"1", "-", "c", "-", "e"}], ",", 
        RowBox[{"1", "-", "uC", "-", "uE"}]}], "]"}]}], " ", "-", " ", 
     RowBox[{"e", " ", 
      RowBox[{"P", "[", 
       RowBox[{"c", ",", "uC"}], "]"}]}]}]}], " ", 
   "}"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vars", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{"c", ",", "e"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"params", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{"uC", ",", "uE"}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814765794072323`*^9, 3.7814766117715178`*^9}, {
   3.781477330698406*^9, 3.781477334652973*^9}, {3.781514028937169*^9, 
   3.781514058092772*^9}, {3.8345112824226017`*^9, 3.834511375199737*^9}, {
   3.834511412186735*^9, 3.834511599555974*^9}, {3.834511816205503*^9, 
   3.834511834067465*^9}, {3.834511895561866*^9, 3.8345119116918583`*^9}, {
   3.834512395209323*^9, 3.834512395702973*^9}, {3.834512440325427*^9, 
   3.834512440487258*^9}, {3.834579440938689*^9, 3.8345794415017757`*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"f039ab92-f90a-48ab-a750-f5331c5afff0"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"c", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", "uC"}], "+", 
    RowBox[{"c", " ", "e", " ", "uC"}], "-", 
    RowBox[{"c", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
    RowBox[{"c", " ", "e", " ", "uE"}]}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "c"}], " ", "e", " ", "uC"}], "-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", "e", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "+", 
    RowBox[{"c", " ", "e", " ", "uE"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", "e", " ", "uE"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.8345116029019823`*^9, 3.834511913686256*^9, 
  3.834512445702838*^9, 3.834579228283432*^9, 3.8422744704807673`*^9},
 CellLabel->"Out[2]=",ExpressionUUID->"7e16a359-969f-4700-82a0-1f622b41ae60"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"numEqns", " ", "=", " ", 
  RowBox[{"Length", "[", "DEs", "]"}]}]], "Input",
 CellChangeTimes->{{3.781474692441131*^9, 3.7814747167748203`*^9}, {
  3.7814765273261337`*^9, 3.781476531800425*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"d0cdc25a-626b-4417-b973-3fb5a7bb112a"],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{3.834511611870761*^9, 3.8345119160760517`*^9, 
  3.834512449287446*^9, 3.8345792298433847`*^9, 3.842274475584813*^9},
 CellLabel->"Out[5]=",ExpressionUUID->"2376c29c-241c-40cd-b9c5-76b123ffc64d"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot direction field", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781477421345943*^9, 
  3.781477441545467*^9}, {3.834578949649909*^9, 
  3.8345789498863773`*^9}},ExpressionUUID->"ca20d7d9-cb7a-458b-a457-\
4d93f93bd2ac"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"StreamPlot", "[", 
    RowBox[{
     RowBox[{"DEs", "/.", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"uC", "\[Rule]", "cig"}], ",", 
        RowBox[{"uE", "\[Rule]", "ecig"}]}], "}"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"c", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"e", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Smoking", ",", " ", "Vaping"}], "}"}]}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"cig", ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"ecig", ",", "0", ",", "1"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, {
   3.7814764997404947`*^9, 3.781476509740304*^9}, {3.78147661985701*^9, 
   3.781476620330504*^9}, 3.781477464922221*^9, {3.834511623489006*^9, 
   3.834511707445655*^9}, 3.8345792153431683`*^9},
 CellLabel->"In[6]:=",ExpressionUUID->"55d4307b-2def-480b-909a-8fa64afb3df7"],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`cig$$ = 0.289, $CellContext`ecig$$ = 0.358, 
    Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, 
    Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 
    1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`cig$$], 0, 1}, {
      Hold[$CellContext`ecig$$], 0, 1}}, Typeset`size$$ = {
    450., {221., 227.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`cig$$ = 0, $CellContext`ecig$$ = 0}, 
      "ControllerVariables" :> {}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> StreamPlot[
        ReplaceAll[$CellContext`DEs, {$CellContext`uC -> $CellContext`cig$$, \
$CellContext`uE -> $CellContext`ecig$$}], {$CellContext`c, 0, 
         1}, {$CellContext`e, 0, 1}, 
        FrameLabel -> {$CellContext`Smoking, $CellContext`Vaping}], 
      "Specifications" :> {{$CellContext`cig$$, 0, 1}, {$CellContext`ecig$$, 
         0, 1}}, "Options" :> {}, "DefaultOptions" :> {}],
     ImageSizeCache->{505., {288., 295.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{
  3.834511708186625*^9, 3.83451191840098*^9, {3.834511969647971*^9, 
   3.83451197062269*^9}, 3.8345120024572372`*^9, 3.83451245177335*^9, {
   3.834579216882131*^9, 3.834579232537354*^9}, 3.8422744785900927`*^9},
 CellLabel->"Out[6]=",ExpressionUUID->"c1841426-9856-4c3e-957d-625ad9a262a7"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Solve for equilibria", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 
  3.7814722277070637`*^9}},ExpressionUUID->"792ae558-5cd9-4cb4-9098-\
f15685b3a8a0"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"equil", " ", "=", " ", 
  RowBox[{
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"DEs", " ", "\[Equal]", " ", 
      RowBox[{"ConstantArray", "[", 
       RowBox[{"0", ",", "numEqns"}], "]"}]}], ",", "vars"}], "]"}], "//", 
   "FullSimplify"}]}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, {
   3.7814764997404947`*^9, 3.781476509740304*^9}, {3.78147661985701*^9, 
   3.781476620330504*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"ec9f6ea3-d58f-43b9-bfb0-6e78c6ee5dc3"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"c", "\[Rule]", "0"}], ",", 
     RowBox[{"e", "\[Rule]", "1"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"c", "\[Rule]", "0"}], ",", 
     RowBox[{"e", "\[Rule]", "0"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"c", "\[Rule]", "1"}], ",", 
     RowBox[{"e", "\[Rule]", "0"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.8345124159976883`*^9, 3.8345124767271433`*^9, 
  3.842274486816431*^9},
 CellLabel->"Out[7]=",ExpressionUUID->"191d46de-c43f-4e33-a5c6-5f23a4a3cb58"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Linearize system (Jacobian)", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}},ExpressionUUID->"0dac2df4-84b9-41e3-93d3-\
7bfc871b7768"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"J", " ", "=", " ", 
  RowBox[{"D", "[", 
   RowBox[{"DEs", ",", 
    RowBox[{"{", "vars", "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, 
   3.781472428807068*^9, {3.7814766355579987`*^9, 3.781476646201336*^9}, {
   3.8345124977788897`*^9, 3.8345125010778427`*^9}, {3.83451350271008*^9, 
   3.834513515332168*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"f0335e49-428b-48c5-824c-82cbc6e4009e"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"-", "c"}], " ", "uC"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", "uC"}], "+", 
      RowBox[{"e", " ", "uC"}], "+", 
      RowBox[{"c", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
      RowBox[{"e", " ", "uE"}]}], ",", 
     RowBox[{
      RowBox[{"c", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
      RowBox[{"c", " ", "uE"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"-", "e"}], " ", "uC"}], "+", 
      RowBox[{"e", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}]}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "c"}], " ", "uC"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "+", 
      RowBox[{"e", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "+", 
      RowBox[{"c", " ", "uE"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", "uE"}], "-", 
      RowBox[{"e", " ", "uE"}]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.834512484043462*^9, 3.8345125025139313`*^9}, {
   3.834513503855713*^9, 3.834513516099617*^9}, 3.842274490466264*^9},
 CellLabel->"Out[8]=",ExpressionUUID->"7c7bc343-4570-46c0-9caa-fa4da1db0314"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"J", "//", "MatrixForm"}], " ", "//", "Simplify"}]], "Input",
 CellChangeTimes->{{3.834513529871863*^9, 3.8345135320437202`*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"22eaa91f-3917-46f7-9d83-baf60c961f48"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       RowBox[{"2", " ", "uC"}], "+", "uE", "-", 
       RowBox[{"2", " ", "c", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}], "-", 
       RowBox[{"e", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "uC", "+", 
          RowBox[{"2", " ", "uE"}]}], ")"}]}]}], 
      RowBox[{
       RowBox[{"-", "c"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "uC", "+", 
         RowBox[{"2", " ", "uE"}]}], ")"}]}]},
     {
      RowBox[{
       RowBox[{"-", "e"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}], 
      RowBox[{
       RowBox[{
        RowBox[{"-", "c"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "e"}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "uC", "+", 
          RowBox[{"2", " ", "uE"}]}], ")"}]}]}]}
    },
    GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.834513518640321*^9, 3.834513532671836*^9}, 
   3.84227449229935*^9},
 CellLabel->
  "Out[9]//MatrixForm=",ExpressionUUID->"2be04c0a-2bf7-4bf1-9ccb-\
1285bc3809c6"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Evaluate Jacobian at equilibria", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 
  3.7814729129771643`*^9}},ExpressionUUID->"6e8de002-011d-4fd8-8fff-\
c3e3d9524dcc"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Jstar", " ", "=", " ", 
  RowBox[{
   RowBox[{"J", "/.", "equil"}], " ", "//", "FullSimplify"}]}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, 
   3.781472428807068*^9, {3.7814729275631857`*^9, 3.781472961260558*^9}, {
   3.781473116004878*^9, 3.781473141513234*^9}, {3.781473253486924*^9, 
   3.7814732759163113`*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"f3d8facb-adf3-4ead-bfaf-960ba1c0f366"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"uC", "-", "uE"}], ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"1", "-", 
        RowBox[{"2", " ", "uC"}], "-", "uE"}], ",", 
       RowBox[{"1", "-", "uC", "-", 
        RowBox[{"2", " ", "uE"}]}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"2", " ", "uC"}], "+", "uE"}], ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "uC", "+", 
        RowBox[{"2", " ", "uE"}]}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"1", "-", 
        RowBox[{"2", " ", "uC"}], "-", "uE"}], ",", 
       RowBox[{"1", "-", "uC", "-", 
        RowBox[{"2", " ", "uE"}]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{
        RowBox[{"-", "uC"}], "+", "uE"}]}], "}"}]}], "}"}]}], "}"}]], "Output",\

 CellChangeTimes->{3.8345125119552107`*^9, 3.834513620140758*^9, 
  3.842274494933042*^9},
 CellLabel->"Out[10]=",ExpressionUUID->"cbfdf903-3001-4670-8722-a3c13c157685"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", " ", "/@", " ", "Jstar"}]], "Input",
 CellChangeTimes->{{3.781513774371643*^9, 3.781513812723462*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"592d3cf7-6d47-4849-8720-ea76a052f7e4"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   TagBox[
    RowBox[{"(", "\[NoBreak]", GridBox[{
       {
        RowBox[{"uC", "-", "uE"}], "0"},
       {
        RowBox[{"1", "-", 
         RowBox[{"2", " ", "uC"}], "-", "uE"}], 
        RowBox[{"1", "-", "uC", "-", 
         RowBox[{"2", " ", "uE"}]}]}
      },
      GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.7]}, 
          Offset[0.27999999999999997`]}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}}], "\[NoBreak]", ")"}],
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]], ",", 
   TagBox[
    RowBox[{"(", "\[NoBreak]", GridBox[{
       {
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "uC"}], "+", "uE"}], "0"},
       {"0", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "uC", "+", 
         RowBox[{"2", " ", "uE"}]}]}
      },
      GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.7]}, 
          Offset[0.27999999999999997`]}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}}], "\[NoBreak]", ")"}],
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]], ",", 
   TagBox[
    RowBox[{"(", "\[NoBreak]", GridBox[{
       {
        RowBox[{"1", "-", 
         RowBox[{"2", " ", "uC"}], "-", "uE"}], 
        RowBox[{"1", "-", "uC", "-", 
         RowBox[{"2", " ", "uE"}]}]},
       {"0", 
        RowBox[{
         RowBox[{"-", "uC"}], "+", "uE"}]}
      },
      GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.7]}, 
          Offset[0.27999999999999997`]}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}}], "\[NoBreak]", ")"}],
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]]}], "}"}]], "Output",
 CellChangeTimes->{3.834512514593965*^9, 3.8345136232722282`*^9, 
  3.84227449703018*^9},
 CellLabel->"Out[11]=",ExpressionUUID->"d67a3d1b-bfa8-4985-948d-677405659664"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Compute eigenvalues", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 3.7814729129771643`*^9}, {
  3.78147344349218*^9, 
  3.781473448224465*^9}},ExpressionUUID->"5c727f8b-5d37-4f37-a49b-\
e9bdeb9dc426"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eigen", " ", "=", " ", 
  RowBox[{
   RowBox[{"Eigenvalues", "/@", "Jstar"}], " ", "//", 
   "FullSimplify"}]}]], "Input",
 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
  3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
  3.781476066235634*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"ca284c07-f2f6-41ad-8388-0ee0d89debcf"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"1", "-", "uC", "-", 
      RowBox[{"2", " ", "uE"}]}], ",", 
     RowBox[{"uC", "-", "uE"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      RowBox[{"2", " ", "uC"}], "+", "uE"}], ",", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "uC", "+", 
      RowBox[{"2", " ", "uE"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"1", "-", 
      RowBox[{"2", " ", "uC"}], "-", "uE"}], ",", 
     RowBox[{
      RowBox[{"-", "uC"}], "+", "uE"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.83451255475968*^9, 3.834513630234558*^9, 
  3.842274499358429*^9},
 CellLabel->"Out[12]=",ExpressionUUID->"83c90568-ee4d-4a92-9ab4-91c05cf4f791"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Compute the trace and determinant", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 3.7814729129771643`*^9}, {
  3.78147344349218*^9, 3.781473448224465*^9}, {3.78147830643638*^9, 
  3.7814783120648193`*^9}, {3.781552785807474*^9, 3.7815527960775223`*^9}, {
  3.834579050006682*^9, 3.834579059173544*^9}, {3.842274512058419*^9, 
  3.84227451225204*^9}},ExpressionUUID->"24cd3001-ea45-47ab-bb96-\
5c6765007add"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"tr", " ", "=", " ", 
  RowBox[{
   RowBox[{"Tr", "/@", "Jstar"}], " ", "//", "FullSimplify"}]}]], "Input",
 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
  3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
  3.781476066235634*^9}, {3.781478323943038*^9, 3.781478329366476*^9}},
 CellLabel->"In[13]:=",ExpressionUUID->"e6436050-4ae4-4e7b-ba6d-22cd79f0763c"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"1", "-", 
    RowBox[{"3", " ", "uE"}]}], ",", 
   RowBox[{
    RowBox[{"-", "2"}], "+", 
    RowBox[{"3", " ", "uC"}], "+", 
    RowBox[{"3", " ", "uE"}]}], ",", 
   RowBox[{"1", "-", 
    RowBox[{"3", " ", "uC"}]}]}], "}"}]], "Output",
 CellChangeTimes->{3.834512576251791*^9, 3.842274515626686*^9},
 CellLabel->"Out[13]=",ExpressionUUID->"09110649-9296-4262-b952-8afd8ec53130"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"det", " ", "=", " ", 
  RowBox[{
   RowBox[{"Det", "/@", "Jstar"}], " ", "//", "FullSimplify"}]}]], "Input",
 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
  3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
  3.781476066235634*^9}, {3.781478323943038*^9, 3.781478329366476*^9}, {
  3.7814783600281353`*^9, 3.7814783639524193`*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"969a80a6-e6b9-4a50-bc52-e27f08dc7073"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"uC", "-", "uE"}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "uC", "+", 
        RowBox[{"2", " ", "uE"}]}], ")"}]}], ")"}]}], ",", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "uC", "+", 
      RowBox[{"2", " ", "uE"}]}], ")"}]}], ",", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"uC", "-", "uE"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{3.834512577623479*^9, 3.842274517298068*^9},
 CellLabel->"Out[14]=",ExpressionUUID->"648b9271-c75f-42f8-ad44-88fffa7aa3b8"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Visualize eigenvalues and stability", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 3.7814729129771643`*^9}, {
  3.78147344349218*^9, 3.781473448224465*^9}, {3.7815141129222813`*^9, 
  3.781514130440056*^9}},ExpressionUUID->"72b4bc04-7676-4749-91f2-\
3fd8b4907447"],

Cell[CellGroupData[{

Cell["Equilibrium is stable if Re[\[Lambda]]<0 for all eigenvalues.", "Item",
 CellChangeTimes->{{3.781476265050036*^9, 3.781476288207602*^9}, {
   3.7814763275108747`*^9, 3.781476331218816*^9}, {3.781476402397532*^9, 
   3.7814764323909683`*^9}, {3.7814766701297197`*^9, 3.781476681448688*^9}, {
   3.781477784547036*^9, 3.78147781265799*^9}, {3.7814778624128103`*^9, 
   3.7814778946608353`*^9}, 3.781514150949218*^9, {3.78151422679467*^9, 
   3.781514290002342*^9}, {3.7815151562346563`*^9, 
   3.781515205910882*^9}},ExpressionUUID->"8aa5c24f-f858-4959-8683-\
e9356406a9f1"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"eigen1", " ", "=", " ", 
   RowBox[{"eigen", "[", 
    RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"r1", " ", "=", " ", 
  RowBox[{"RegionPlot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen1", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "0"}], " ", "&&", 
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen1", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "<", "0"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"uC", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"uE", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"uC", ",", "uE"}], "}"}]}], ",", " ", 
    RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",", 
    RowBox[{"BoundaryStyle", "\[Rule]", "Blue"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
   3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
   3.781476066235634*^9}, 3.7815141879058037`*^9, {3.781514299900753*^9, 
   3.7815143389630337`*^9}, {3.7815147790959682`*^9, 3.781514886796404*^9}, {
   3.7815149248973227`*^9, 3.78151495836263*^9}, {3.781515003778427*^9, 
   3.781515016404488*^9}, {3.7815150850408707`*^9, 3.7815150882864017`*^9}, {
   3.781515314103553*^9, 3.781515321347374*^9}, {3.834514231169278*^9, 
   3.83451430577695*^9}, {3.8345143807965803`*^9, 3.834514381094407*^9}, {
   3.8422745397564497`*^9, 3.842274541424534*^9}, {3.842274702075144*^9, 
   3.8422747102614107`*^9}},
 CellLabel->"In[32]:=",ExpressionUUID->"f56a7601-52fa-4249-8237-485543004148"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJxtmHlUU3cexSNFW6rBjlrFgNatNlLcKKHqjH6RIFZlNDhKFSH2aBdqUJRR
adBOOSpxqVqLQa1FS9QyroQqrQUjLkXF4IJVCARkS14IgbAEXGthHufM/cVD
6j/3fE4k7977fu/7vjB8Wdz8T9wEAkFqD4GgS7nHJ2RqQxF9/2Njm0RtoKC5
c0ZwkgeMo/fcieQk+fT4xa5rKoeR+tPfOlQOPWP8PBg/D8b3gPMLR4Vvc5yk
8AXj9xtklTT3t0FxWx25jHE9MK4HxvXAuB4Y1wPjuuDXTp3y1EyPnVayMCZA
bahmfsDwA4YfMPyA4QcMP2D4AcMPGL7A8Lelx+RDOlUt8weGPzD8geEPDH9g
+APDHxj+wPAHhk8w/Pqf8ejFSczMLxh+wfALhl8w/ILhFwy/YPgFwy8YfsHw
DYb/6kXGVUKOY/7B8A+GfzD8g+EfDP9g+AfDPxj+wfAPhn8wcoCRZ7f7KYNE
XcfygJEHjDxg5AEjDxh5wMgDRh4w8oCRB4w8YOQBIxcY+f6RtSFILq1n+cDI
B0Y+MPKBkQ+MfGDkAyMfGPnAyAdGPjDygZEPjJxg5LUtCTuhcthYXjDygpEX
jLxg5AUjLxh5wcgLRl4w8oKRF4y8YOQFIy8YucHIf+DVIf20mkaWH4z8YOQH
Iz8Y+cHID0Z+MPKDkR+M/GDkByM/GPnByA9GfjB6AKOP0HP2RIOsifUBRh9g
9AFGH2D0AUYfYPQBRh9g9AFGH2D0AUYfYPQBRh9g9AFGH2D0AkY/rfI8U2dn
M+sHjH7A6AeMfsDoB4x+wOgHjH7A6AeMfsDoB4x+wOgHjH7A6AeMfsDoB4ye
wKv2nb1rkN2hD8IvRsmlZdT0zPtrIVfM2G3y2O1qw2U6dlw+TqupYP3gcyg+
11ZMe3OH4yB1dHb9q2J943MoPsf9AkPPnCqKEWtrmC75UHNPJjcx9XglfopS
yDE9nxl8VKOzMP0ksn8fvcLKtH8v81qHyMb0yk/ZD0X6BqZx0cmhUqWdqY9H
hFYhbmba/f6jJyjOAbR7PzgX0O594JxAu/eBcwPt3g/OEbR7XzhX0O794ZxB
u/eJcwft3i/OIXvPdusb5xLavX+cU2j3+4FzC8V9GXTsXIRSWEz587Zc0uhK
mMb/sUCsV5Qyfeu/b3/rEBmZ3p7/+JlIX840seP6MqnyIVPxyf2FCnEVU5tv
gEivqGZ6NSltrlxaw/RAiftmh6iW6Wq/ledVDqeGbipuEOlNTIeWTh2m1ZiZ
to/NWCBVckxvbfbcbpBZmB4pW39RIa5jmji+qrWz06my5Jmj1QYrU3G5NlKs
rWfaMcHrG53KxrRElfSbTN7A9HSF9YlZ0sh0i3+4n1JoZxq5LecjIefUiZXD
UzW6Jqa9AnbclKibmeJ5STmSmKNT/c6eAzCek1/r4vwC1TfYcwBm790BxYWc
JJs9B2A8J5/uDUqVS2+y+9qmWH8oUH2e3UecE/j6/GmtXSa/z84J5jv88cPI
RyAoYucG7wf4yzL/2apTXWXnBu8X+PMcUnvGIjnGzg3mBb4f18f3Y/52n69g
5MP1MI+7z1cw8uP6+H70edlW86lWc8vl+9DnCVFOUKksj/18UUWaQv4Xfx8o
rfLMMJiLKcTd+6hhj3Pu+Udl71H+xd8L/N/PShDz+fH/HKF5Ier9Bjq8RLdW
3Oyci15/XBCXtbj+fWD37vBZSr4P/L+5h5cO5B4Z6fnZ6FDlPytd5ufxRwKr
ZEEl+SZ39uaKq13m56bR96KkZdVUOqForTC51mV+Tvhdc02xrZZUFekPJQFm
l/lZ+WX8OPX7ZgrYtiZUbnKdnzvHSPfrLBzVvhesVaXUuczPKcX9BVxqHe2p
6uelnV7vMj+tSeYY4Yx6mva1KcnQ4jo/9/n9fE/SbqPGwOz6zh8aXeZnSGny
FPnRRjpYu2W+eF6Ty/x0bI44qprP8+6FF2QdzufV+PlIuZK/33u8Pu3h3fUe
W39lnZa/nw1xr95X8L2vWt6iy2rl980bx4/p+J5Xuns9V/nfpg9ea2mVZJdS
8Jhzy+SmYjqo33JOZykmt7lvf7Y1/hIdWfrtb6p+FWR9/cl73i/tD2lJCwK0
08uoffjoZOEMfo8IWxHziiCFOn7xTzWsrqKnkSajttW5Pzw1POopSK8gmfJG
hPxoBdsf2uL5txh//ud33f6OKlrQdZsF1bRImvZAlllDbgN6/tnJvyd7NcdO
VUabKOtS8SfipSbKPjg1Q9OHI7ki447Mk6NloZ599Rcs1GdQwvvKixZ6w1GZ
4FhhpdyrM9M1sVbKO6StFvHvs5hVXh56bxvFzkqaJb3ZQANF9Wsc/HtM9Eh2
VvGFnfKv5RhFiXa6nj7cW/1OM8Wv2SGVjnHuW1E1M6Vq/vn3DZg2VfvSPrFv
0ewp2ppil/0i9Ht7ucJazJ7XXyfVD+BullOLtqG3IKjMZd9Id+w+owstY/tY
vz07LBL+eQ57dddD2eEKl/3jxawJXwl/rKCMLnxR5bJ/nOxaM3pUk1v5Z0HK
MzUu+8eijZO/NWTVUJZq0glNlMll/+glfv15J9+33N+jn74357J/ZN83LhP3
5ahPZVmiI9fisn8s++pUoSzPQrnbT5pEK6wu+8cb724MUK60UoxkQ5h0sM1l
/8grCTuk8bHRwJo5PysKGlz2j9hNQ3rpCxsof6fPUHWC3WX/EI1rWuXg7+fq
SXaVbrRzLxQvHmrS8ed5an7gCinfe8uDzYWqFP73pfULfbhUAymfTJzcxs+v
77xL/67me57t9/YhQ4uRjtx+bJe0G2lCP/Fdj3fL2fnN2Hx9lXheJb0YdeCS
an4lJdnfS9eVVJNfXlqBorSaxn3n7sHx88r4+cqJ6q21VBFSvEbIz6Ntb047
qAs0046WH40SM0eBVzxf4fg5NCltvVS+t47MsVUKIT9vLDNnnlYF11OKV9YD
SUg9qdsGvalttVFQftJUeZuNgn+wfmlIb6SmuPAM1ZFGapmdY+nk50ea94i+
2vAmyipcKJHvbaaotLxoYarzvXrWva3VnOp832dKEowK/v2K9/nru7xO60TO
9/nH5pwvhfz7uWXzBKn2pXmTsvjMXPWTe3Qw6q3twp4Gij9ZTt6BzvlzfmFx
BLfuOrWfuxwpX26kj2f/2zcwu4DS3S4tkd8vI1HAYm5pm3M+7Wv5YVjIs59I
1nuZn/bKQ3KUld/g1v1Czxd/NFY70bkfFJTlGUS+zRS5fLhGleLMhf0A/mz/
2virbuB9Cjk1ZIdwXQnbD+CP/T7+/+uLnv1827D6LjkimhySUaUkuZI/Xs8/
76bY1AMO/v05p/zeAEGgk72+8t7vyC2h1cPsdbK8Esr8fsJW9ZMrNDfzm3zV
pnLqGOU+OPyl83bdt7FTtKKUht66KFOuLKW1JboY6WAj6RN252h8jNTn2Wm3
Gc80dLznxH2G6koa+c2fx7xfOo8jNu26pygopy9GLh2hLyynO2PlU9QJD2n0
3fFfOxIf0saycUd1o6voQaKgXeRbReqrncmOFdUUtDcmUxNbTbFeMbmq4Bry
qi+qlYTUUPCqoiYRv/82TZs8SO9dS6L8SSO1rbV0Xa2ZI2+rpZbBmg+lN02U
ZvNIcvB7cEGcx05DupnWBsVnq46Y6fC1NZcVX3A0e5/RKkrkaL23sb1znoVG
NAYP0YZbKGxN8Bj1O3X0dPqpcOmYOhp142S0mN9/i/b3VxkEVnru0z9FV2Kl
DPuGXEWplX6P33BdlllP/5Gamzr5Pfh4gem5OdlGC74LG6neaqOkoWHjldEN
5Nec/aF4aQNFrM1eLuT3X7cZQ3bqAhtpnN7ngKaPnYwHky/LPO3kPiz5lsRs
p6wWe7uZ34Mr1tkF+gtNpAqNGKO82ET/A6HeXmk=
    "], {{
      {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1], EdgeForm[None], 
       GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJxNmHecV9URxe9dii4gK0WK1EVUiiViUMounRVYQIQliAY0FGm6BGFRpISl
CAgq1d4rIcaSGIMlhhBjbzFGjbGXGGs0FsgKknM+9/s+b/84O++WmTkz976Z
99viKXPHVhaEED6JIdSR/K/wlXCzcIhQT6gvdBPaCx2EvcIe4SPhG+Fr4Vbh
eeEB4RahMXoHCd8J3wp3CC2ERth+Udgp3C78IOwXPhX2CQ2EZkJznhsKhyGt
/71QI3ws/A+enm+JbftvhSwSivHt9XZCU+y3FZowbl/LZ0f8WacTerZ3mjBI
GCwcwVxroQu2bPsl4SFhu9BGOBQfLwuPCDuEo4TDWT8aaf2u2Mhy3Rwe3Rmb
11DhZKGX8GP07f8Y1h3rcfA2x1eER4VfCccz11noiV5X7HXDTwn7fiScxHo3
/HXHT2/ksUIfpH32RVq/h3AksZ6INNfXhF3CXUJ/4QT2noLdPuS3J/4HEafH
Q5gz3zI4WWcYen3hUkz8w5lzTKPxNwD0gNdApH2MYG+p8LqwW/i1MEroh/6p
6FvvWmGtsE4Yw5z5XiOsYW0idsuFKmGaMF2o4Cwdx0+I33GMZ87jN4THhLuF
SfBwHFOEscI44afCSNYms26OZyHN62yk7+/PkNafgE/n6XSkuU7FdgV8K+B1
Hv7MZSY6jm8GtjyexdwZwrnCmehUojcZ20OIfy5z5juf2KaSo/Hk5hykffyc
vY5pHtIxnY+0/gJsTCPOwfhbiF3b2ygsFZYJvyA2c7yAdce0iDgd0zvCk8K9
wkXMzRaWo+c4VsDPvHz+VfisZt1rK1k331VIx70aae4XI6vgOIdcLkOa61ps
m+865IXCZuIxr/eEp4X7hA3wXixcLizB9qZaOdiCXjW2ZhD/VuYc31VwNcdL
2Oc8rUfaxzb2OtYrkI71SqT1r8aG8/SW8Lhwj3Ad8dj2BDj4PNxz3FvcV1w7
XU9cX27At+P7QHhW+I1wI3OXCneG1KMKQ+ob7hnuJe5L7jmuz+5j7ifuN+4n
7iON8NeY8e08u6/9MqSe5vp+W0g9s25I/fMjbLinuZ81xa7reTPGrvNtQ96r
st5l6b6znfki+HlvG/Tbo5/1j6Ow3YE167nOdwp572xOPC0YGx3R8Vwx45bo
ebwDO62Jsw25N5cu+HbddM0biG/3iq7ouQ8ch4/urNmua3TWr45hLeslrudZ
r+oc8v7Ro5a/E/Hfk7H9ncTYPk5mbLu9GHeCS+a7N2tZ37LvcfoO2io8g63B
xOZe5b7gnuh+5T7QD04DGGe9ZADcS7DbDn5d4DQEu+4rrrWl8BvKmnm59pfB
aRhj2xrOuATd4cRQhv71Ib036+HnnjASTu4Fo7HlWuramtXNecTpvjAGHq7R
We0cy5p9T0Df/l3nTw95v5kc8l7ifXezpxwe7gNnwmsS49HoTSKGCnyXwbEi
5D3srJD3sLOx7To8k72u7Vm/spwa8v4xHV7e65pyD7pncG5jiMP8XGPn4M81
ey45mILvcdieQgznoeMYKhlPRq8y5L3KtvoT8yj053MGrpsbOLusx8zGzwL2
OQbX/CpicO3N+tUFrJm36/oSeLiuLwt5v1kZ8l6yFH7L2We+1YzNdQXjeeit
IH+L8D0LnotC3sNWhbyHrYaH+81l6Ph+Zv3Kcm3I+4fXsl61GH4b0Z9G3Avg
vok1c3K/2cr6GnwvxPYaYtuMjmPbwrgavS0h71XbQt6rriCGqxjb7tWMsx7j
/vIUvBfj09+F1xCP38msp13H2oGQ+kX2+8p9zT3FNdY9xfX+d8z791mMqX+5
h3weUr9xr/l9SD3Juo15LiKXN+Dbv0X+FlIN3wBnrz0YUo/x/j+FVI9ci0bw
7Hf2DyHVYtfh43l2DXRt38H8n0N6z73fd+Fe8vAEz77DrfBlDg+j6+eJ6Pod
/EtI76SfZ/Fs3S9C6t+F8L+PnN/IGdyEnyfx6xzsJK5n2OP9Hwq/Zf9z5Ne5
/VK4n+e/kpcicvYwPF3/d5Gf7eS0NXE8Qh7+Tr46k6dHydWrwh959nk+wHm9
wPn5+R+cQSm5300+HevT8Hc+Hicnb5Iv5+ptcu35d0O6i87D+8Rv3RJiKCUf
z5IH+3+RnP2TsyznTB7D/v3kq4CcPV8rjhfgv5PcFRHrXcTrXL5EzpzLl8mZ
8/QKedqN73Jy9iq6u8hLKTZfw66/r/x/igjfN+DrfLxVK09vkxPn5p2Q35F3
yY/z9F6tHL9PrnwOr3MOzt8HtXL2IXn7d0jfeObwSUjfjf7/yV7WnKsGzPv7
67OQvjMPDul/F9bx+9uE+abk+01y/gM6hfj5Cl/1GNdn/Qv2mNd/Qrrf9v8l
HG5Gtx4+P8eva6Lrk2uZ76/vru+tv4G/IZaPGdfB5l587+G5AO7fEtcdPDch
Zn9j+3uzIKb6tQ8O32N7ruZmxhSHgzuAjtdqsOH5GjjXoFsXP9/hN+OU1dE9
cHNO7bMZeTiA7mfk91P8HSA/NfjL5vzs/a675l/I89ec7T7m/Xyu5IyY4rPd
/fh37Ps5o/34zZ4LyMlsyXNiyrHtFeDLv238fenvzY6a2xbTHVwpOSimGtRX
8qKY3sv5krNiuidV+jOH/+ltlmwX030fKLkipveykWR1TO/uwZJLY3qnF0pW
xpTb3pIXxlQD6kou4vm0mOzaZlvJTTG9W4dLXh7T+zdU8uKY3vVDJVfHdNe6
Sy6PqaY3jikexzJYclVMdaK95JaY3kvHV48Yu0guian/2ecY/NrnaPxabyy6
l0qOjOndst5B6DqOPsTiuEuJ3TmuIM/rJUfEVJ/WSQ6LqSbVkexEnhvqzwRi
OVqyhD7tWIcQr3M8gDz3k1wWU91dI1kWU81zPhpiZ5TkZTHVs1MlN8ZUw8ol
N8RU82yjEDsNJDtzXq1iitnxdpPsH9N7ab2W6NaXPIK7YT+t8eV7cQh34xLJ
w2KqhT6TIs7F59mEM10r2Sym2jxP8vyYatB0yWLhXyHdnbrk2T6r8Hus5MSY
fkMfKdkrpnfUsbYh3uYx5d05d3wLibFFTGfjc2kaUx6dQ/9u6xBTzR4eUwzm
77Oaz3mdEhNvc/ZdGM99+D/TwHpN
          "]], PolygonBox[CompressedData["
1:eJwllMVuFVAURc+hQArB3WWAE9zdWqwtUNyhuLVASBgw4FsIBMIYd7fi7i4h
BAshjFk7e7DSlZeX3nP32fd1rqgqr6wREcVQE3pkxBg4ht+FpXgzOI8/g1H4
LviL74Ud+HA4jP+GIfhsuIl/hC14d9iH/4AleBM4iz+BuXgX2BP+ziK8AZzE
H8BKvB1cx9/BMrwFXMRfwAi8K9wOz7AcbwWX8VewAK8Hx/F7sBhvBKfxR1CB
d4Bq/AMMw3fCn/BMK/A2cBV/A2Px3fAP3w/z8EI4it+BVXgnuIV/gtX4GviC
H4LydEbKRmcOxmfBDfw91E9noLufgJbpO+kul6AEHwivw3dsnZ5Rs12BcXhv
uB++8wS8DzwMZ1onnYmy0I4bpjNRFqdge3pmzfoLGqd3pl2dgTJ8ELwNZ7IJ
3wzf8AJom96ZdnUN2qczVra6Y8d0RspGHdmAb4SvOH9iZjoTZaGdr8e3wk/8
O8zBh8LncMbr0h1TtzRDJV6V/m5huBsj4Ui4I9PwAfAy3KEpeH94Hu74fLxb
utsHwmf3SrtmWIj3TO/+IPTFi9LZPA5nMyN9d2Wks6and6cztbvS9K60w374
ZDiHPw3PMhUuhGdqnu68uq7PaqU7p67pTupi3fRd1clJ6Zk0izpenD5D/1tv
rmn6Test60x1Y2J69+rI6PRvgN6+Or02/cb1trUjZb8t3Q3tQL8FtdO712+C
3k5Bepd6Q+ri+HSX1cn/piCTNQ==
          "]]}]}, {}, {}, {}, {}}, 
     {RGBColor[0, 0, 1], AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwt0ldbDgAAhuGvyyxF02ioNEhJSdJQhNIQCtFCSlKKkoyMaEi0EyFOHftP
fke3i4P7+H0PnuT2/vq+oEAg8JtfhBNLGrmUUkszPYwwwyxzzLPAIkss84kV
PvOFVb7yje+s8YOf/3f/kOxEJx0M8JACLrNMErdJZIkGjlHPIntpJ4EFLpHP
ReaJ5xZxzHGBo9QxSyw32cMM58mjlo/s5ga7+EANR6hmmp20EcN7qsjlHFNE
00oU76gkhwomiaSFCCY4y2HOME44zexgjNNkU85bttNEGG84RRYnGSWUawTz
klIyOMELtnKVTTyjmHSKeMpGnlDII+5wj24G2ccGhjlOGlfYwggHuE8XvfQw
RAqbeU4J+2lkG6/J5AF36edvrI9JJYRXlHGQ6xwK+tfxOqkfPwo=
       "]]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, Automatic},
  BaseStyle->{FontFamily -> "Arial"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox["uE", TraditionalForm], None}, {
     FormBox["uC", TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 1}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.834514278930832*^9, 3.834514306303125*^9}, 
   3.8345143962456903`*^9, {3.842274520323626*^9, 3.842274542028345*^9}, 
   3.842274710830452*^9},
 Background->None,
 CellLabel->"Out[33]=",ExpressionUUID->"313a5eee-569e-40b1-a654-e775860f404d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"eigen2", " ", "=", " ", 
   RowBox[{"eigen", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"r2", " ", "=", " ", 
  RowBox[{"RegionPlot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen2", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "0"}], " ", "&&", 
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen2", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "<", "0"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"uC", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"uE", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"uC", ",", "uE"}], "}"}]}], ",", " ", 
    RowBox[{"PlotStyle", "\[Rule]", "Orange"}], ",", 
    RowBox[{"BoundaryStyle", "\[Rule]", "Orange"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.834514372808457*^9, 3.834514390636292*^9}, {
  3.842274544405828*^9, 3.8422745456513*^9}, {3.84227463328498*^9, 
  3.8422746401381397`*^9}, {3.842274687781416*^9, 3.842274692152658*^9}},
 CellLabel->"In[34]:=",ExpressionUUID->"c7e86c4e-3b40-43e0-b4ff-f51cde07196a"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJx1VX9M1GUYP0kKKDAsY4N+LIR2DgRSj5IpDwj+GPLj+GFrCrfFilUQ1xhZ
sLXZ0ENJGtoJRgJiYALqKTALPAMUDzsCDuX4Hr+R+maJCZweASElW5+X7d2+
f9yefd573uf5PJ/neZ/vq8nquPfsZDLZ6OPfonWoqXEpD00Lhm1t94o9aK0m
4OhrbupcayPDSQWdu0VFK8PPkeuCxmpkWJyuUmoFE8Mh0Ts9RUUPw99W3nug
0AoMT8/nX9dY+xmOTfAvEpTDJMWLt+AHDH7A4AcMfsDgBwx+wOAnpQvPg7fg
Aww+wOADDD7A4AMMPlJ94XXgefAWfIDBBxh8gMFHag543fm6+by8RX5g5AdG
fqm54/vM68zXyeflLfJLzTU/R3wfeR35uvg8vOXfDT+n/JzwfeJ1498dP+f8
nPHvsuKMyk9XPki6weBVedZiho8/b24XFfW0I/ZKoiqsj+w2rj2kFZpZHGD8
/+MdtW+gto3uz3p86SyaKb2wtktQdjKewPj/6KnsBr3mpuQ5sFQeYCnewKiT
94NOiAOdkZefG6m6cE9KJ+ThdQav7AVDcljWELm8NHbud0UFvfK99xGrez9d
+O3RlF5zlVpj9jeV63vpcdNelMlMzOIcfriHOIgrry5qT5WP0P6FjJJ8ax51
xE3PuhsHqMq9IcSi/Iky/kmQG1Mt1Hz3doqu/Bdyq6h7K8vZzCzO4Yd7iIO4
yIP5AubPwWvzlfVzWuEC4+20qeyas3id1aWPXHlArutmFufwwz3EQVy2N/73
R3/gj37AH/PGW7zXD2bG/lKqbjE9+DkHxjykfB1yTBX2M9OLf0fAmI8HqXtL
ArU/MD2l/IDhh3tSPIDBA7xQF+pEXahzrjZpW1bUMNOpdI8+Uz7Rz3QMX+7x
nVAgMJ2/2mX5V9FjYnkRv/DtiCDdbTOdj/JpkPu3sjqwvyZ140/LQvrIsaXI
eEyoZu8B+yryqfwhZekgnc4xpMtjhill27DzCrGaIny9S4TJfnpyYt4xd10r
TfbktGuOChT0WqhcNJuYxTn8cA9xEPdM9IJM9vEIbW9b7Z1nPUKRfyf661wH
qfps/Atr1E2MN+rYe9GrSm8YYO/cWnY5SVVnoavrdd1CQAfTAbq4784RUt+3
sL1RuMP9sPMuM7O4hzh/nM3q1qzrYX0FD/B6d/6jROO4kc0JeKOOyc3dPi5i
I9s3eHenF+H8CNMZupscDNnO3iNsL0F39AFxEBd9Qx/BA7zQd8wB6kBd0AG6
QEfoCh7gFbdId2GEYgp8W1ym3WlGsNnLTg7SI/vLz4bPNtKJfQkbdKF99OHo
WLRoM9KWNXXJql/NdO/l+gjtGz3M4hx+uIc4iIs8yAvdoKNnks1V/GRo6Tsz
5fONEDBExW86fZqbcYlu2CeEaN362fuLKM5r09zso4YsJ7vc2huUdr42Payh
d+k7t+VPmyKzl7Tj273EVbfYOfwQB3GRFzwwx56xqdplshqaO7X2uNA8RMEx
nV2i7SLNbN0ToPtigJ5ob1OJtha2P+KnD90RTM3sXdS/frImbGMrle48bNAk
99MzB9L8AjMNdOJ+g0q12kJqeVSlXN3F9klT/Ao/j4hO9s4K2lZmhfWZKPzh
3YeK5QJtmvqsUT/dzc7hh3P4oW7oAB7gBR7gBR2hK+qGDqgbOqAv6FNcyfwy
2egwJV1yfOegtZx9n4LavBwq9cVsr2HPYe9hD2IvYk9CF+ik/XxsVrFhad6g
O/qgL6Fz8oml+QMP8Crbl9lhHF+aR1ic4x7iIA/y/gfzqUZS
    "], {{
      {RGBColor[1, 0.5, 0], AbsoluteThickness[1.6], Opacity[1], EdgeForm[
       None], GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJwNw1dTEgAAAGD3ZCMb2bJlylKG7D1VphNnuRUEQRDEh67rUcuyeXVd156X
mS0rW5alPXief6B/kd/dR45O+qIlRUVF/06WngQUl4BKy4AlpeCycmhFJaS8
AlZZ5TPU4KAwFBCEBoExYAgWAkXUAuBV1RAsJpuBExFIEhJFRqEpaAwVg6Vh
cWYTHgkAMghEJpHEIpHZZAqHQs3naLwGOqEOwaczBAymkMmyWbmSRp6Ux5fx
BdoWpZjNSal0coGwEBOrJFJlk0QtlfkdioRcoxCK6PWEoEtdi6grAQJKQcAy
MKgcAq6EQYsBtWGPtgIKEbHYACQCiEKC0CgwBm2AV1fBYU6/whHQLGaNCCKh
jlCPJBHRFDKKTMJQKSYaLsehEpkMEotJZrMoHHY9vQFLo2ZwsKW8hc7nMQR8
plDAEgnZYpFVzGmUNBEYdL5MKpDLhAp5TCFaLtiblC0SlbKB1yhVqxJqmUMj
T2mVLp1qMuj26Fu5Vk6LtrW7wzoV9na2GfAmXJsR3tVuarfQ+s34GkM1FIcd
sHGjvpqegCPuVKhdqt6QS+tpVTjkHXZx0q2e8DsNbXBTO83aIZ7u8ae92r6I
hyeVzPQGH68NrD8dbBs1dozbbd1WX6fBE9F7w1p/jyOZ0Ji7TLN9YVdIl0iq
55Oa5LxbHBNNzYbjseYru865ePPIKf3kTDAWV6hSyolpvzwhi88528cswV5X
ZEQfGtbBMzBajorL4I1ZeCaLzy7iOTmuvSDOL3Fzea4lTyssN8cKzZd2bNvn
zTsXbemULpXWhvs8gSHN6WjNq03LQlp3/CDyesv+csOYXvAe3Q8d3gusfRlw
B9VP1gfPHUfefxg9exS6tufee+veuG389Hls/533z++pM4eBy8+Gbh34fv2c
/PF94vqLkdWP/Q9Xohe2zbtvnDf2vQc3fY9W+68+H175G926a//6bXzzjuU/
dnTQtQ==
          "]], 
         PolygonBox[{{150, 143, 93, 71}, {201, 187, 75, 109}, {147, 140, 91, 
          97}, {199, 154, 74, 100}, {171, 156, 102, 87}, {195, 151, 71, 99}, {
          189, 203, 110, 77}, {180, 141, 66, 103}, {168, 223, 130, 72}, {213, 
          145, 95, 113}, {214, 146, 96, 114}, {215, 147, 97, 115}, {153, 144, 
          94, 74}, {163, 142, 92, 82}, {178, 137, 58, 101}, {173, 160, 104, 
          91}, {166, 222, 127, 69}, {170, 224, 132, 75}, {179, 139, 63, 
          102}, {182, 174, 98, 106}, {172, 158, 103, 89}, {212, 164, 82, 
          98}, {145, 136, 87, 95}, {191, 209, 112, 81}, {221, 161, 81, 124}, {
          222, 165, 28, 127}, {217, 157, 77, 118}, {219, 159, 79, 121}, {205, 
          214, 114, 111}, {133, 192, 107, 83}, {223, 167, 44, 130}, {135, 200,
           109, 85}, {134, 196, 108, 84}, {137, 217, 118, 58}, {193, 183, 69, 
          107}, {190, 206, 111, 79}, {202, 213, 113, 110}, {181, 162, 105, 
          104}, {142, 133, 83, 92}, {141, 221, 124, 66}, {139, 219, 121, 
          63}, {146, 138, 89, 96}, {208, 215, 115, 112}, {165, 175, 99, 28}, {
          169, 177, 101, 58}, {143, 134, 84, 93}, {167, 176, 100, 44}, {144, 
          135, 85, 94}, {197, 185, 72, 108}, {224, 169, 58, 
          132}}]}]}, {}, {}, {}, {}}, 
     {RGBColor[1, 0.5, 0], AbsoluteThickness[1.6], 
      LineBox[{2, 1, 11, 20, 29, 37, 45, 52, 59, 64, 67, 105, 162, 181, 160, 
       173, 220, 140, 147, 215, 208, 209, 191, 210, 161, 221, 141, 180, 158, 
       172, 218, 138, 146, 214, 205, 206, 190, 207, 159, 219, 139, 179, 156, 
       171, 216, 136, 145, 213, 202, 203, 189, 204, 157, 217, 137, 178, 155, 
       177, 169, 224, 170, 188, 187, 201, 200, 135, 144, 153, 198, 154, 199, 
       152, 176, 167, 223, 168, 186, 185, 197, 196, 134, 143, 150, 194, 151, 
       195, 149, 175, 165, 222, 166, 184, 183, 193, 192, 133, 142, 163, 211, 
       164, 212, 148, 174, 182, 106, 10, 9, 8, 7, 6, 5, 4, 3, 2}]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, Automatic},
  BaseStyle->{FontFamily -> "Arial"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox["uE", TraditionalForm], None}, {
     FormBox["uC", TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 1}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.834514393223419*^9, {3.8422745231164713`*^9, 3.842274546137476*^9}, 
   3.8422746410222263`*^9, {3.842274692872795*^9, 3.842274713104623*^9}},
 Background->None,
 CellLabel->"Out[35]=",ExpressionUUID->"b168c029-3e40-4119-bb00-4ee6276a91ae"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"eigen3", " ", "=", " ", 
   RowBox[{"eigen", "[", 
    RowBox[{"[", "3", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"r3", " ", "=", " ", 
  RowBox[{"RegionPlot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen3", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "0"}], " ", "&&", 
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen3", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "<", "0"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"uC", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"uE", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"uC", ",", "uE"}], "}"}]}], ",", " ", 
    RowBox[{"PlotStyle", "\[Rule]", "Green"}], ",", 
    RowBox[{"BoundaryStyle", "\[Rule]", "Green"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.834514489270554*^9, 3.834514494944057*^9}, {
  3.842274549048967*^9, 3.842274550913122*^9}, {3.84227471512286*^9, 
  3.84227472092936*^9}},
 CellLabel->"In[36]:=",ExpressionUUID->"068e9cfd-920e-4a4b-8c13-9ff540e8dd83"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJxdmHlUU9cWxqNFLWpoq61iAOtUGy0qKqHqq24kiE/lYfApVdTYpR2oQVCe
Q4P2laUSZ2sxqKWoxIGnqIQqrRUjTqgYHKsQEuYMEAJhCDhRS95lLb8TFv7j
+n3se863z93nnH0zdHnsvK+683i8E914vI7/ixZE+St1lfT2mTMequnR07Z2
m3xYozAynnDOvadFZGZcudAQw7dYGO91O6MTKWsYf5a1MVAqrmVsWxx6WuGw
MT7Uy6efWlXPOOSCPV4naWDcLM01OZ2NjMPnjzuok5RTXsGI8O2ODIJfMPyC
4RcMv2D4BcMvGH7B8AuGXzD8guE37ObA2G2OHOYXDL9g+AXDLxh+wfALhl8w
/ILhFwy/4Oev99xSOAy0dN+DSIsoj/kHwz8Y/sHwD4Z/MPyD4R8M/2D4B8M/
GP7B8N+f3mtXOLTMPxj+wfAPhn8w/IPhHwz/YPgHwz8Y/sHwD/7lZH2LSKkj
y/PTEqXuEcsHjHzAyAeMfMDIB4x8wMgHjHzAyAeMfMDIB4x8wMgnMGzOMIvo
KcsHjHzAyAeMfMDIB4x8wMgHjHzAyAeMfMDIB4x8uvpHXvAPhn8w/IPhHwz/
YPgHwz8Y/sHwD4Z/MPx39Ys84BcMv2D4BcMvGH7B8AuGXzD8guEXDL9d/cE3
/IHhDwx/YPgDwx8Y/sDwB4Y/MPx19QOf8AOGHzD8gOEHDD9g+AHDDxh+us4P
X5gfjPnBmB+M+cGYH4z5wZi/63zwgfnAmA+M+cCYD4z5wJiv6/iYF+ODMT4Y
44MxPhjjdx0P82A8MMYDYzwwxuv6PMbF82A8D8bzXeMxDuLBiO/6dzyHv3dl
xLU7O/5VkLp02gc7HSmM0X+dOCUdq1aVUvfJY3Yoddfon+FXlkjFeoo5cP6h
TvKAccMrr118SyEhHjrGg37uzKMoobqK6Ys/Vz2WSE1Md38rboqcb2H6xcyg
4ypNNdO/iuzfVyuzMr1/T/Nah8DG9Ou/ZpcJtHVMj12aGCKW25nu7R6hlgkb
mY71gA6f6I+QB+5r5In7AecsdMRDx7kGHeNDx7kCHesCHfscOtYJOvYhdKwb
dOwb6FhH6Kh76FhX6Khj6Fhn6KhT6Fh36Kw/eqPjPUBHXULHe4GOOoWO98T6
8Llbr6o0RTTwxIUIOb+Q4v6aL9TKipn+4f8++skhMDD9/rznrwTaEqbHt99e
LpaXMV2YcbBAJqxgum20v0Arq2T6jYTUMKm4iumHity2OARGpq/2XXVR4TAy
PWRzYZ1Aa2L64OKpQ9QqM9Nbx6TPF8stTL+3xWOHTlLN9GP69Vdkwhqmx4+r
aHY6a5guSZw5UqmzMl1Yoo4UqmuZ3u7n+aNGYWN6kSLhpkRax/SzpdYXZlE9
07dOCPeV8+1Mj9x+6Qu+xc708eVDk1WaBqb39N95V6RsZDrq9ND7hQUWUbbr
nn3DWNcW2frDAcqLrO7/qIn1DVDeYfsCjPf09f7AZKn4LttHScfiL2kUf7J9
BsZ7//al0S6RPmHj43sB7znL/HezRnGDcdfzDf0s6oY7DL15PBfj/EP+qEvk
j+9B1I2Hj/FctegEY5y34KlXJrYpdVnMT+/Pjt7kW26x+TSh/RKF6sfsvMK8
mA/+sd5grDfyxXrAP9YbjPVGvjhX4BfnDhj7APlhvtOCS4HFklw23jVb1ddq
1T0Wv7U97vAex07mH+ct1gPPwy/OX6wPxoM/3FfYhxj/synDzqs6/X5QaUyo
lHT63lXsqvDgdfr+Crlz6oSGi49Z0aTJas6gutheT2Qd5/f66+vUXPw+z6+7
eXHxhm+HS+VcfPvvE5J1q7n7LXRl1Fu8JDq27Kebin7cfRn20Tfb4q66+uol
2fvk3PM3bsX24nWut7ebmkXZxbTKzbNNMeE+q+dHpakyKTf+NzEPROpO9Zai
3XpBU11IQaMuLJeaCql16MhE/gw9pSbM91dP15NEfidCeryUXuqe9eClldL8
jtfFq6R5HWXRztXJ+z3+dnLrtVCc+lSSWUVZVwu/Ei4zUc/G6KnypSaSytIf
SDwslJ0yNV3V10J9B274VH6lmpaHeLyjvVxNOTdmpqmirfSuo3yDY6WVomI8
3bVeNso9rK4UcOf7AEHtGgd3jkfPSpglvltHebcuGQTxdhI8k5yXfWenuDU7
xeJRjXQ7baiX8mNXX+3512WhvimDfHze00s67Y+9e8NnyTvdx6iHl5Emg7o5
hyZ8mrVB2Ol+xv619n4x0Yvj4gqPdJ3ZdV/jnh7tP22quhMfWDh7irqqkK1/
sJvXcd0+7vvqF3uJzFrI4hwhucHKgzoWh/ma1HV9eIF69r6PLNasFTYaKM2x
95wmRM/iwo4sG2B5ZmBxyCe0154yyZFSth5t55eGyP9VTq9n+f3AP1nK4k49
41lF88tZHPZPesefX1ew/Tk60dnHUlhJGR3bo5urf9g88vESsb6SxWH/dC/5
JlB+ror1DcV+j9byE420cNPkn3RZVSzO70/VLdl2I4tDX5GlmHRatcTE+gpF
aVqZyN9MPYW925xcfSGu/Pu4scpPzSwOfYd0gns/bR8L6zv8t68JkZq4Onxi
WC58x9Wf7B4lPqipdsWhL+lbro935FSzvsQ4MUitSKqh5T+cKZDkVrO4KYX9
eZbkGhaHviVnR4ZJwNUz+pZ9Ff081dNr6d1PNvnLV1lZnDXBHMWfUev6HehN
XxMl2hgqHmRjfc20XaYEXRO3L4pCD6u8Xf3PAd/fHotaXXHoewZUzflNll/H
+p76gOxa59F6it7s01NbUMfigosTp0iP17M49EV5u70HKzfYWV+UYtw6Tzi3
gQRjG2Ic3P5DnGNLxHHFvAYWh3tr9SS7QjOykVriuFXm6sSvn/Ch+ycl0372
Kv6Hkqs7+Yvxk1u4/Tc1L2ClmKtj4aLBJk3H71vrF3hbknXU9HRLgSJJR8fu
P7eLWg002/ejw7omA70eceiqYl45pW+5HSOcW06+uan5suJKSrBPTNMUVXLn
6Krxym1GGvuzm7uFq7ftH0xL0QSYqTS4cA2fq4+A6x5vWbj3vLPppEFktpA5
ukLG597fpNT1Yun+GkryzHoqCq6l6pkzzyqCaikwL2GqtMVGypaBH6ibbdQQ
G56uOFZPQUet3+vS6inVa9g76vAGapp9qdrJrQ/yxf3xpfnS93zu/kK/0rbo
izHq8WXk0Jfcsaz7nXrv8TyrEbj6k7TuVxdLn+jpy9n/GR2QnU+Zog0GGdd/
oB8JPuOzk7+uiGz/3vSHZsATdr8I/BdZlrVkkKTPcl/19TI60HR0SPCrXyks
88c8xeYSyvzFb5vyxXV239UPHnJCyt0XuI/iMkrIKyCPWi9ci5SuMNDFBYUR
lnW3yRHR4BCNKCbBq9/u61Y/ZP3Dwk0nbSru/sD91bTFT6zm7peUJR/u4Pfg
+qdF58KULx7TktTcpfzkRsoqWCCS7m+kfH2uTjDa1d+g74hcMVSlSGqkUz3G
H9BVllPfV2e7z3ilYve3W9+iekNTiut7c++Cy5L2Rjrv1tJsTnaNg/oa/uPf
J7w61Vf7CLdB4Vx9maKTDzm481t0PW+cljuvwXNKHr/PC3hKg+9dkchXFdPt
0fVOwcpi+m74smHaghIatnnPY1l+CT2N57UKRlfQJv3Y45qRFeRZ+8goCq6i
aM+oHEVQFd1WquZIW4wkyJs0XN1spLWBcdmKY2bKj3XfrUsz07D6IB91eDWt
9zK0OudW06OD/RU6npVG3MlYKuT67f+KzQ1Orq/+M27jbUlmLfk2Zn8uXFZH
CYNDx8mX1pEhJfGaxMNOY7Xeh1R97aQIiRglv9JApevsPO3lBlo9xF4jyS0i
zx+8Djpyiki7Ye8llbeB1hZposSDDDTy4bhdjvgyejBGOkW5oYwC90dlqqIr
SXnDmehYWUkN0yYP1HoZKSjmUYOA+95ItbknOLjviqZBqs/Fd000+4DBKoi3
0JFba67JvrPQy+lnwsWjaih0TdAo5cc1lG7fmCMrtlKbd/8kTZGV5v8cOly5
zUan8k1t5kQbdZ/hs1sTUE8Ra7NX8Lnvgawme6uZ6/vdhiTeE5nt9H+h3dwe

    "], {{
      {RGBColor[0, 1, 0], AbsoluteThickness[1.6], Opacity[1], EdgeForm[None], 
       GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJxNmHucl0UVxs8sIggLXlBA1BTUFJCrgLC6CKLckZsIosg1TUIR3ZQVVjBA
EFbRBV25uobl3exCRkpERkaWmpqZmppaCWUUFabU2nk83/cz+8f5PDPzzpzn
zJl5z/P+fu1nzB13TYmZ1SezRo5NHQ9xu9/bjRz1rIVbS4u25hwOHuJ2BNjY
7UjwULejwCZuxzNPz08A9fwYt8Pcmrn9ya3UrY1bW9riPbYBfztQ/MeB8tsK
nqZuR4Py+wV49Px01iiWk5gn/kN9j/WOW9z+B//Jbqc0iOVUULF8EVQs3RvM
PREe+e0Ij/baCVQsnUHNPQNULHvcOrh1detGW35Pg0dx94BHsfQEFcuZoOY2
8b0ktwe83Ysx7bs3qLj6gIrrLFBxWQr+s936MqYYz2kQ1yB8yUd/9q+4zgUV
10BiEv8ktxFuI93OY0zry/Gn9efjT7FcAIr/r25lbkPchhGXYhkOluO7nFhG
goplFDjA7UJQcQ3Fp3yNZkxxjWVvimUMY+qPY0xxTYRbnN8grvFuE/CrGC8G
Ne8i5gyFvyd5uYRYFeNKt6+5LXG7lFgV12WgYpkCKsbLQcX1TTimuk2jrViu
Yo3mPsiz2W4z2YPOZBaoWL4ETna7AlQsV4KK5cug/E6HR3udAcrvV+DS84fo
V7hdzTzxXwOKfy4o/mtB8c8DxX8dKP4FPNPYHHjk96twaewGUJw3guKcD4qz
EhTnTaD8Xg+P8rcQHo09DEeV2820xbOUmDR3EWPiXAyK8xawknOuhPMRfC2z
OPub8HUr4/J3O/3lbivwK3+3gfK3ivVLOTfdKd2hR1lbzfpFrL+jgb/VoPzd
Ceo+PgZ3jdsa2lq3FtS6u0Gtuwe8y+1x1tTSX0mM9zIuH+tZr3XrGFP/CeZt
cHvS7T63OreNjGvut+hv5rnam7B1+C7xOrbfogbq/dO7OthCu6Qd0pJCw9Qv
NEz9QsPULzRMfWlLK/qFxhR91fMT6J9IX9olvWhjWcPaWtawYy1rWDv4jqdf
aKT60rBj4PuzhUa2pn9SAz7V6k7wSSNOtqxbp1jWrVPhO42++E6nL76O9Atd
Uc2SNkoXO1jWLfHttdCtLjxTPe8Kd3f64u5Bv9CtHpZ1sSPcveArtKqXZa3q
zXrV9nLLutWHOPrSl+5Jw8qIQxpxtmV96k9M5TxTLAN4pjWq00PgVv0fCPcg
+oU+SQ+KGq36+JGFVg1mDwPxK+5h+C10a5hl3RpuWatGWtaqUfgYTV/cY+gr
lrH0Fcc4+oX+ye8Lblst7r24VasLLRKOt6xVE+DQPiYT00SeyZ90YRIxTaY/
gjHN+7fbH92+bVmrlBfVGNUK1Y6X3J6yrFfTLDRL3NPpF1oyHb8z6RdaVfRV
z+dY1pKrLWvVFZa16kpiuYr+y27bLOuh9GI2fHPoF1o4y7IWFv258BVacqNl
rbrWslbNg+96+q+6PW1ZD6UdFZa1qtAS1eebLWvhXMtaNd+yVlXCsYC+OBbS
f81tu2W9kr8qy1rYUJ8WWdanxZb1qeiv4AzFsYRn4lhK/3W3HZZ1TFqzzLI+
rYBjOc+0fqVljZE2SDPecNtpWa90X6ot65PuUPHdM8WyVq0mjlX4e8vtWcva
Ja2psaxbayxrlbTmbbddlvVK+lLL3HX0NW89/XfdnrPQqEL/1jJPGiTNes9t
t4VGac0mnsmkVdIs/a7Sbyrpwwduz1tonX5zfd0+L1/2mds+i98jjX3ggMVv
MdU7vfOqAd+3eL/1Lt/pcwa4PePtH1i8Z8U34YO0f0hbd1738CHaz9DWndSd
eZj2z8hNLXt8gvaPOfNlnNejtH9Cu5r8P0b7p7SLb4nHaf8IruKb6hHaP4dL
+XyRfWqPvyCvhe4/SfuXtOssfsd9x0If9Tvouxba9Cu371nUZ3FuZ48b8bmB
nD0F16/J41Ryto32K+RxNjl7mvZvyGMFZ7IVrt+yzyr2t4P278ij8vYOua4l
fzsZf5OcVpO/Z2n/npzWkKfnWKu87mL8D+RxA3t6kX29Tx43k7PnaW8lR4MZ
e59x5f4F9qQ87cancvMSOdtGvoo79jI5Uz5eJU/byVEF+XiNc9hBjqrIx+vk
aSc5WkaOXyHPyscb5Ek5eKtBbt4kP7vIUQ17fY/YleO3yZNy8y75202+NpDL
d8inNE13qw78gJw15r8Rvb+fWLyb+iY7LAXq++tDi2+11sxRu7nF3fyQ8U8t
3nF9P+lu/o38/53x9szZS1t3+S8W9/mgxfdt8c0hHIKPjxg7SLufuFPwdMXH
XvzsgasD+/jYos4ohn9Y1J9/wtUPHrXLQD17APwX7Y/x1Rz+g6zbh98tzPkP
XAfIUXPy9An53Mf+2xPvp8S7H5/9mKf1OotmKc6g+E4V6jtUtbSesRL+Gylj
fj1nUc+8LfB+xvzSFLlTUS5lbRf2rzF932mshP9bDnA/FE9LxxYNar2wBeO6
P3KrbybdszqLscb857afWMuITXWtNfdH7S3wNiEerW2KT2lIU/wU97I1+VC7
Ddgs5b03I2+6Y03Yyx7uTRewlLaRk67EUEp+VmqMmtuIvZcQW0vycJvbGdTo
Foxr+WTtKcW7qO9D1RXVq8sVS4qatNjxyBTnfoPjLSnu10zHzvg83PGoFHep
M3waX+64wu2/3h6YQi9Vn7oQt2Iud1ydooYtdVyS4tv2Ase7UtSti6XRKerE
hY51KWrPcMf7UtStwY6bU9S2/o4bU9Swexwnpag9wxxrUtS8UY5rUtS28xw3
paiX/RzvSFFruzmuSlEnlKdLyZVimUg89zuOT1H/FNdoYlMsQ4hHuZxKPts6
Lkjxns9wnJfi7vRwrIZLexrBvo5zXJTiPe/teHsKfVBc3Zm/Qe9ECk1Y73hW
Cv1Z53hmCp3R/gaxR+XmXPLTyfHWFHWhLIUv+RmXYm/aV58UfuWzwrEqRW1Q
vN2IYazj2hT6oLyOJLe6Fx25G7pHR3OXWjm2SVFnZqW4E7oP81PcJ83XvTiH
u6G7cD73QefTlzNSPnqRkyNS+JXP6xwrU9Q2cbaDV+c/lDswxbE2hV7d7XhR
Cr1amMKX/Ex3vDeFFml/Y9jjhBRrPv+GcLwMP5o7jfmXpLh/uns9U5yHzuL/
YYKdTw==
          "]], PolygonBox[CompressedData["
1:eJwlk9dvzmEYht9HnDm2FV9tacQqam9qqxZtrdZerdh7E/+G8Q/Yau+9914l
hGgIwoED15374Equ5Pvye9/3fu4nU15ZUFErpZQPtaF+pLQWvuMZmIv3hCr8
PnTBm8At/DDMxwfCKfwxtMTHw0X8FSzAh8IZ/Cm0xovhKv4OluFL4RP+Gebh
feEE/hDq4TPhAL4XKvCJcAV/C3XxctiD74ZF+Bi4gL+EdvgO+IPXgSX4NLiB
f4CF+Ag4hz+HOXg3OIrfhUq8BK7h76EVXgSX8TfQFi+Dm/hHWIwXwCX8NbTB
p8J1vBpm453D2d2GFvhoOI+/gK7h/+i3I9AM7wPH8QdQGD5D39YdJoczURbK
tAO+Ar7iB6Ejvgpq8EPQP5yxsj0JvcMz1mx1Rr/wmTpLM+gRzkRZHINeeFO4
l9yJQXhzeJTcAXWnAfxI7tAWfDPk4d9gRngGyl6ZbQtnqOx+wxp8dXhW6sS6
cAc0e31zejhTZakZ5oQ7pO6oI1PCM9OslPkGfH24y1+gNNxBdU8znRSeqWap
Tg0Od1pdPg3Dwx1Wd89CI3wj/Ex+k7JdGX6bMm6Izwpnvy95VzqFs9fObMW3
h9+6H7Lx/PC3nyXPbnn4rpqhsh0QnpUyzuBDwnd7At3xLLiT3Nld4TP07X+w
M7wj2o2/MDbcOXVNOzIuvDPaFe3shPAOa3fV4VHhHdFuqKMjw3fWXbUzw8J3
0l2049q19uGztHObwpkoi1/J2TUOuzJUl/LC3VKn9JbccPf1pv/zDKGd
          
          "]]}]}, {}, {}, {}, {}}, 
     {RGBColor[0, 1, 0], AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwt0tlSDQAAgOETZWnGCBVCGyrKvhQp2ihZctBiCVFKG0W0qpRSpEQqpcQL
dOXlXPQ14+Kb+R/gT6hqDjaFBAKBv2wSq73EC1FFFO85yCAH2Mpm3tBOC62k
8YFUHhNNJzvoYDvbiKCLbl7yisOMcIgn7OQtu+ghhki20EsfbbzmKKMc4SPH
2M0JZjhONSf5xCn2kM4sp6khgzHOEEsmPzjLU87xmSziOM8c2dRygXFyiCeP
eXKpI58JCkjgEj+5yDMK+UIRiRSzwGXqucIkV9nLdRa5RgMlfOUG+7jJL4I0
cotv3GY/ZSxRyuoE5UxRQRJ3+c0dmrnHd+6TzAP+UMlzHjLNI1J4xwD9DDPE
P28ts8gC88wyxQSjDNBFK3VUEqSADNKIJ4pwNrKB9awjjFDWsub/1yuenkGm

       "]]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, Automatic},
  BaseStyle->{FontFamily -> "Arial"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox["uE", TraditionalForm], None}, {
     FormBox["uC", TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 1}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.834514495672421*^9, {3.842274525902257*^9, 3.842274551576386*^9}, 
   3.842274722766551*^9},
 Background->None,
 CellLabel->"Out[37]=",ExpressionUUID->"853e50f5-f1a6-4178-9cc5-1de55872be0c"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"r1", ",", "r2", ",", "r3"}], "]"}]], "Input",
 CellChangeTimes->{{3.842274554540286*^9, 3.842274564219358*^9}},
 CellLabel->"In[38]:=",ExpressionUUID->"667ba917-4317-4dec-b2df-542439b636e8"],

Cell[BoxData[
 GraphicsBox[{{GraphicsComplexBox[CompressedData["
1:eJxtmHlUU3cexSNFW6rBjlrFgNatNlLcKKHqjH6RIFZlNDhKFSH2aBdqUJRR
adBOOSpxqVqLQa1FS9QyroQqrQUjLkXF4IJVCARkS14IgbAEXGthHufM/cVD
6j/3fE4k7977fu/7vjB8Wdz8T9wEAkFqD4GgS7nHJ2RqQxF9/2Njm0RtoKC5
c0ZwkgeMo/fcieQk+fT4xa5rKoeR+tPfOlQOPWP8PBg/D8b3gPMLR4Vvc5yk
8AXj9xtklTT3t0FxWx25jHE9MK4HxvXAuB4Y1wPjuuDXTp3y1EyPnVayMCZA
bahmfsDwA4YfMPyA4QcMP2D4AcMPGL7A8Lelx+RDOlUt8weGPzD8geEPDH9g
+APDHxj+wPAHhk8w/Pqf8ejFSczMLxh+wfALhl8w/ILhFwy/YPgFwy8YfsHw
DYb/6kXGVUKOY/7B8A+GfzD8g+EfDP9g+AfDPxj+wfAPhn8wcoCRZ7f7KYNE
XcfygJEHjDxg5AEjDxh5wMgDRh4w8oCRB4w8YOQBIxcY+f6RtSFILq1n+cDI
B0Y+MPKBkQ+MfGDkAyMfGPnAyAdGPjDygZEPjJxg5LUtCTuhcthYXjDygpEX
jLxg5AUjLxh5wcgLRl4w8oKRF4y8YOQFIy8YucHIf+DVIf20mkaWH4z8YOQH
Iz8Y+cHID0Z+MPKDkR+M/GDkByM/GPnByA9GfjB6AKOP0HP2RIOsifUBRh9g
9AFGH2D0AUYfYPQBRh9g9AFGH2D0AUYfYPQBRh9g9AFGH2D0AkY/rfI8U2dn
M+sHjH7A6AeMfsDoB4x+wOgHjH7A6AeMfsDoB4x+wOgHjH7A6AeMfsDoB4ye
wKv2nb1rkN2hD8IvRsmlZdT0zPtrIVfM2G3y2O1qw2U6dlw+TqupYP3gcyg+
11ZMe3OH4yB1dHb9q2J943MoPsf9AkPPnCqKEWtrmC75UHNPJjcx9XglfopS
yDE9nxl8VKOzMP0ksn8fvcLKtH8v81qHyMb0yk/ZD0X6BqZx0cmhUqWdqY9H
hFYhbmba/f6jJyjOAbR7PzgX0O594JxAu/eBcwPt3g/OEbR7XzhX0O794ZxB
u/eJcwft3i/OIXvPdusb5xLavX+cU2j3+4FzC8V9GXTsXIRSWEz587Zc0uhK
mMb/sUCsV5Qyfeu/b3/rEBmZ3p7/+JlIX840seP6MqnyIVPxyf2FCnEVU5tv
gEivqGZ6NSltrlxaw/RAiftmh6iW6Wq/ledVDqeGbipuEOlNTIeWTh2m1ZiZ
to/NWCBVckxvbfbcbpBZmB4pW39RIa5jmji+qrWz06my5Jmj1QYrU3G5NlKs
rWfaMcHrG53KxrRElfSbTN7A9HSF9YlZ0sh0i3+4n1JoZxq5LecjIefUiZXD
UzW6Jqa9AnbclKibmeJ5STmSmKNT/c6eAzCek1/r4vwC1TfYcwBm790BxYWc
JJs9B2A8J5/uDUqVS2+y+9qmWH8oUH2e3UecE/j6/GmtXSa/z84J5jv88cPI
RyAoYucG7wf4yzL/2apTXWXnBu8X+PMcUnvGIjnGzg3mBb4f18f3Y/52n69g
5MP1MI+7z1cw8uP6+H70edlW86lWc8vl+9DnCVFOUKksj/18UUWaQv4Xfx8o
rfLMMJiLKcTd+6hhj3Pu+Udl71H+xd8L/N/PShDz+fH/HKF5Ier9Bjq8RLdW
3Oyci15/XBCXtbj+fWD37vBZSr4P/L+5h5cO5B4Z6fnZ6FDlPytd5ufxRwKr
ZEEl+SZ39uaKq13m56bR96KkZdVUOqForTC51mV+Tvhdc02xrZZUFekPJQFm
l/lZ+WX8OPX7ZgrYtiZUbnKdnzvHSPfrLBzVvhesVaXUuczPKcX9BVxqHe2p
6uelnV7vMj+tSeYY4Yx6mva1KcnQ4jo/9/n9fE/SbqPGwOz6zh8aXeZnSGny
FPnRRjpYu2W+eF6Ty/x0bI44qprP8+6FF2QdzufV+PlIuZK/33u8Pu3h3fUe
W39lnZa/nw1xr95X8L2vWt6iy2rl980bx4/p+J5Xuns9V/nfpg9ea2mVZJdS
8Jhzy+SmYjqo33JOZykmt7lvf7Y1/hIdWfrtb6p+FWR9/cl73i/tD2lJCwK0
08uoffjoZOEMfo8IWxHziiCFOn7xTzWsrqKnkSajttW5Pzw1POopSK8gmfJG
hPxoBdsf2uL5txh//ud33f6OKlrQdZsF1bRImvZAlllDbgN6/tnJvyd7NcdO
VUabKOtS8SfipSbKPjg1Q9OHI7ki447Mk6NloZ599Rcs1GdQwvvKixZ6w1GZ
4FhhpdyrM9M1sVbKO6StFvHvs5hVXh56bxvFzkqaJb3ZQANF9Wsc/HtM9Eh2
VvGFnfKv5RhFiXa6nj7cW/1OM8Wv2SGVjnHuW1E1M6Vq/vn3DZg2VfvSPrFv
0ewp2ppil/0i9Ht7ucJazJ7XXyfVD+BullOLtqG3IKjMZd9Id+w+owstY/tY
vz07LBL+eQ57dddD2eEKl/3jxawJXwl/rKCMLnxR5bJ/nOxaM3pUk1v5Z0HK
MzUu+8eijZO/NWTVUJZq0glNlMll/+glfv15J9+33N+jn74357J/ZN83LhP3
5ahPZVmiI9fisn8s++pUoSzPQrnbT5pEK6wu+8cb724MUK60UoxkQ5h0sM1l
/8grCTuk8bHRwJo5PysKGlz2j9hNQ3rpCxsof6fPUHWC3WX/EI1rWuXg7+fq
SXaVbrRzLxQvHmrS8ed5an7gCinfe8uDzYWqFP73pfULfbhUAymfTJzcxs+v
77xL/67me57t9/YhQ4uRjtx+bJe0G2lCP/Fdj3fL2fnN2Hx9lXheJb0YdeCS
an4lJdnfS9eVVJNfXlqBorSaxn3n7sHx88r4+cqJ6q21VBFSvEbIz6Ntb047
qAs0046WH40SM0eBVzxf4fg5NCltvVS+t47MsVUKIT9vLDNnnlYF11OKV9YD
SUg9qdsGvalttVFQftJUeZuNgn+wfmlIb6SmuPAM1ZFGapmdY+nk50ea94i+
2vAmyipcKJHvbaaotLxoYarzvXrWva3VnOp832dKEowK/v2K9/nru7xO60TO
9/nH5pwvhfz7uWXzBKn2pXmTsvjMXPWTe3Qw6q3twp4Gij9ZTt6BzvlzfmFx
BLfuOrWfuxwpX26kj2f/2zcwu4DS3S4tkd8vI1HAYm5pm3M+7Wv5YVjIs59I
1nuZn/bKQ3KUld/g1v1Czxd/NFY70bkfFJTlGUS+zRS5fLhGleLMhf0A/mz/
2virbuB9Cjk1ZIdwXQnbD+CP/T7+/+uLnv1827D6LjkimhySUaUkuZI/Xs8/
76bY1AMO/v05p/zeAEGgk72+8t7vyC2h1cPsdbK8Esr8fsJW9ZMrNDfzm3zV
pnLqGOU+OPyl83bdt7FTtKKUht66KFOuLKW1JboY6WAj6RN252h8jNTn2Wm3
Gc80dLznxH2G6koa+c2fx7xfOo8jNu26pygopy9GLh2hLyynO2PlU9QJD2n0
3fFfOxIf0saycUd1o6voQaKgXeRbReqrncmOFdUUtDcmUxNbTbFeMbmq4Bry
qi+qlYTUUPCqoiYRv/82TZs8SO9dS6L8SSO1rbV0Xa2ZI2+rpZbBmg+lN02U
ZvNIcvB7cEGcx05DupnWBsVnq46Y6fC1NZcVX3A0e5/RKkrkaL23sb1znoVG
NAYP0YZbKGxN8Bj1O3X0dPqpcOmYOhp142S0mN9/i/b3VxkEVnru0z9FV2Kl
DPuGXEWplX6P33BdlllP/5Gamzr5Pfh4gem5OdlGC74LG6neaqOkoWHjldEN
5Nec/aF4aQNFrM1eLuT3X7cZQ3bqAhtpnN7ngKaPnYwHky/LPO3kPiz5lsRs
p6wWe7uZ34Mr1tkF+gtNpAqNGKO82ET/A6HeXmk=
     "], {{
       {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1], EdgeForm[None],
         GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJxNmHecV9URxe9dii4gK0WK1EVUiiViUMounRVYQIQliAY0FGm6BGFRpISl
CAgq1d4rIcaSGIMlhhBjbzFGjbGXGGs0FsgKknM+9/s+b/84O++WmTkz976Z
99viKXPHVhaEED6JIdSR/K/wlXCzcIhQT6gvdBPaCx2EvcIe4SPhG+Fr4Vbh
eeEB4RahMXoHCd8J3wp3CC2ERth+Udgp3C78IOwXPhX2CQ2EZkJznhsKhyGt
/71QI3ws/A+enm+JbftvhSwSivHt9XZCU+y3FZowbl/LZ0f8WacTerZ3mjBI
GCwcwVxroQu2bPsl4SFhu9BGOBQfLwuPCDuEo4TDWT8aaf2u2Mhy3Rwe3Rmb
11DhZKGX8GP07f8Y1h3rcfA2x1eER4VfCccz11noiV5X7HXDTwn7fiScxHo3
/HXHT2/ksUIfpH32RVq/h3AksZ6INNfXhF3CXUJ/4QT2noLdPuS3J/4HEafH
Q5gz3zI4WWcYen3hUkz8w5lzTKPxNwD0gNdApH2MYG+p8LqwW/i1MEroh/6p
6FvvWmGtsE4Yw5z5XiOsYW0idsuFKmGaMF2o4Cwdx0+I33GMZ87jN4THhLuF
SfBwHFOEscI44afCSNYms26OZyHN62yk7+/PkNafgE/n6XSkuU7FdgV8K+B1
Hv7MZSY6jm8GtjyexdwZwrnCmehUojcZ20OIfy5z5juf2KaSo/Hk5hykffyc
vY5pHtIxnY+0/gJsTCPOwfhbiF3b2ygsFZYJvyA2c7yAdce0iDgd0zvCk8K9
wkXMzRaWo+c4VsDPvHz+VfisZt1rK1k331VIx70aae4XI6vgOIdcLkOa61ps
m+865IXCZuIxr/eEp4X7hA3wXixcLizB9qZaOdiCXjW2ZhD/VuYc31VwNcdL
2Oc8rUfaxzb2OtYrkI71SqT1r8aG8/SW8Lhwj3Ad8dj2BDj4PNxz3FvcV1w7
XU9cX27At+P7QHhW+I1wI3OXCneG1KMKQ+ob7hnuJe5L7jmuz+5j7ifuN+4n
7iON8NeY8e08u6/9MqSe5vp+W0g9s25I/fMjbLinuZ81xa7reTPGrvNtQ96r
st5l6b6znfki+HlvG/Tbo5/1j6Ow3YE167nOdwp572xOPC0YGx3R8Vwx45bo
ebwDO62Jsw25N5cu+HbddM0biG/3iq7ouQ8ch4/urNmua3TWr45hLeslrudZ
r+oc8v7Ro5a/E/Hfk7H9ncTYPk5mbLu9GHeCS+a7N2tZ37LvcfoO2io8g63B
xOZe5b7gnuh+5T7QD04DGGe9ZADcS7DbDn5d4DQEu+4rrrWl8BvKmnm59pfB
aRhj2xrOuATd4cRQhv71Ib036+HnnjASTu4Fo7HlWuramtXNecTpvjAGHq7R
We0cy5p9T0Df/l3nTw95v5kc8l7ifXezpxwe7gNnwmsS49HoTSKGCnyXwbEi
5D3srJD3sLOx7To8k72u7Vm/spwa8v4xHV7e65pyD7pncG5jiMP8XGPn4M81
ey45mILvcdieQgznoeMYKhlPRq8y5L3KtvoT8yj053MGrpsbOLusx8zGzwL2
OQbX/CpicO3N+tUFrJm36/oSeLiuLwt5v1kZ8l6yFH7L2We+1YzNdQXjeeit
IH+L8D0LnotC3sNWhbyHrYaH+81l6Ph+Zv3Kcm3I+4fXsl61GH4b0Z9G3Avg
vok1c3K/2cr6GnwvxPYaYtuMjmPbwrgavS0h71XbQt6rriCGqxjb7tWMsx7j
/vIUvBfj09+F1xCP38msp13H2oGQ+kX2+8p9zT3FNdY9xfX+d8z791mMqX+5
h3weUr9xr/l9SD3Juo15LiKXN+Dbv0X+FlIN3wBnrz0YUo/x/j+FVI9ci0bw
7Hf2DyHVYtfh43l2DXRt38H8n0N6z73fd+Fe8vAEz77DrfBlDg+j6+eJ6Pod
/EtI76SfZ/Fs3S9C6t+F8L+PnN/IGdyEnyfx6xzsJK5n2OP9Hwq/Zf9z5Ne5
/VK4n+e/kpcicvYwPF3/d5Gf7eS0NXE8Qh7+Tr46k6dHydWrwh959nk+wHm9
wPn5+R+cQSm5300+HevT8Hc+Hicnb5Iv5+ptcu35d0O6i87D+8Rv3RJiKCUf
z5IH+3+RnP2TsyznTB7D/v3kq4CcPV8rjhfgv5PcFRHrXcTrXL5EzpzLl8mZ
8/QKedqN73Jy9iq6u8hLKTZfw66/r/x/igjfN+DrfLxVK09vkxPn5p2Q35F3
yY/z9F6tHL9PrnwOr3MOzt8HtXL2IXn7d0jfeObwSUjfjf7/yV7WnKsGzPv7
67OQvjMPDul/F9bx+9uE+abk+01y/gM6hfj5Cl/1GNdn/Qv2mNd/Qrrf9v8l
HG5Gtx4+P8eva6Lrk2uZ76/vru+tv4G/IZaPGdfB5l587+G5AO7fEtcdPDch
Zn9j+3uzIKb6tQ8O32N7ruZmxhSHgzuAjtdqsOH5GjjXoFsXP9/hN+OU1dE9
cHNO7bMZeTiA7mfk91P8HSA/NfjL5vzs/a675l/I89ec7T7m/Xyu5IyY4rPd
/fh37Ps5o/34zZ4LyMlsyXNiyrHtFeDLv238fenvzY6a2xbTHVwpOSimGtRX
8qKY3sv5krNiuidV+jOH/+ltlmwX030fKLkipveykWR1TO/uwZJLY3qnF0pW
xpTb3pIXxlQD6kou4vm0mOzaZlvJTTG9W4dLXh7T+zdU8uKY3vVDJVfHdNe6
Sy6PqaY3jikexzJYclVMdaK95JaY3kvHV48Yu0guian/2ecY/NrnaPxabyy6
l0qOjOndst5B6DqOPsTiuEuJ3TmuIM/rJUfEVJ/WSQ6LqSbVkexEnhvqzwRi
OVqyhD7tWIcQr3M8gDz3k1wWU91dI1kWU81zPhpiZ5TkZTHVs1MlN8ZUw8ol
N8RU82yjEDsNJDtzXq1iitnxdpPsH9N7ab2W6NaXPIK7YT+t8eV7cQh34xLJ
w2KqhT6TIs7F59mEM10r2Sym2jxP8vyYatB0yWLhXyHdnbrk2T6r8Hus5MSY
fkMfKdkrpnfUsbYh3uYx5d05d3wLibFFTGfjc2kaUx6dQ/9u6xBTzR4eUwzm
77Oaz3mdEhNvc/ZdGM99+D/TwHpN
           "]], PolygonBox[CompressedData["
1:eJwllMVuFVAURc+hQArB3WWAE9zdWqwtUNyhuLVASBgw4FsIBMIYd7fi7i4h
BAshjFk7e7DSlZeX3nP32fd1rqgqr6wREcVQE3pkxBg4ht+FpXgzOI8/g1H4
LviL74Ud+HA4jP+GIfhsuIl/hC14d9iH/4AleBM4iz+BuXgX2BP+ziK8AZzE
H8BKvB1cx9/BMrwFXMRfwAi8K9wOz7AcbwWX8VewAK8Hx/F7sBhvBKfxR1CB
d4Bq/AMMw3fCn/BMK/A2cBV/A2Px3fAP3w/z8EI4it+BVXgnuIV/gtX4GviC
H4LydEbKRmcOxmfBDfw91E9noLufgJbpO+kul6AEHwivw3dsnZ5Rs12BcXhv
uB++8wS8DzwMZ1onnYmy0I4bpjNRFqdge3pmzfoLGqd3pl2dgTJ8ELwNZ7IJ
3wzf8AJom96ZdnUN2qczVra6Y8d0RspGHdmAb4SvOH9iZjoTZaGdr8e3wk/8
O8zBh8LncMbr0h1TtzRDJV6V/m5huBsj4Ui4I9PwAfAy3KEpeH94Hu74fLxb
utsHwmf3SrtmWIj3TO/+IPTFi9LZPA5nMyN9d2Wks6and6cztbvS9K60w374
ZDiHPw3PMhUuhGdqnu68uq7PaqU7p67pTupi3fRd1clJ6Zk0izpenD5D/1tv
rmn6Test60x1Y2J69+rI6PRvgN6+Or02/cb1trUjZb8t3Q3tQL8FtdO712+C
3k5Bepd6Q+ri+HSX1cn/piCTNQ==
           "]]}]}, {}, {}, {}, {}}, 
      {RGBColor[0, 0, 1], AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwt0ldbDgAAhuGvyyxF02ioNEhJSdJQhNIQCtFCSlKKkoyMaEi0EyFOHftP
fke3i4P7+H0PnuT2/vq+oEAg8JtfhBNLGrmUUkszPYwwwyxzzLPAIkss84kV
PvOFVb7yje+s8YOf/3f/kOxEJx0M8JACLrNMErdJZIkGjlHPIntpJ4EFLpHP
ReaJ5xZxzHGBo9QxSyw32cMM58mjlo/s5ga7+EANR6hmmp20EcN7qsjlHFNE
00oU76gkhwomiaSFCCY4y2HOME44zexgjNNkU85bttNEGG84RRYnGSWUawTz
klIyOMELtnKVTTyjmHSKeMpGnlDII+5wj24G2ccGhjlOGlfYwggHuE8XvfQw
RAqbeU4J+2lkG6/J5AF36edvrI9JJYRXlHGQ6xwK+tfxOqkfPwo=
        "]]}}], {}}, {
    GraphicsComplexBox[CompressedData["
1:eJx1VX9M1GUYP0kKKDAsY4N+LIR2DgRSj5IpDwj+GPLj+GFrCrfFilUQ1xhZ
sLXZ0ENJGtoJRgJiYALqKTALPAMUDzsCDuX4Hr+R+maJCZweASElW5+X7d2+
f9yefd573uf5PJ/neZ/vq8nquPfsZDLZ6OPfonWoqXEpD00Lhm1t94o9aK0m
4OhrbupcayPDSQWdu0VFK8PPkeuCxmpkWJyuUmoFE8Mh0Ts9RUUPw99W3nug
0AoMT8/nX9dY+xmOTfAvEpTDJMWLt+AHDH7A4AcMfsDgBwx+wOAnpQvPg7fg
Aww+wOADDD7A4AMMPlJ94XXgefAWfIDBBxh8gMFHag543fm6+by8RX5g5AdG
fqm54/vM68zXyeflLfJLzTU/R3wfeR35uvg8vOXfDT+n/JzwfeJ1498dP+f8
nPHvsuKMyk9XPki6weBVedZiho8/b24XFfW0I/ZKoiqsj+w2rj2kFZpZHGD8
/+MdtW+gto3uz3p86SyaKb2wtktQdjKewPj/6KnsBr3mpuQ5sFQeYCnewKiT
94NOiAOdkZefG6m6cE9KJ+ThdQav7AVDcljWELm8NHbud0UFvfK99xGrez9d
+O3RlF5zlVpj9jeV63vpcdNelMlMzOIcfriHOIgrry5qT5WP0P6FjJJ8ax51
xE3PuhsHqMq9IcSi/Iky/kmQG1Mt1Hz3doqu/Bdyq6h7K8vZzCzO4Yd7iIO4
yIP5AubPwWvzlfVzWuEC4+20qeyas3id1aWPXHlArutmFufwwz3EQVy2N/73
R3/gj37AH/PGW7zXD2bG/lKqbjE9+DkHxjykfB1yTBX2M9OLf0fAmI8HqXtL
ArU/MD2l/IDhh3tSPIDBA7xQF+pEXahzrjZpW1bUMNOpdI8+Uz7Rz3QMX+7x
nVAgMJ2/2mX5V9FjYnkRv/DtiCDdbTOdj/JpkPu3sjqwvyZ140/LQvrIsaXI
eEyoZu8B+yryqfwhZekgnc4xpMtjhill27DzCrGaIny9S4TJfnpyYt4xd10r
TfbktGuOChT0WqhcNJuYxTn8cA9xEPdM9IJM9vEIbW9b7Z1nPUKRfyf661wH
qfps/Atr1E2MN+rYe9GrSm8YYO/cWnY5SVVnoavrdd1CQAfTAbq4784RUt+3
sL1RuMP9sPMuM7O4hzh/nM3q1qzrYX0FD/B6d/6jROO4kc0JeKOOyc3dPi5i
I9s3eHenF+H8CNMZupscDNnO3iNsL0F39AFxEBd9Qx/BA7zQd8wB6kBd0AG6
QEfoCh7gFbdId2GEYgp8W1ym3WlGsNnLTg7SI/vLz4bPNtKJfQkbdKF99OHo
WLRoM9KWNXXJql/NdO/l+gjtGz3M4hx+uIc4iIs8yAvdoKNnks1V/GRo6Tsz
5fONEDBExW86fZqbcYlu2CeEaN362fuLKM5r09zso4YsJ7vc2huUdr42Payh
d+k7t+VPmyKzl7Tj273EVbfYOfwQB3GRFzwwx56xqdplshqaO7X2uNA8RMEx
nV2i7SLNbN0ToPtigJ5ob1OJtha2P+KnD90RTM3sXdS/frImbGMrle48bNAk
99MzB9L8AjMNdOJ+g0q12kJqeVSlXN3F9klT/Ao/j4hO9s4K2lZmhfWZKPzh
3YeK5QJtmvqsUT/dzc7hh3P4oW7oAB7gBR7gBR2hK+qGDqgbOqAv6FNcyfwy
2egwJV1yfOegtZx9n4LavBwq9cVsr2HPYe9hD2IvYk9CF+ik/XxsVrFhad6g
O/qgL6Fz8oml+QMP8Crbl9lhHF+aR1ic4x7iIA/y/gfzqUZS
     "], {{
       {RGBColor[1, 0.5, 0], AbsoluteThickness[1.6], Opacity[1], EdgeForm[
        None], GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJwNw1dTEgAAAGD3ZCMb2bJlylKG7D1VphNnuRUEQRDEh67rUcuyeXVd156X
mS0rW5alPXief6B/kd/dR45O+qIlRUVF/06WngQUl4BKy4AlpeCycmhFJaS8
AlZZ5TPU4KAwFBCEBoExYAgWAkXUAuBV1RAsJpuBExFIEhJFRqEpaAwVg6Vh
cWYTHgkAMghEJpHEIpHZZAqHQs3naLwGOqEOwaczBAymkMmyWbmSRp6Ux5fx
BdoWpZjNSal0coGwEBOrJFJlk0QtlfkdioRcoxCK6PWEoEtdi6grAQJKQcAy
MKgcAq6EQYsBtWGPtgIKEbHYACQCiEKC0CgwBm2AV1fBYU6/whHQLGaNCCKh
jlCPJBHRFDKKTMJQKSYaLsehEpkMEotJZrMoHHY9vQFLo2ZwsKW8hc7nMQR8
plDAEgnZYpFVzGmUNBEYdL5MKpDLhAp5TCFaLtiblC0SlbKB1yhVqxJqmUMj
T2mVLp1qMuj26Fu5Vk6LtrW7wzoV9na2GfAmXJsR3tVuarfQ+s34GkM1FIcd
sHGjvpqegCPuVKhdqt6QS+tpVTjkHXZx0q2e8DsNbXBTO83aIZ7u8ae92r6I
hyeVzPQGH68NrD8dbBs1dozbbd1WX6fBE9F7w1p/jyOZ0Ji7TLN9YVdIl0iq
55Oa5LxbHBNNzYbjseYru865ePPIKf3kTDAWV6hSyolpvzwhi88528cswV5X
ZEQfGtbBMzBajorL4I1ZeCaLzy7iOTmuvSDOL3Fzea4lTyssN8cKzZd2bNvn
zTsXbemULpXWhvs8gSHN6WjNq03LQlp3/CDyesv+csOYXvAe3Q8d3gusfRlw
B9VP1gfPHUfefxg9exS6tufee+veuG389Hls/533z++pM4eBy8+Gbh34fv2c
/PF94vqLkdWP/Q9Xohe2zbtvnDf2vQc3fY9W+68+H175G926a//6bXzzjuU/
dnTQtQ==
           "]], 
          PolygonBox[{{150, 143, 93, 71}, {201, 187, 75, 109}, {147, 140, 91, 
           97}, {199, 154, 74, 100}, {171, 156, 102, 87}, {195, 151, 71, 
           99}, {189, 203, 110, 77}, {180, 141, 66, 103}, {168, 223, 130, 
           72}, {213, 145, 95, 113}, {214, 146, 96, 114}, {215, 147, 97, 
           115}, {153, 144, 94, 74}, {163, 142, 92, 82}, {178, 137, 58, 
           101}, {173, 160, 104, 91}, {166, 222, 127, 69}, {170, 224, 132, 
           75}, {179, 139, 63, 102}, {182, 174, 98, 106}, {172, 158, 103, 
           89}, {212, 164, 82, 98}, {145, 136, 87, 95}, {191, 209, 112, 81}, {
           221, 161, 81, 124}, {222, 165, 28, 127}, {217, 157, 77, 118}, {219,
            159, 79, 121}, {205, 214, 114, 111}, {133, 192, 107, 83}, {223, 
           167, 44, 130}, {135, 200, 109, 85}, {134, 196, 108, 84}, {137, 217,
            118, 58}, {193, 183, 69, 107}, {190, 206, 111, 79}, {202, 213, 
           113, 110}, {181, 162, 105, 104}, {142, 133, 83, 92}, {141, 221, 
           124, 66}, {139, 219, 121, 63}, {146, 138, 89, 96}, {208, 215, 115, 
           112}, {165, 175, 99, 28}, {169, 177, 101, 58}, {143, 134, 84, 
           93}, {167, 176, 100, 44}, {144, 135, 85, 94}, {197, 185, 72, 
           108}, {224, 169, 58, 132}}]}]}, {}, {}, {}, {}}, 
      {RGBColor[1, 0.5, 0], AbsoluteThickness[1.6], 
       LineBox[{2, 1, 11, 20, 29, 37, 45, 52, 59, 64, 67, 105, 162, 181, 160, 
        173, 220, 140, 147, 215, 208, 209, 191, 210, 161, 221, 141, 180, 158, 
        172, 218, 138, 146, 214, 205, 206, 190, 207, 159, 219, 139, 179, 156, 
        171, 216, 136, 145, 213, 202, 203, 189, 204, 157, 217, 137, 178, 155, 
        177, 169, 224, 170, 188, 187, 201, 200, 135, 144, 153, 198, 154, 199, 
        152, 176, 167, 223, 168, 186, 185, 197, 196, 134, 143, 150, 194, 151, 
        195, 149, 175, 165, 222, 166, 184, 183, 193, 192, 133, 142, 163, 211, 
        164, 212, 148, 174, 182, 106, 10, 9, 8, 7, 6, 5, 4, 3, 2}]}}], {}}, {
    GraphicsComplexBox[CompressedData["
1:eJxdmHlUU9cWxqNFLWpoq61iAOtUGy0qKqHqq24kiE/lYfApVdTYpR2oQVCe
Q4P2laUSZ2sxqKWoxIGnqIQqrRUjTqgYHKsQEuYMEAJhCDhRS95lLb8TFv7j
+n3se863z93nnH0zdHnsvK+683i8E914vI7/ixZE+St1lfT2mTMequnR07Z2
m3xYozAynnDOvadFZGZcudAQw7dYGO91O6MTKWsYf5a1MVAqrmVsWxx6WuGw
MT7Uy6efWlXPOOSCPV4naWDcLM01OZ2NjMPnjzuok5RTXsGI8O2ODIJfMPyC
4RcMv2D4BcMvGH7B8AuGXzD8guE37ObA2G2OHOYXDL9g+AXDLxh+wfALhl8w
/ILhFwy/4Oev99xSOAy0dN+DSIsoj/kHwz8Y/sHwD4Z/MPyD4R8M/2D4B8M/
GP7B8N+f3mtXOLTMPxj+wfAPhn8w/IPhHwz/YPgHwz8Y/sHwD/7lZH2LSKkj
y/PTEqXuEcsHjHzAyAeMfMDIB4x8wMgHjHzAyAeMfMDIB4x8wMgnMGzOMIvo
KcsHjHzAyAeMfMDIB4x8wMgHjHzAyAeMfMDIB4x8uvpHXvAPhn8w/IPhHwz/
YPgHwz8Y/sHwD4Z/MPx39Ys84BcMv2D4BcMvGH7B8AuGXzD8guEXDL9d/cE3
/IHhDwx/YPgDwx8Y/sDwB4Y/MPx19QOf8AOGHzD8gOEHDD9g+AHDDxh+us4P
X5gfjPnBmB+M+cGYH4z5wZi/63zwgfnAmA+M+cCYD4z5wJiv6/iYF+ODMT4Y
44MxPhjjdx0P82A8MMYDYzwwxuv6PMbF82A8D8bzXeMxDuLBiO/6dzyHv3dl
xLU7O/5VkLp02gc7HSmM0X+dOCUdq1aVUvfJY3Yoddfon+FXlkjFeoo5cP6h
TvKAccMrr118SyEhHjrGg37uzKMoobqK6Ys/Vz2WSE1Md38rboqcb2H6xcyg
4ypNNdO/iuzfVyuzMr1/T/Nah8DG9Ou/ZpcJtHVMj12aGCKW25nu7R6hlgkb
mY71gA6f6I+QB+5r5In7AecsdMRDx7kGHeNDx7kCHesCHfscOtYJOvYhdKwb
dOwb6FhH6Kh76FhX6Khj6Fhn6KhT6Fh36Kw/eqPjPUBHXULHe4GOOoWO98T6
8Llbr6o0RTTwxIUIOb+Q4v6aL9TKipn+4f8++skhMDD9/rznrwTaEqbHt99e
LpaXMV2YcbBAJqxgum20v0Arq2T6jYTUMKm4iumHity2OARGpq/2XXVR4TAy
PWRzYZ1Aa2L64OKpQ9QqM9Nbx6TPF8stTL+3xWOHTlLN9GP69Vdkwhqmx4+r
aHY6a5guSZw5UqmzMl1Yoo4UqmuZ3u7n+aNGYWN6kSLhpkRax/SzpdYXZlE9
07dOCPeV8+1Mj9x+6Qu+xc708eVDk1WaBqb39N95V6RsZDrq9ND7hQUWUbbr
nn3DWNcW2frDAcqLrO7/qIn1DVDeYfsCjPf09f7AZKn4LttHScfiL2kUf7J9
BsZ7//al0S6RPmHj43sB7znL/HezRnGDcdfzDf0s6oY7DL15PBfj/EP+qEvk
j+9B1I2Hj/FctegEY5y34KlXJrYpdVnMT+/Pjt7kW26x+TSh/RKF6sfsvMK8
mA/+sd5grDfyxXrAP9YbjPVGvjhX4BfnDhj7APlhvtOCS4HFklw23jVb1ddq
1T0Wv7U97vAex07mH+ct1gPPwy/OX6wPxoM/3FfYhxj/synDzqs6/X5QaUyo
lHT63lXsqvDgdfr+Crlz6oSGi49Z0aTJas6gutheT2Qd5/f66+vUXPw+z6+7
eXHxhm+HS+VcfPvvE5J1q7n7LXRl1Fu8JDq27Kebin7cfRn20Tfb4q66+uol
2fvk3PM3bsX24nWut7ebmkXZxbTKzbNNMeE+q+dHpakyKTf+NzEPROpO9Zai
3XpBU11IQaMuLJeaCql16MhE/gw9pSbM91dP15NEfidCeryUXuqe9eClldL8
jtfFq6R5HWXRztXJ+z3+dnLrtVCc+lSSWUVZVwu/Ei4zUc/G6KnypSaSytIf
SDwslJ0yNV3V10J9B274VH6lmpaHeLyjvVxNOTdmpqmirfSuo3yDY6WVomI8
3bVeNso9rK4UcOf7AEHtGgd3jkfPSpglvltHebcuGQTxdhI8k5yXfWenuDU7
xeJRjXQ7baiX8mNXX+3512WhvimDfHze00s67Y+9e8NnyTvdx6iHl5Emg7o5
hyZ8mrVB2Ol+xv619n4x0Yvj4gqPdJ3ZdV/jnh7tP22quhMfWDh7irqqkK1/
sJvXcd0+7vvqF3uJzFrI4hwhucHKgzoWh/ma1HV9eIF69r6PLNasFTYaKM2x
95wmRM/iwo4sG2B5ZmBxyCe0154yyZFSth5t55eGyP9VTq9n+f3AP1nK4k49
41lF88tZHPZPesefX1ew/Tk60dnHUlhJGR3bo5urf9g88vESsb6SxWH/dC/5
JlB+ror1DcV+j9byE420cNPkn3RZVSzO70/VLdl2I4tDX5GlmHRatcTE+gpF
aVqZyN9MPYW925xcfSGu/Pu4scpPzSwOfYd0gns/bR8L6zv8t68JkZq4Onxi
WC58x9Wf7B4lPqipdsWhL+lbro935FSzvsQ4MUitSKqh5T+cKZDkVrO4KYX9
eZbkGhaHviVnR4ZJwNUz+pZ9Ff081dNr6d1PNvnLV1lZnDXBHMWfUev6HehN
XxMl2hgqHmRjfc20XaYEXRO3L4pCD6u8Xf3PAd/fHotaXXHoewZUzflNll/H
+p76gOxa59F6it7s01NbUMfigosTp0iP17M49EV5u70HKzfYWV+UYtw6Tzi3
gQRjG2Ic3P5DnGNLxHHFvAYWh3tr9SS7QjOykVriuFXm6sSvn/Ch+ycl0372
Kv6Hkqs7+Yvxk1u4/Tc1L2ClmKtj4aLBJk3H71vrF3hbknXU9HRLgSJJR8fu
P7eLWg002/ejw7omA70eceiqYl45pW+5HSOcW06+uan5suJKSrBPTNMUVXLn
6Krxym1GGvuzm7uFq7ftH0xL0QSYqTS4cA2fq4+A6x5vWbj3vLPppEFktpA5
ukLG597fpNT1Yun+GkryzHoqCq6l6pkzzyqCaikwL2GqtMVGypaBH6ibbdQQ
G56uOFZPQUet3+vS6inVa9g76vAGapp9qdrJrQ/yxf3xpfnS93zu/kK/0rbo
izHq8WXk0Jfcsaz7nXrv8TyrEbj6k7TuVxdLn+jpy9n/GR2QnU+Zog0GGdd/
oB8JPuOzk7+uiGz/3vSHZsATdr8I/BdZlrVkkKTPcl/19TI60HR0SPCrXyks
88c8xeYSyvzFb5vyxXV239UPHnJCyt0XuI/iMkrIKyCPWi9ci5SuMNDFBYUR
lnW3yRHR4BCNKCbBq9/u61Y/ZP3Dwk0nbSru/sD91bTFT6zm7peUJR/u4Pfg
+qdF58KULx7TktTcpfzkRsoqWCCS7m+kfH2uTjDa1d+g74hcMVSlSGqkUz3G
H9BVllPfV2e7z3ilYve3W9+iekNTiut7c++Cy5L2Rjrv1tJsTnaNg/oa/uPf
J7w61Vf7CLdB4Vx9maKTDzm481t0PW+cljuvwXNKHr/PC3hKg+9dkchXFdPt
0fVOwcpi+m74smHaghIatnnPY1l+CT2N57UKRlfQJv3Y45qRFeRZ+8goCq6i
aM+oHEVQFd1WquZIW4wkyJs0XN1spLWBcdmKY2bKj3XfrUsz07D6IB91eDWt
9zK0OudW06OD/RU6npVG3MlYKuT67f+KzQ1Orq/+M27jbUlmLfk2Zn8uXFZH
CYNDx8mX1pEhJfGaxMNOY7Xeh1R97aQIiRglv9JApevsPO3lBlo9xF4jyS0i
zx+8Djpyiki7Ye8llbeB1hZposSDDDTy4bhdjvgyejBGOkW5oYwC90dlqqIr
SXnDmehYWUkN0yYP1HoZKSjmUYOA+95ItbknOLjviqZBqs/Fd000+4DBKoi3
0JFba67JvrPQy+lnwsWjaih0TdAo5cc1lG7fmCMrtlKbd/8kTZGV5v8cOly5
zUan8k1t5kQbdZ/hs1sTUE8Ra7NX8Lnvgawme6uZ6/vdhiTeE5nt9H+h3dwe

     "], {{
       {RGBColor[0, 1, 0], AbsoluteThickness[1.6], Opacity[1], EdgeForm[None],
         GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJxNmHucl0UVxs8sIggLXlBA1BTUFJCrgLC6CKLckZsIosg1TUIR3ZQVVjBA
EFbRBV25uobl3exCRkpERkaWmpqZmppaCWUUFabU2nk83/cz+8f5PDPzzpzn
zJl5z/P+fu1nzB13TYmZ1SezRo5NHQ9xu9/bjRz1rIVbS4u25hwOHuJ2BNjY
7UjwULejwCZuxzNPz08A9fwYt8Pcmrn9ya3UrY1bW9riPbYBfztQ/MeB8tsK
nqZuR4Py+wV49Px01iiWk5gn/kN9j/WOW9z+B//Jbqc0iOVUULF8EVQs3RvM
PREe+e0Ij/baCVQsnUHNPQNULHvcOrh1detGW35Pg0dx94BHsfQEFcuZoOY2
8b0ktwe83Ysx7bs3qLj6gIrrLFBxWQr+s936MqYYz2kQ1yB8yUd/9q+4zgUV
10BiEv8ktxFuI93OY0zry/Gn9efjT7FcAIr/r25lbkPchhGXYhkOluO7nFhG
goplFDjA7UJQcQ3Fp3yNZkxxjWVvimUMY+qPY0xxTYRbnN8grvFuE/CrGC8G
Ne8i5gyFvyd5uYRYFeNKt6+5LXG7lFgV12WgYpkCKsbLQcX1TTimuk2jrViu
Yo3mPsiz2W4z2YPOZBaoWL4ETna7AlQsV4KK5cug/E6HR3udAcrvV+DS84fo
V7hdzTzxXwOKfy4o/mtB8c8DxX8dKP4FPNPYHHjk96twaewGUJw3guKcD4qz
EhTnTaD8Xg+P8rcQHo09DEeV2820xbOUmDR3EWPiXAyK8xawknOuhPMRfC2z
OPub8HUr4/J3O/3lbivwK3+3gfK3ivVLOTfdKd2hR1lbzfpFrL+jgb/VoPzd
Ceo+PgZ3jdsa2lq3FtS6u0Gtuwe8y+1x1tTSX0mM9zIuH+tZr3XrGFP/CeZt
cHvS7T63OreNjGvut+hv5rnam7B1+C7xOrbfogbq/dO7OthCu6Qd0pJCw9Qv
NEz9QsPULzRMfWlLK/qFxhR91fMT6J9IX9olvWhjWcPaWtawYy1rWDv4jqdf
aKT60rBj4PuzhUa2pn9SAz7V6k7wSSNOtqxbp1jWrVPhO42++E6nL76O9Atd
Uc2SNkoXO1jWLfHttdCtLjxTPe8Kd3f64u5Bv9CtHpZ1sSPcveArtKqXZa3q
zXrV9nLLutWHOPrSl+5Jw8qIQxpxtmV96k9M5TxTLAN4pjWq00PgVv0fCPcg
+oU+SQ+KGq36+JGFVg1mDwPxK+5h+C10a5hl3RpuWatGWtaqUfgYTV/cY+gr
lrH0Fcc4+oX+ye8Lblst7r24VasLLRKOt6xVE+DQPiYT00SeyZ90YRIxTaY/
gjHN+7fbH92+bVmrlBfVGNUK1Y6X3J6yrFfTLDRL3NPpF1oyHb8z6RdaVfRV
z+dY1pKrLWvVFZa16kpiuYr+y27bLOuh9GI2fHPoF1o4y7IWFv258BVacqNl
rbrWslbNg+96+q+6PW1ZD6UdFZa1qtAS1eebLWvhXMtaNd+yVlXCsYC+OBbS
f81tu2W9kr8qy1rYUJ8WWdanxZb1qeiv4AzFsYRn4lhK/3W3HZZ1TFqzzLI+
rYBjOc+0fqVljZE2SDPecNtpWa90X6ot65PuUPHdM8WyVq0mjlX4e8vtWcva
Ja2psaxbayxrlbTmbbddlvVK+lLL3HX0NW89/XfdnrPQqEL/1jJPGiTNes9t
t4VGac0mnsmkVdIs/a7Sbyrpwwduz1tonX5zfd0+L1/2mds+i98jjX3ggMVv
MdU7vfOqAd+3eL/1Lt/pcwa4PePtH1i8Z8U34YO0f0hbd1738CHaz9DWndSd
eZj2z8hNLXt8gvaPOfNlnNejtH9Cu5r8P0b7p7SLb4nHaf8IruKb6hHaP4dL
+XyRfWqPvyCvhe4/SfuXtOssfsd9x0If9Tvouxba9Cu371nUZ3FuZ48b8bmB
nD0F16/J41Ryto32K+RxNjl7mvZvyGMFZ7IVrt+yzyr2t4P278ij8vYOua4l
fzsZf5OcVpO/Z2n/npzWkKfnWKu87mL8D+RxA3t6kX29Tx43k7PnaW8lR4MZ
e59x5f4F9qQ87cancvMSOdtGvoo79jI5Uz5eJU/byVEF+XiNc9hBjqrIx+vk
aSc5WkaOXyHPyscb5Ek5eKtBbt4kP7vIUQ17fY/YleO3yZNy8y75202+NpDL
d8inNE13qw78gJw15r8Rvb+fWLyb+iY7LAXq++tDi2+11sxRu7nF3fyQ8U8t
3nF9P+lu/o38/53x9szZS1t3+S8W9/mgxfdt8c0hHIKPjxg7SLufuFPwdMXH
XvzsgasD+/jYos4ohn9Y1J9/wtUPHrXLQD17APwX7Y/x1Rz+g6zbh98tzPkP
XAfIUXPy9An53Mf+2xPvp8S7H5/9mKf1OotmKc6g+E4V6jtUtbSesRL+Gylj
fj1nUc+8LfB+xvzSFLlTUS5lbRf2rzF932mshP9bDnA/FE9LxxYNar2wBeO6
P3KrbybdszqLscb857afWMuITXWtNfdH7S3wNiEerW2KT2lIU/wU97I1+VC7
Ddgs5b03I2+6Y03Yyx7uTRewlLaRk67EUEp+VmqMmtuIvZcQW0vycJvbGdTo
Foxr+WTtKcW7qO9D1RXVq8sVS4qatNjxyBTnfoPjLSnu10zHzvg83PGoFHep
M3waX+64wu2/3h6YQi9Vn7oQt2Iud1ydooYtdVyS4tv2Ase7UtSti6XRKerE
hY51KWrPcMf7UtStwY6bU9S2/o4bU9Swexwnpag9wxxrUtS8UY5rUtS28xw3
paiX/RzvSFFruzmuSlEnlKdLyZVimUg89zuOT1H/FNdoYlMsQ4hHuZxKPts6
Lkjxns9wnJfi7vRwrIZLexrBvo5zXJTiPe/teHsKfVBc3Zm/Qe9ECk1Y73hW
Cv1Z53hmCp3R/gaxR+XmXPLTyfHWFHWhLIUv+RmXYm/aV58UfuWzwrEqRW1Q
vN2IYazj2hT6oLyOJLe6Fx25G7pHR3OXWjm2SVFnZqW4E7oP81PcJ83XvTiH
u6G7cD73QefTlzNSPnqRkyNS+JXP6xwrU9Q2cbaDV+c/lDswxbE2hV7d7XhR
Cr1amMKX/Ex3vDeFFml/Y9jjhBRrPv+GcLwMP5o7jfmXpLh/uns9U5yHzuL/
YYKdTw==
           "]], PolygonBox[CompressedData["
1:eJwlk9dvzmEYht9HnDm2FV9tacQqam9qqxZtrdZerdh7E/+G8Q/Yau+9914l
hGgIwoED15374Equ5Pvye9/3fu4nU15ZUFErpZQPtaF+pLQWvuMZmIv3hCr8
PnTBm8At/DDMxwfCKfwxtMTHw0X8FSzAh8IZ/Cm0xovhKv4OluFL4RP+Gebh
feEE/hDq4TPhAL4XKvCJcAV/C3XxctiD74ZF+Bi4gL+EdvgO+IPXgSX4NLiB
f4CF+Ag4hz+HOXg3OIrfhUq8BK7h76EVXgSX8TfQFi+Dm/hHWIwXwCX8NbTB
p8J1vBpm453D2d2GFvhoOI+/gK7h/+i3I9AM7wPH8QdQGD5D39YdJoczURbK
tAO+Ar7iB6Ejvgpq8EPQP5yxsj0JvcMz1mx1Rr/wmTpLM+gRzkRZHINeeFO4
l9yJQXhzeJTcAXWnAfxI7tAWfDPk4d9gRngGyl6ZbQtnqOx+wxp8dXhW6sS6
cAc0e31zejhTZakZ5oQ7pO6oI1PCM9OslPkGfH24y1+gNNxBdU8znRSeqWap
Tg0Od1pdPg3Dwx1Wd89CI3wj/Ex+k7JdGX6bMm6Izwpnvy95VzqFs9fObMW3
h9+6H7Lx/PC3nyXPbnn4rpqhsh0QnpUyzuBDwnd7At3xLLiT3Nld4TP07X+w
M7wj2o2/MDbcOXVNOzIuvDPaFe3shPAOa3fV4VHhHdFuqKMjw3fWXbUzw8J3
0l2049q19uGztHObwpkoi1/J2TUOuzJUl/LC3VKn9JbccPf1pv/zDKGd
           
           "]]}]}, {}, {}, {}, {}}, 
      {RGBColor[0, 1, 0], AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwt0tlSDQAAgOETZWnGCBVCGyrKvhQp2ihZctBiCVFKG0W0qpRSpEQqpcQL
dOXlXPQ14+Kb+R/gT6hqDjaFBAKBv2wSq73EC1FFFO85yCAH2Mpm3tBOC62k
8YFUHhNNJzvoYDvbiKCLbl7yisOMcIgn7OQtu+ghhki20EsfbbzmKKMc4SPH
2M0JZjhONSf5xCn2kM4sp6khgzHOEEsmPzjLU87xmSziOM8c2dRygXFyiCeP
eXKpI58JCkjgEj+5yDMK+UIRiRSzwGXqucIkV9nLdRa5RgMlfOUG+7jJL4I0
cotv3GY/ZSxRyuoE5UxRQRJ3+c0dmrnHd+6TzAP+UMlzHjLNI1J4xwD9DDPE
P28ts8gC88wyxQSjDNBFK3VUEqSADNKIJ4pwNrKB9awjjFDWsub/1yuenkGm

        "]]}}], {}}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, Automatic},
  BaseStyle->{FontFamily -> "Arial"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox["uE", TraditionalForm], None}, {
     FormBox["uC", TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 1}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.842274565075407*^9, 3.842274725161994*^9},
 Background->None,
 CellLabel->"Out[38]=",ExpressionUUID->"0cc1c95b-36c0-4b1b-9751-bbb2e00791a6"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Transform E-> p*Edata to represent over-inclusive vaping definition\
\>", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.781476201147489*^9, 3.7814762303612967`*^9}, {3.7814765934921217`*^9, 
  3.781476599052796*^9}, {3.842274795047114*^9, 3.842274838823551*^9}, {
  3.842274876881468*^9, 
  3.842274907797729*^9}},ExpressionUUID->"08f23630-6205-47fe-bef3-\
851508e50197"],

Cell[CellGroupData[{

Cell["Special case: \[Alpha] = 1", "Item",
 CellChangeTimes->{{3.781476265050036*^9, 3.781476288207602*^9}, {
   3.7814763275108747`*^9, 3.781476331218816*^9}, {3.781476402397532*^9, 
   3.7814764323909683`*^9}, {3.7814766701297197`*^9, 3.781476681448688*^9}, {
   3.781477784547036*^9, 3.78147781265799*^9}, {3.7814778624128103`*^9, 
   3.7814778946608353`*^9}, 3.781514150949218*^9, {3.78151422679467*^9, 
   3.781514290002342*^9}, {3.7815151562346563`*^9, 3.781515205910882*^9}, {
   3.834579124092229*^9, 
   3.834579130126539*^9}},ExpressionUUID->"dfe2d88b-f43b-4b63-b949-\
281ef52c42cf"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"DEsdata", " ", "=", " ", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", 
       RowBox[{"P", "[", 
        RowBox[{"c", ",", "uC"}], "]"}]}], " ", "+", " ", 
      RowBox[{"e", " ", 
       RowBox[{"P", "[", 
        RowBox[{"c", ",", "uC"}], "]"}]}], " ", "-", " ", 
      RowBox[{"c", " ", 
       RowBox[{"P", "[", 
        RowBox[{
         RowBox[{"1", "-", "c", "-", "e"}], ",", 
         RowBox[{"1", "-", "uC", "-", "uE"}]}], "]"}]}], " ", "-", " ", 
      RowBox[{"c", " ", 
       RowBox[{"P", "[", 
        RowBox[{"e", ",", "uE"}], "]"}]}]}], " ", ",", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", "e"}], ")"}], " ", 
       RowBox[{"P", "[", 
        RowBox[{"e", ",", "uE"}], "]"}]}], " ", "+", " ", 
      RowBox[{"c", " ", 
       RowBox[{"P", "[", 
        RowBox[{"e", ",", "uE"}], "]"}]}], " ", "-", " ", 
      RowBox[{"e", " ", 
       RowBox[{"P", "[", 
        RowBox[{
         RowBox[{"1", "-", "c", "-", "e"}], ",", 
         RowBox[{"1", "-", "uC", "-", "uE"}]}], "]"}]}], " ", "-", " ", 
      RowBox[{"e", " ", 
       RowBox[{"P", "[", 
        RowBox[{"c", ",", "uC"}], "]"}]}]}]}], " ", "}"}], " ", "/.", 
   RowBox[{"{", 
    RowBox[{"e", "\[Rule]", 
     RowBox[{"p", " ", "edata"}]}], "}"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"varsdata", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{"c", ",", "edata"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"paramsdata", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{"uC", ",", "uE", ",", "p"}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814765794072323`*^9, 3.7814766117715178`*^9}, {
   3.781477330698406*^9, 3.781477334652973*^9}, {3.781514028937169*^9, 
   3.781514058092772*^9}, {3.8345112824226017`*^9, 3.834511375199737*^9}, {
   3.834511412186735*^9, 3.834511599555974*^9}, {3.834511816205503*^9, 
   3.834511834067465*^9}, {3.834511895561866*^9, 3.8345119116918583`*^9}, {
   3.834512395209323*^9, 3.834512395702973*^9}, {3.834512440325427*^9, 
   3.834512440487258*^9}, {3.834579440938689*^9, 3.8345794415017757`*^9}, {
   3.8422748968881893`*^9, 3.842274953746295*^9}, {3.842275002865492*^9, 
   3.842275009536105*^9}, {3.842275127988654*^9, 3.842275131110653*^9}},
 CellLabel->"In[53]:=",ExpressionUUID->"60605db5-98e4-425b-998f-80caa384fb53"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"c", " ", "edata", " ", "p", " ", "uC"}], "+", 
    RowBox[{"c", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", 
       RowBox[{"edata", " ", "p"}]}], ")"}], " ", "uC"}], "-", 
    RowBox[{"c", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", 
       RowBox[{"edata", " ", "p"}]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
    RowBox[{"c", " ", "edata", " ", "p", " ", "uE"}]}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "c"}], " ", "edata", " ", "p", " ", "uC"}], "-", 
    RowBox[{"edata", " ", "p", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", 
       RowBox[{"edata", " ", "p"}]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "+", 
    RowBox[{"c", " ", "edata", " ", "p", " ", "uE"}], "+", 
    RowBox[{"edata", " ", "p", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", "c", "-", 
       RowBox[{"edata", " ", "p"}]}], ")"}], " ", "uE"}]}]}], "}"}]], "Output",\

 CellChangeTimes->{
  3.8345116029019823`*^9, 3.834511913686256*^9, 3.834512445702838*^9, 
   3.834579228283432*^9, 3.8422744704807673`*^9, 3.842274955447542*^9, {
   3.842275004588724*^9, 3.84227501016873*^9}, 3.842275131818899*^9},
 CellLabel->"Out[53]=",ExpressionUUID->"8eac1b0e-f123-4191-95fd-d15f9f7f0272"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"numEqns", " ", "=", " ", 
  RowBox[{"Length", "[", "DEs", "]"}]}]], "Input",
 CellChangeTimes->{{3.781474692441131*^9, 3.7814747167748203`*^9}, {
  3.7814765273261337`*^9, 3.781476531800425*^9}},
 CellLabel->"In[56]:=",ExpressionUUID->"c6a2e8ce-9921-4772-9e4c-9aaf537add02"],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{
  3.834511611870761*^9, 3.8345119160760517`*^9, 3.834512449287446*^9, 
   3.8345792298433847`*^9, 3.842274475584813*^9, 3.8422749603116827`*^9, {
   3.842275006443632*^9, 3.8422750148265*^9}, 3.8422751335012608`*^9},
 CellLabel->"Out[56]=",ExpressionUUID->"f1794859-77c8-48d4-88e3-e4d85584e52f"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot direction field", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781477421345943*^9, 
  3.781477441545467*^9}, {3.834578949649909*^9, 
  3.8345789498863773`*^9}},ExpressionUUID->"854c782d-2b8b-444b-a73a-\
b526a1e58fec"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"StreamPlot", "[", 
    RowBox[{
     RowBox[{"DEsdata", "/.", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"uC", "\[Rule]", "cig"}], ",", 
        RowBox[{"uE", "\[Rule]", "ecig"}], ",", 
        RowBox[{"p", "\[Rule]", " ", "\[Rho]"}]}], "}"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"c", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"edata", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Smoking", ",", " ", "Vaping"}], "}"}]}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"cig", ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"ecig", ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Rho]", ",", "0.1", ",", "0.9"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, {
   3.7814764997404947`*^9, 3.781476509740304*^9}, {3.78147661985701*^9, 
   3.781476620330504*^9}, 3.781477464922221*^9, {3.834511623489006*^9, 
   3.834511707445655*^9}, 3.8345792153431683`*^9, {3.8422749779959297`*^9, 
   3.842274981344179*^9}, {3.8422750240322647`*^9, 3.842275057352314*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"5e9f4592-8b37-4f66-b7f2-a07dcded9e65"],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`cig$$ = 0.329, $CellContext`ecig$$ = 
    0.244, $CellContext`\[Rho]$$ = 0.736, Typeset`show$$ = True, 
    Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`cig$$], 0, 1}, {
      Hold[$CellContext`ecig$$], 0, 1}, {
      Hold[$CellContext`\[Rho]$$], 0.1, 0.9}}, Typeset`size$$ = {
    450., {221., 227.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`cig$$ = 0, $CellContext`ecig$$ = 
        0, $CellContext`\[Rho]$$ = 0.1}, "ControllerVariables" :> {}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> StreamPlot[
        ReplaceAll[$CellContext`DEsdata, {$CellContext`uC -> \
$CellContext`cig$$, $CellContext`uE -> $CellContext`ecig$$, $CellContext`p -> \
$CellContext`\[Rho]$$}], {$CellContext`c, 0, 1}, {$CellContext`edata, 0, 1}, 
        FrameLabel -> {$CellContext`Smoking, $CellContext`Vaping}], 
      "Specifications" :> {{$CellContext`cig$$, 0, 1}, {$CellContext`ecig$$, 
         0, 1}, {$CellContext`\[Rho]$$, 0.1, 0.9}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{505., {303., 310.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{3.842275057677845*^9},
 CellLabel->"Out[52]=",ExpressionUUID->"c243464a-18b0-45b2-a6c0-f99eff34106b"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Solve for equilibria", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 
  3.7814722277070637`*^9}},ExpressionUUID->"6e04284b-303c-4a0c-842f-\
7bbae89eb16b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"equildata", " ", "=", " ", 
  RowBox[{
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"DEsdata", " ", "\[Equal]", " ", 
      RowBox[{"ConstantArray", "[", 
       RowBox[{"0", ",", "numEqns"}], "]"}]}], ",", "varsdata"}], "]"}], "//",
    "FullSimplify"}]}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, {
   3.7814764997404947`*^9, 3.781476509740304*^9}, {3.78147661985701*^9, 
   3.781476620330504*^9}, {3.842275118904078*^9, 3.8422751462292023`*^9}},
 CellLabel->"In[57]:=",ExpressionUUID->"1ae1de4c-6ed7-4458-8a86-d2722aeea5de"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"c", "\[Rule]", "0"}], ",", 
     RowBox[{"edata", "\[Rule]", 
      FractionBox["1", "p"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"c", "\[Rule]", "0"}], ",", 
     RowBox[{"edata", "\[Rule]", "0"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"c", "\[Rule]", "1"}], ",", 
     RowBox[{"edata", "\[Rule]", "0"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.8345124159976883`*^9, 3.8345124767271433`*^9, 
  3.842274486816431*^9, 3.842275146859703*^9},
 CellLabel->"Out[57]=",ExpressionUUID->"d30075e1-90ce-48f7-a9c5-d056d0302794"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Linearize system (Jacobian)", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}},ExpressionUUID->"d1a22f61-c8d3-47a1-9758-\
bea908fb52dd"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Jdata", " ", "=", " ", 
  RowBox[{"D", "[", 
   RowBox[{"DEsdata", ",", 
    RowBox[{"{", "varsdata", "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, 
   3.781472428807068*^9, {3.7814766355579987`*^9, 3.781476646201336*^9}, {
   3.8345124977788897`*^9, 3.8345125010778427`*^9}, {3.83451350271008*^9, 
   3.834513515332168*^9}, {3.842275175539514*^9, 3.842275180825841*^9}},
 CellLabel->"In[58]:=",ExpressionUUID->"b5ec988d-29df-4e49-b816-498bf993d405"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"-", "c"}], " ", "uC"}], "+", 
      RowBox[{"edata", " ", "p", " ", "uC"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", 
         RowBox[{"edata", " ", "p"}]}], ")"}], " ", "uC"}], "+", 
      RowBox[{"c", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", 
         RowBox[{"edata", " ", "p"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
      RowBox[{"edata", " ", "p", " ", "uE"}]}], ",", 
     RowBox[{
      RowBox[{"c", " ", "p", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
      RowBox[{"c", " ", "p", " ", "uE"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"-", "edata"}], " ", "p", " ", "uC"}], "+", 
      RowBox[{"edata", " ", "p", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}]}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "c"}], " ", "p", " ", "uC"}], "+", 
      RowBox[{"edata", " ", 
       SuperscriptBox["p", "2"], " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "-", 
      RowBox[{"p", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", 
         RowBox[{"edata", " ", "p"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "uC", "-", "uE"}], ")"}]}], "+", 
      RowBox[{"c", " ", "p", " ", "uE"}], "-", 
      RowBox[{"edata", " ", 
       SuperscriptBox["p", "2"], " ", "uE"}], "+", 
      RowBox[{"p", " ", 
       RowBox[{"(", 
        RowBox[{"1", "-", "c", "-", 
         RowBox[{"edata", " ", "p"}]}], ")"}], " ", "uE"}]}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.834512484043462*^9, 3.8345125025139313`*^9}, {
   3.834513503855713*^9, 3.834513516099617*^9}, 3.842274490466264*^9, 
   3.8422751814102707`*^9},
 CellLabel->"Out[58]=",ExpressionUUID->"0c17ad1a-0f92-4b11-973e-7a3e8a837c04"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Jdata", "//", "MatrixForm"}], " ", "//", "Simplify"}]], "Input",
 CellChangeTimes->{{3.834513529871863*^9, 3.8345135320437202`*^9}, {
  3.8422751837399178`*^9, 3.842275184119329*^9}},
 CellLabel->"In[59]:=",ExpressionUUID->"68341285-01db-4a82-ab17-14188cd97bd0"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       RowBox[{"2", " ", "uC"}], "+", "uE", "-", 
       RowBox[{"2", " ", "c", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}], "-", 
       RowBox[{"edata", " ", "p", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "uC", "+", 
          RowBox[{"2", " ", "uE"}]}], ")"}]}]}], 
      RowBox[{
       RowBox[{"-", "c"}], " ", "p", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "uC", "+", 
         RowBox[{"2", " ", "uE"}]}], ")"}]}]},
     {
      RowBox[{
       RowBox[{"-", "edata"}], " ", "p", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}], 
      RowBox[{"p", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "c"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{"2", " ", "edata", " ", "p"}]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "uC", "+", 
            RowBox[{"2", " ", "uE"}]}], ")"}]}]}], ")"}]}]}
    },
    GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.834513518640321*^9, 3.834513532671836*^9}, 
   3.84227449229935*^9, 3.8422751863894987`*^9},
 CellLabel->
  "Out[59]//MatrixForm=",ExpressionUUID->"97402d29-f18b-49b4-a545-\
606ba8e61535"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Evaluate Jacobian at equilibria", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 
  3.7814729129771643`*^9}},ExpressionUUID->"f80e7fa8-76db-458f-84cb-\
46a26f53cc60"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Jstardata", " ", "=", " ", 
  RowBox[{
   RowBox[{"Jdata", "/.", "equildata"}], " ", "//", 
   "FullSimplify"}]}]], "Input",
 CellChangeTimes->{
  3.781472180031621*^9, {3.7814722370203133`*^9, 3.781472263174114*^9}, 
   3.781472428807068*^9, {3.7814729275631857`*^9, 3.781472961260558*^9}, {
   3.781473116004878*^9, 3.781473141513234*^9}, {3.781473253486924*^9, 
   3.7814732759163113`*^9}, {3.842275192649678*^9, 3.842275197006742*^9}},
 CellLabel->"In[60]:=",ExpressionUUID->"c8cf73a1-4ff4-48c6-8851-4b7969d6d44d"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"uC", "-", "uE"}], ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"1", "-", 
        RowBox[{"2", " ", "uC"}], "-", "uE"}], ",", 
       RowBox[{
        RowBox[{"-", "p"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "uC", "+", 
          RowBox[{"2", " ", "uE"}]}], ")"}]}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"2", " ", "uC"}], "+", "uE"}], ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"p", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "uC", "+", 
          RowBox[{"2", " ", "uE"}]}], ")"}]}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"1", "-", 
        RowBox[{"2", " ", "uC"}], "-", "uE"}], ",", 
       RowBox[{
        RowBox[{"-", "p"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "uC", "+", 
          RowBox[{"2", " ", "uE"}]}], ")"}]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"p", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "uC"}], "+", "uE"}], ")"}]}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.8345125119552107`*^9, 3.834513620140758*^9, 
  3.842274494933042*^9, 3.8422751977123632`*^9},
 CellLabel->"Out[60]=",ExpressionUUID->"c6ff30d0-d33a-478e-9f5a-3888f2dca2f7"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", " ", "/@", " ", "Jstardata"}]], "Input",
 CellChangeTimes->{{3.781513774371643*^9, 3.781513812723462*^9}, {
  3.842275201416206*^9, 3.842275202929487*^9}},
 CellLabel->"In[61]:=",ExpressionUUID->"abefd052-73e1-4d38-83c5-06a05ccbdaf1"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   TagBox[
    RowBox[{"(", "\[NoBreak]", GridBox[{
       {
        RowBox[{"uC", "-", "uE"}], "0"},
       {
        RowBox[{"1", "-", 
         RowBox[{"2", " ", "uC"}], "-", "uE"}], 
        RowBox[{
         RowBox[{"-", "p"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "uC", "+", 
           RowBox[{"2", " ", "uE"}]}], ")"}]}]}
      },
      GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.7]}, 
          Offset[0.27999999999999997`]}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}}], "\[NoBreak]", ")"}],
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]], ",", 
   TagBox[
    RowBox[{"(", "\[NoBreak]", GridBox[{
       {
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "uC"}], "+", "uE"}], "0"},
       {"0", 
        RowBox[{"p", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "uC", "+", 
           RowBox[{"2", " ", "uE"}]}], ")"}]}]}
      },
      GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.7]}, 
          Offset[0.27999999999999997`]}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}}], "\[NoBreak]", ")"}],
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]], ",", 
   TagBox[
    RowBox[{"(", "\[NoBreak]", GridBox[{
       {
        RowBox[{"1", "-", 
         RowBox[{"2", " ", "uC"}], "-", "uE"}], 
        RowBox[{
         RowBox[{"-", "p"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "uC", "+", 
           RowBox[{"2", " ", "uE"}]}], ")"}]}]},
       {"0", 
        RowBox[{"p", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "uC"}], "+", "uE"}], ")"}]}]}
      },
      GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.7]}, 
          Offset[0.27999999999999997`]}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}}], "\[NoBreak]", ")"}],
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]]}], "}"}]], "Output",
 CellChangeTimes->{3.834512514593965*^9, 3.8345136232722282`*^9, 
  3.84227449703018*^9, 3.8422752034392967`*^9},
 CellLabel->"Out[61]=",ExpressionUUID->"30839153-fd82-4a63-a016-e74593002774"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Compute eigenvalues", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 3.7814729129771643`*^9}, {
  3.78147344349218*^9, 
  3.781473448224465*^9}},ExpressionUUID->"e39c5931-fac6-4efb-beb0-\
441a6e3d99ea"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eigendata", " ", "=", " ", 
  RowBox[{
   RowBox[{"Eigenvalues", "/@", "Jstardata"}], " ", "//", 
   "FullSimplify"}]}]], "Input",
 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
  3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
  3.781476066235634*^9}, {3.8422752128037024`*^9, 3.8422752283022118`*^9}},
 CellLabel->"In[62]:=",ExpressionUUID->"ce0bb349-5cb2-4f0b-83a4-be5de002e619"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"uC", "-", "uE"}], ",", 
     RowBox[{
      RowBox[{"-", "p"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "uC", "+", 
        RowBox[{"2", " ", "uE"}]}], ")"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      RowBox[{"2", " ", "uC"}], "+", "uE"}], ",", 
     RowBox[{"p", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "uC", "+", 
        RowBox[{"2", " ", "uE"}]}], ")"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"1", "-", 
      RowBox[{"2", " ", "uC"}], "-", "uE"}], ",", 
     RowBox[{"p", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "uC"}], "+", "uE"}], ")"}]}]}], "}"}]}], "}"}]], "Output",\

 CellChangeTimes->{3.83451255475968*^9, 3.834513630234558*^9, 
  3.842274499358429*^9, 3.842275228888631*^9},
 CellLabel->"Out[62]=",ExpressionUUID->"56b0a5e8-cffb-4f25-8eea-17001a3fd636"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Compute the trace and determinant", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 3.7814729129771643`*^9}, {
  3.78147344349218*^9, 3.781473448224465*^9}, {3.78147830643638*^9, 
  3.7814783120648193`*^9}, {3.781552785807474*^9, 3.7815527960775223`*^9}, {
  3.834579050006682*^9, 3.834579059173544*^9}, {3.842274512058419*^9, 
  3.84227451225204*^9}},ExpressionUUID->"5baf18f4-4931-481f-ad99-\
f5aaddf7e5b1"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"trdata", " ", "=", " ", 
  RowBox[{
   RowBox[{"Tr", "/@", "Jstardata"}], " ", "//", "FullSimplify"}]}]], "Input",\

 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
  3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
  3.781476066235634*^9}, {3.781478323943038*^9, 3.781478329366476*^9}, {
  3.8422752360215187`*^9, 3.842275239132975*^9}},
 CellLabel->"In[63]:=",ExpressionUUID->"c28a2f38-4007-43cf-b4c9-f755f30d66d4"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"uC", "-", "uE", "-", 
    RowBox[{"p", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "uC", "+", 
       RowBox[{"2", " ", "uE"}]}], ")"}]}]}], ",", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    RowBox[{"2", " ", "uC"}], "+", "uE", "+", 
    RowBox[{"p", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "uC", "+", 
       RowBox[{"2", " ", "uE"}]}], ")"}]}]}], ",", 
   RowBox[{"1", "-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"2", "+", "p"}], ")"}], " ", "uC"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "p"}], ")"}], " ", "uE"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.834512576251791*^9, 3.842274515626686*^9, 
  3.842275239779957*^9},
 CellLabel->"Out[63]=",ExpressionUUID->"4f37fb08-8ea6-403d-b301-89ecf38f4878"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"detdata", " ", "=", " ", 
  RowBox[{
   RowBox[{"Det", "/@", "Jstardata"}], " ", "//", "FullSimplify"}]}]], "Input",\

 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
  3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
  3.781476066235634*^9}, {3.781478323943038*^9, 3.781478329366476*^9}, {
  3.7814783600281353`*^9, 3.7814783639524193`*^9}, {3.842275241532503*^9, 
  3.8422752443353252`*^9}},
 CellLabel->"In[64]:=",ExpressionUUID->"5517b3df-8a28-4367-999e-17e479c8c77a"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "p"}], " ", 
    RowBox[{"(", 
     RowBox[{"uC", "-", "uE"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "uC", "+", 
      RowBox[{"2", " ", "uE"}]}], ")"}]}], ",", 
   RowBox[{"p", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "uC", "+", 
      RowBox[{"2", " ", "uE"}]}], ")"}]}], ",", 
   RowBox[{"p", " ", 
    RowBox[{"(", 
     RowBox[{"uC", "-", "uE"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      RowBox[{"2", " ", "uC"}], "+", "uE"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{3.834512577623479*^9, 3.842274517298068*^9, 
  3.842275244997375*^9},
 CellLabel->"Out[64]=",ExpressionUUID->"fb871d7b-af88-4454-9324-07fd9b1edabf"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Visualize eigenvalues and stability", "Section",
 CellChangeTimes->{{3.7814721374536943`*^9, 3.7814721459008713`*^9}, {
  3.7814722240590878`*^9, 3.7814722277070637`*^9}, {3.781472402298065*^9, 
  3.781472410204132*^9}, {3.781472903279071*^9, 3.7814729129771643`*^9}, {
  3.78147344349218*^9, 3.781473448224465*^9}, {3.7815141129222813`*^9, 
  3.781514130440056*^9}},ExpressionUUID->"1badad71-e6c0-4744-b8c1-\
3bf58078cf37"],

Cell[CellGroupData[{

Cell["Equilibrium is stable if Re[\[Lambda]]<0 for all eigenvalues.", "Item",
 CellChangeTimes->{{3.781476265050036*^9, 3.781476288207602*^9}, {
   3.7814763275108747`*^9, 3.781476331218816*^9}, {3.781476402397532*^9, 
   3.7814764323909683`*^9}, {3.7814766701297197`*^9, 3.781476681448688*^9}, {
   3.781477784547036*^9, 3.78147781265799*^9}, {3.7814778624128103`*^9, 
   3.7814778946608353`*^9}, 3.781514150949218*^9, {3.78151422679467*^9, 
   3.781514290002342*^9}, {3.7815151562346563`*^9, 
   3.781515205910882*^9}},ExpressionUUID->"75a5b281-a86d-4e96-91f2-\
c41dbc6d5f85"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"eigen1data", " ", "=", " ", 
   RowBox[{
    RowBox[{"eigendata", "[", 
     RowBox[{"[", "1", "]"}], "]"}], " ", "/.", 
    RowBox[{"{", 
     RowBox[{"p", "\[Rule]", "0.1"}], "}"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"r1data", " ", "=", " ", 
  RowBox[{"RegionPlot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen1data", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "0"}], " ", "&&", 
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen1data", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "<", "0"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"uC", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"uE", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"uC", ",", "uE"}], "}"}]}], ",", " ", 
    RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",", 
    RowBox[{"BoundaryStyle", "\[Rule]", "Blue"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.78147591960296*^9, 3.781475926646345*^9}, {
   3.781476027151567*^9, 3.781476028652916*^9}, {3.781476062425387*^9, 
   3.781476066235634*^9}, 3.7815141879058037`*^9, {3.781514299900753*^9, 
   3.7815143389630337`*^9}, {3.7815147790959682`*^9, 3.781514886796404*^9}, {
   3.7815149248973227`*^9, 3.78151495836263*^9}, {3.781515003778427*^9, 
   3.781515016404488*^9}, {3.7815150850408707`*^9, 3.7815150882864017`*^9}, {
   3.781515314103553*^9, 3.781515321347374*^9}, {3.834514231169278*^9, 
   3.83451430577695*^9}, {3.8345143807965803`*^9, 3.834514381094407*^9}, {
   3.8422745397564497`*^9, 3.842274541424534*^9}, {3.842274702075144*^9, 
   3.8422747102614107`*^9}, {3.84227526770434*^9, 3.842275306898538*^9}, {
   3.842275341996676*^9, 3.8422753421318398`*^9}, {3.842275411247903*^9, 
   3.8422754119265957`*^9}},ExpressionUUID->"033cc442-147b-4f00-aba2-\
3b966159717d"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJxtmHlUU3cexSNFW6rBjlrFgNatNlLcKKHqjH6RIFZlNDhKFSH2aBdqUJRR
adBOOSpxqVqLQa1FS9QyroQqrQUjLkXF4IJVCARkS14IgbAEXGthHufM/cVD
6j/3fE4k7977fu/7vjB8Wdz8T9wEAkFqD4GgS7nHJ2RqQxF9/2Njm0RtoKC5
c0ZwkgeMo/fcieQk+fT4xa5rKoeR+tPfOlQOPWP8PBg/D8b3gPMLR4Vvc5yk
8AXj9xtklTT3t0FxWx25jHE9MK4HxvXAuB4Y1wPjuuDXTp3y1EyPnVayMCZA
bahmfsDwA4YfMPyA4QcMP2D4AcMPGL7A8Lelx+RDOlUt8weGPzD8geEPDH9g
+APDHxj+wPAHhk8w/Pqf8ejFSczMLxh+wfALhl8w/ILhFwy/YPgFwy8YfsHw
DYb/6kXGVUKOY/7B8A+GfzD8g+EfDP9g+AfDPxj+wfAPhn8wcoCRZ7f7KYNE
XcfygJEHjDxg5AEjDxh5wMgDRh4w8oCRB4w8YOQBIxcY+f6RtSFILq1n+cDI
B0Y+MPKBkQ+MfGDkAyMfGPnAyAdGPjDygZEPjJxg5LUtCTuhcthYXjDygpEX
jLxg5AUjLxh5wcgLRl4w8oKRF4y8YOQFIy8YucHIf+DVIf20mkaWH4z8YOQH
Iz8Y+cHID0Z+MPKDkR+M/GDkByM/GPnByA9GfjB6AKOP0HP2RIOsifUBRh9g
9AFGH2D0AUYfYPQBRh9g9AFGH2D0AUYfYPQBRh9g9AFGH2D0AkY/rfI8U2dn
M+sHjH7A6AeMfsDoB4x+wOgHjH7A6AeMfsDoB4x+wOgHjH7A6AeMfsDoB4ye
wKv2nb1rkN2hD8IvRsmlZdT0zPtrIVfM2G3y2O1qw2U6dlw+TqupYP3gcyg+
11ZMe3OH4yB1dHb9q2J943MoPsf9AkPPnCqKEWtrmC75UHNPJjcx9XglfopS
yDE9nxl8VKOzMP0ksn8fvcLKtH8v81qHyMb0yk/ZD0X6BqZx0cmhUqWdqY9H
hFYhbmba/f6jJyjOAbR7PzgX0O594JxAu/eBcwPt3g/OEbR7XzhX0O794ZxB
u/eJcwft3i/OIXvPdusb5xLavX+cU2j3+4FzC8V9GXTsXIRSWEz587Zc0uhK
mMb/sUCsV5Qyfeu/b3/rEBmZ3p7/+JlIX840seP6MqnyIVPxyf2FCnEVU5tv
gEivqGZ6NSltrlxaw/RAiftmh6iW6Wq/ledVDqeGbipuEOlNTIeWTh2m1ZiZ
to/NWCBVckxvbfbcbpBZmB4pW39RIa5jmji+qrWz06my5Jmj1QYrU3G5NlKs
rWfaMcHrG53KxrRElfSbTN7A9HSF9YlZ0sh0i3+4n1JoZxq5LecjIefUiZXD
UzW6Jqa9AnbclKibmeJ5STmSmKNT/c6eAzCek1/r4vwC1TfYcwBm790BxYWc
JJs9B2A8J5/uDUqVS2+y+9qmWH8oUH2e3UecE/j6/GmtXSa/z84J5jv88cPI
RyAoYucG7wf4yzL/2apTXWXnBu8X+PMcUnvGIjnGzg3mBb4f18f3Y/52n69g
5MP1MI+7z1cw8uP6+H70edlW86lWc8vl+9DnCVFOUKksj/18UUWaQv4Xfx8o
rfLMMJiLKcTd+6hhj3Pu+Udl71H+xd8L/N/PShDz+fH/HKF5Ier9Bjq8RLdW
3Oyci15/XBCXtbj+fWD37vBZSr4P/L+5h5cO5B4Z6fnZ6FDlPytd5ufxRwKr
ZEEl+SZ39uaKq13m56bR96KkZdVUOqForTC51mV+Tvhdc02xrZZUFekPJQFm
l/lZ+WX8OPX7ZgrYtiZUbnKdnzvHSPfrLBzVvhesVaXUuczPKcX9BVxqHe2p
6uelnV7vMj+tSeYY4Yx6mva1KcnQ4jo/9/n9fE/SbqPGwOz6zh8aXeZnSGny
FPnRRjpYu2W+eF6Ty/x0bI44qprP8+6FF2QdzufV+PlIuZK/33u8Pu3h3fUe
W39lnZa/nw1xr95X8L2vWt6iy2rl980bx4/p+J5Xuns9V/nfpg9ea2mVZJdS
8Jhzy+SmYjqo33JOZykmt7lvf7Y1/hIdWfrtb6p+FWR9/cl73i/tD2lJCwK0
08uoffjoZOEMfo8IWxHziiCFOn7xTzWsrqKnkSajttW5Pzw1POopSK8gmfJG
hPxoBdsf2uL5txh//ud33f6OKlrQdZsF1bRImvZAlllDbgN6/tnJvyd7NcdO
VUabKOtS8SfipSbKPjg1Q9OHI7ki447Mk6NloZ599Rcs1GdQwvvKixZ6w1GZ
4FhhpdyrM9M1sVbKO6StFvHvs5hVXh56bxvFzkqaJb3ZQANF9Wsc/HtM9Eh2
VvGFnfKv5RhFiXa6nj7cW/1OM8Wv2SGVjnHuW1E1M6Vq/vn3DZg2VfvSPrFv
0ewp2ppil/0i9Ht7ucJazJ7XXyfVD+BullOLtqG3IKjMZd9Id+w+owstY/tY
vz07LBL+eQ57dddD2eEKl/3jxawJXwl/rKCMLnxR5bJ/nOxaM3pUk1v5Z0HK
MzUu+8eijZO/NWTVUJZq0glNlMll/+glfv15J9+33N+jn74357J/ZN83LhP3
5ahPZVmiI9fisn8s++pUoSzPQrnbT5pEK6wu+8cb724MUK60UoxkQ5h0sM1l
/8grCTuk8bHRwJo5PysKGlz2j9hNQ3rpCxsof6fPUHWC3WX/EI1rWuXg7+fq
SXaVbrRzLxQvHmrS8ed5an7gCinfe8uDzYWqFP73pfULfbhUAymfTJzcxs+v
77xL/67me57t9/YhQ4uRjtx+bJe0G2lCP/Fdj3fL2fnN2Hx9lXheJb0YdeCS
an4lJdnfS9eVVJNfXlqBorSaxn3n7sHx88r4+cqJ6q21VBFSvEbIz6Ntb047
qAs0046WH40SM0eBVzxf4fg5NCltvVS+t47MsVUKIT9vLDNnnlYF11OKV9YD
SUg9qdsGvalttVFQftJUeZuNgn+wfmlIb6SmuPAM1ZFGapmdY+nk50ea94i+
2vAmyipcKJHvbaaotLxoYarzvXrWva3VnOp832dKEowK/v2K9/nru7xO60TO
9/nH5pwvhfz7uWXzBKn2pXmTsvjMXPWTe3Qw6q3twp4Gij9ZTt6BzvlzfmFx
BLfuOrWfuxwpX26kj2f/2zcwu4DS3S4tkd8vI1HAYm5pm3M+7Wv5YVjIs59I
1nuZn/bKQ3KUld/g1v1Czxd/NFY70bkfFJTlGUS+zRS5fLhGleLMhf0A/mz/
2virbuB9Cjk1ZIdwXQnbD+CP/T7+/+uLnv1827D6LjkimhySUaUkuZI/Xs8/
76bY1AMO/v05p/zeAEGgk72+8t7vyC2h1cPsdbK8Esr8fsJW9ZMrNDfzm3zV
pnLqGOU+OPyl83bdt7FTtKKUht66KFOuLKW1JboY6WAj6RN252h8jNTn2Wm3
Gc80dLznxH2G6koa+c2fx7xfOo8jNu26pygopy9GLh2hLyynO2PlU9QJD2n0
3fFfOxIf0saycUd1o6voQaKgXeRbReqrncmOFdUUtDcmUxNbTbFeMbmq4Bry
qi+qlYTUUPCqoiYRv/82TZs8SO9dS6L8SSO1rbV0Xa2ZI2+rpZbBmg+lN02U
ZvNIcvB7cEGcx05DupnWBsVnq46Y6fC1NZcVX3A0e5/RKkrkaL23sb1znoVG
NAYP0YZbKGxN8Bj1O3X0dPqpcOmYOhp142S0mN9/i/b3VxkEVnru0z9FV2Kl
DPuGXEWplX6P33BdlllP/5Gamzr5Pfh4gem5OdlGC74LG6neaqOkoWHjldEN
5Nec/aF4aQNFrM1eLuT3X7cZQ3bqAhtpnN7ngKaPnYwHky/LPO3kPiz5lsRs
p6wWe7uZ34Mr1tkF+gtNpAqNGKO82ET/A6HeXmk=
    "], {{
      {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1], EdgeForm[None], 
       GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJxNmHecV9URxe9dii4gK0WK1EVUiiViUMounRVYQIQliAY0FGm6BGFRpISl
CAgq1d4rIcaSGIMlhhBjbzFGjbGXGGs0FsgKknM+9/s+b/84O++WmTkz976Z
99viKXPHVhaEED6JIdSR/K/wlXCzcIhQT6gvdBPaCx2EvcIe4SPhG+Fr4Vbh
eeEB4RahMXoHCd8J3wp3CC2ERth+Udgp3C78IOwXPhX2CQ2EZkJznhsKhyGt
/71QI3ws/A+enm+JbftvhSwSivHt9XZCU+y3FZowbl/LZ0f8WacTerZ3mjBI
GCwcwVxroQu2bPsl4SFhu9BGOBQfLwuPCDuEo4TDWT8aaf2u2Mhy3Rwe3Rmb
11DhZKGX8GP07f8Y1h3rcfA2x1eER4VfCccz11noiV5X7HXDTwn7fiScxHo3
/HXHT2/ksUIfpH32RVq/h3AksZ6INNfXhF3CXUJ/4QT2noLdPuS3J/4HEafH
Q5gz3zI4WWcYen3hUkz8w5lzTKPxNwD0gNdApH2MYG+p8LqwW/i1MEroh/6p
6FvvWmGtsE4Yw5z5XiOsYW0idsuFKmGaMF2o4Cwdx0+I33GMZ87jN4THhLuF
SfBwHFOEscI44afCSNYms26OZyHN62yk7+/PkNafgE/n6XSkuU7FdgV8K+B1
Hv7MZSY6jm8GtjyexdwZwrnCmehUojcZ20OIfy5z5juf2KaSo/Hk5hykffyc
vY5pHtIxnY+0/gJsTCPOwfhbiF3b2ygsFZYJvyA2c7yAdce0iDgd0zvCk8K9
wkXMzRaWo+c4VsDPvHz+VfisZt1rK1k331VIx70aae4XI6vgOIdcLkOa61ps
m+865IXCZuIxr/eEp4X7hA3wXixcLizB9qZaOdiCXjW2ZhD/VuYc31VwNcdL
2Oc8rUfaxzb2OtYrkI71SqT1r8aG8/SW8Lhwj3Ad8dj2BDj4PNxz3FvcV1w7
XU9cX27At+P7QHhW+I1wI3OXCneG1KMKQ+ob7hnuJe5L7jmuz+5j7ifuN+4n
7iON8NeY8e08u6/9MqSe5vp+W0g9s25I/fMjbLinuZ81xa7reTPGrvNtQ96r
st5l6b6znfki+HlvG/Tbo5/1j6Ow3YE167nOdwp572xOPC0YGx3R8Vwx45bo
ebwDO62Jsw25N5cu+HbddM0biG/3iq7ouQ8ch4/urNmua3TWr45hLeslrudZ
r+oc8v7Ro5a/E/Hfk7H9ncTYPk5mbLu9GHeCS+a7N2tZ37LvcfoO2io8g63B
xOZe5b7gnuh+5T7QD04DGGe9ZADcS7DbDn5d4DQEu+4rrrWl8BvKmnm59pfB
aRhj2xrOuATd4cRQhv71Ib036+HnnjASTu4Fo7HlWuramtXNecTpvjAGHq7R
We0cy5p9T0Df/l3nTw95v5kc8l7ifXezpxwe7gNnwmsS49HoTSKGCnyXwbEi
5D3srJD3sLOx7To8k72u7Vm/spwa8v4xHV7e65pyD7pncG5jiMP8XGPn4M81
ey45mILvcdieQgznoeMYKhlPRq8y5L3KtvoT8yj053MGrpsbOLusx8zGzwL2
OQbX/CpicO3N+tUFrJm36/oSeLiuLwt5v1kZ8l6yFH7L2We+1YzNdQXjeeit
IH+L8D0LnotC3sNWhbyHrYaH+81l6Ph+Zv3Kcm3I+4fXsl61GH4b0Z9G3Avg
vok1c3K/2cr6GnwvxPYaYtuMjmPbwrgavS0h71XbQt6rriCGqxjb7tWMsx7j
/vIUvBfj09+F1xCP38msp13H2oGQ+kX2+8p9zT3FNdY9xfX+d8z791mMqX+5
h3weUr9xr/l9SD3Juo15LiKXN+Dbv0X+FlIN3wBnrz0YUo/x/j+FVI9ci0bw
7Hf2DyHVYtfh43l2DXRt38H8n0N6z73fd+Fe8vAEz77DrfBlDg+j6+eJ6Pod
/EtI76SfZ/Fs3S9C6t+F8L+PnN/IGdyEnyfx6xzsJK5n2OP9Hwq/Zf9z5Ne5
/VK4n+e/kpcicvYwPF3/d5Gf7eS0NXE8Qh7+Tr46k6dHydWrwh959nk+wHm9
wPn5+R+cQSm5300+HevT8Hc+Hicnb5Iv5+ptcu35d0O6i87D+8Rv3RJiKCUf
z5IH+3+RnP2TsyznTB7D/v3kq4CcPV8rjhfgv5PcFRHrXcTrXL5EzpzLl8mZ
8/QKedqN73Jy9iq6u8hLKTZfw66/r/x/igjfN+DrfLxVK09vkxPn5p2Q35F3
yY/z9F6tHL9PrnwOr3MOzt8HtXL2IXn7d0jfeObwSUjfjf7/yV7WnKsGzPv7
67OQvjMPDul/F9bx+9uE+abk+01y/gM6hfj5Cl/1GNdn/Qv2mNd/Qrrf9v8l
HG5Gtx4+P8eva6Lrk2uZ76/vru+tv4G/IZaPGdfB5l587+G5AO7fEtcdPDch
Zn9j+3uzIKb6tQ8O32N7ruZmxhSHgzuAjtdqsOH5GjjXoFsXP9/hN+OU1dE9
cHNO7bMZeTiA7mfk91P8HSA/NfjL5vzs/a675l/I89ec7T7m/Xyu5IyY4rPd
/fh37Ps5o/34zZ4LyMlsyXNiyrHtFeDLv238fenvzY6a2xbTHVwpOSimGtRX
8qKY3sv5krNiuidV+jOH/+ltlmwX030fKLkipveykWR1TO/uwZJLY3qnF0pW
xpTb3pIXxlQD6kou4vm0mOzaZlvJTTG9W4dLXh7T+zdU8uKY3vVDJVfHdNe6
Sy6PqaY3jikexzJYclVMdaK95JaY3kvHV48Yu0guian/2ecY/NrnaPxabyy6
l0qOjOndst5B6DqOPsTiuEuJ3TmuIM/rJUfEVJ/WSQ6LqSbVkexEnhvqzwRi
OVqyhD7tWIcQr3M8gDz3k1wWU91dI1kWU81zPhpiZ5TkZTHVs1MlN8ZUw8ol
N8RU82yjEDsNJDtzXq1iitnxdpPsH9N7ab2W6NaXPIK7YT+t8eV7cQh34xLJ
w2KqhT6TIs7F59mEM10r2Sym2jxP8vyYatB0yWLhXyHdnbrk2T6r8Hus5MSY
fkMfKdkrpnfUsbYh3uYx5d05d3wLibFFTGfjc2kaUx6dQ/9u6xBTzR4eUwzm
77Oaz3mdEhNvc/ZdGM99+D/TwHpN
          "]], PolygonBox[CompressedData["
1:eJwllMVuFVAURc+hQArB3WWAE9zdWqwtUNyhuLVASBgw4FsIBMIYd7fi7i4h
BAshjFk7e7DSlZeX3nP32fd1rqgqr6wREcVQE3pkxBg4ht+FpXgzOI8/g1H4
LviL74Ud+HA4jP+GIfhsuIl/hC14d9iH/4AleBM4iz+BuXgX2BP+ziK8AZzE
H8BKvB1cx9/BMrwFXMRfwAi8K9wOz7AcbwWX8VewAK8Hx/F7sBhvBKfxR1CB
d4Bq/AMMw3fCn/BMK/A2cBV/A2Px3fAP3w/z8EI4it+BVXgnuIV/gtX4GviC
H4LydEbKRmcOxmfBDfw91E9noLufgJbpO+kul6AEHwivw3dsnZ5Rs12BcXhv
uB++8wS8DzwMZ1onnYmy0I4bpjNRFqdge3pmzfoLGqd3pl2dgTJ8ELwNZ7IJ
3wzf8AJom96ZdnUN2qczVra6Y8d0RspGHdmAb4SvOH9iZjoTZaGdr8e3wk/8
O8zBh8LncMbr0h1TtzRDJV6V/m5huBsj4Ui4I9PwAfAy3KEpeH94Hu74fLxb
utsHwmf3SrtmWIj3TO/+IPTFi9LZPA5nMyN9d2Wks6and6cztbvS9K60w374
ZDiHPw3PMhUuhGdqnu68uq7PaqU7p67pTupi3fRd1clJ6Zk0izpenD5D/1tv
rmn6Test60x1Y2J69+rI6PRvgN6+Or02/cb1trUjZb8t3Q3tQL8FtdO712+C
3k5Bepd6Q+ri+HSX1cn/piCTNQ==
          "]]}]}, {}, {}, {}, {}}, 
     {RGBColor[0, 0, 1], AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwt0ldbDgAAhuGvyyxF02ioNEhJSdJQhNIQCtFCSlKKkoyMaEi0EyFOHftP
fke3i4P7+H0PnuT2/vq+oEAg8JtfhBNLGrmUUkszPYwwwyxzzLPAIkss84kV
PvOFVb7yje+s8YOf/3f/kOxEJx0M8JACLrNMErdJZIkGjlHPIntpJ4EFLpHP
ReaJ5xZxzHGBo9QxSyw32cMM58mjlo/s5ga7+EANR6hmmp20EcN7qsjlHFNE
00oU76gkhwomiaSFCCY4y2HOME44zexgjNNkU85bttNEGG84RRYnGSWUawTz
klIyOMELtnKVTTyjmHSKeMpGnlDII+5wj24G2ccGhjlOGlfYwggHuE8XvfQw
RAqbeU4J+2lkG6/J5AF36edvrI9JJYRXlHGQ6xwK+tfxOqkfPwo=
       "]]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, Automatic},
  BaseStyle->{FontFamily -> "Arial"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox["uE", TraditionalForm], None}, {
     FormBox["uC", TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 1}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.834514278930832*^9, 3.834514306303125*^9}, 
   3.8345143962456903`*^9, {3.842274520323626*^9, 3.842274542028345*^9}, 
   3.842274710830452*^9, {3.8422753010596867`*^9, 3.842275307445259*^9}, 
   3.842275342511567*^9},
 Background->None,
 CellLabel->"Out[70]=",ExpressionUUID->"e58b5636-c999-4422-9be0-5f425259cbbd"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"eigen2data", " ", "=", " ", 
   RowBox[{
    RowBox[{"eigendata", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"p", "\[Rule]", "0.1"}], "}"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"r2data", " ", "=", " ", 
  RowBox[{"RegionPlot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen2data", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "0"}], " ", "&&", 
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen2data", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "<", "0"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"uC", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"uE", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"uC", ",", "uE"}], "}"}]}], ",", " ", 
    RowBox[{"PlotStyle", "\[Rule]", "Orange"}], ",", 
    RowBox[{"BoundaryStyle", "\[Rule]", "Orange"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.834514372808457*^9, 3.834514390636292*^9}, {
  3.842274544405828*^9, 3.8422745456513*^9}, {3.84227463328498*^9, 
  3.8422746401381397`*^9}, {3.842274687781416*^9, 3.842274692152658*^9}, {
  3.842275375165738*^9, 
  3.84227541436786*^9}},ExpressionUUID->"ce31771c-ea40-40a8-a381-\
a1f411569b3c"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJx1VX9M1GUYP0kKKDAsY4N+LIR2DgRSj5IpDwj+GPLj+GFrCrfFilUQ1xhZ
sLXZ0ENJGtoJRgJiYALqKTALPAMUDzsCDuX4Hr+R+maJCZweASElW5+X7d2+
f9yefd573uf5PJ/neZ/vq8nquPfsZDLZ6OPfonWoqXEpD00Lhm1t94o9aK0m
4OhrbupcayPDSQWdu0VFK8PPkeuCxmpkWJyuUmoFE8Mh0Ts9RUUPw99W3nug
0AoMT8/nX9dY+xmOTfAvEpTDJMWLt+AHDH7A4AcMfsDgBwx+wOAnpQvPg7fg
Aww+wOADDD7A4AMMPlJ94XXgefAWfIDBBxh8gMFHag543fm6+by8RX5g5AdG
fqm54/vM68zXyeflLfJLzTU/R3wfeR35uvg8vOXfDT+n/JzwfeJ1498dP+f8
nPHvsuKMyk9XPki6weBVedZiho8/b24XFfW0I/ZKoiqsj+w2rj2kFZpZHGD8
/+MdtW+gto3uz3p86SyaKb2wtktQdjKewPj/6KnsBr3mpuQ5sFQeYCnewKiT
94NOiAOdkZefG6m6cE9KJ+ThdQav7AVDcljWELm8NHbud0UFvfK99xGrez9d
+O3RlF5zlVpj9jeV63vpcdNelMlMzOIcfriHOIgrry5qT5WP0P6FjJJ8ax51
xE3PuhsHqMq9IcSi/Iky/kmQG1Mt1Hz3doqu/Bdyq6h7K8vZzCzO4Yd7iIO4
yIP5AubPwWvzlfVzWuEC4+20qeyas3id1aWPXHlArutmFufwwz3EQVy2N/73
R3/gj37AH/PGW7zXD2bG/lKqbjE9+DkHxjykfB1yTBX2M9OLf0fAmI8HqXtL
ArU/MD2l/IDhh3tSPIDBA7xQF+pEXahzrjZpW1bUMNOpdI8+Uz7Rz3QMX+7x
nVAgMJ2/2mX5V9FjYnkRv/DtiCDdbTOdj/JpkPu3sjqwvyZ140/LQvrIsaXI
eEyoZu8B+yryqfwhZekgnc4xpMtjhill27DzCrGaIny9S4TJfnpyYt4xd10r
TfbktGuOChT0WqhcNJuYxTn8cA9xEPdM9IJM9vEIbW9b7Z1nPUKRfyf661wH
qfps/Atr1E2MN+rYe9GrSm8YYO/cWnY5SVVnoavrdd1CQAfTAbq4784RUt+3
sL1RuMP9sPMuM7O4hzh/nM3q1qzrYX0FD/B6d/6jROO4kc0JeKOOyc3dPi5i
I9s3eHenF+H8CNMZupscDNnO3iNsL0F39AFxEBd9Qx/BA7zQd8wB6kBd0AG6
QEfoCh7gFbdId2GEYgp8W1ym3WlGsNnLTg7SI/vLz4bPNtKJfQkbdKF99OHo
WLRoM9KWNXXJql/NdO/l+gjtGz3M4hx+uIc4iIs8yAvdoKNnks1V/GRo6Tsz
5fONEDBExW86fZqbcYlu2CeEaN362fuLKM5r09zso4YsJ7vc2huUdr42Payh
d+k7t+VPmyKzl7Tj273EVbfYOfwQB3GRFzwwx56xqdplshqaO7X2uNA8RMEx
nV2i7SLNbN0ToPtigJ5ob1OJtha2P+KnD90RTM3sXdS/frImbGMrle48bNAk
99MzB9L8AjMNdOJ+g0q12kJqeVSlXN3F9klT/Ao/j4hO9s4K2lZmhfWZKPzh
3YeK5QJtmvqsUT/dzc7hh3P4oW7oAB7gBR7gBR2hK+qGDqgbOqAv6FNcyfwy
2egwJV1yfOegtZx9n4LavBwq9cVsr2HPYe9hD2IvYk9CF+ik/XxsVrFhad6g
O/qgL6Fz8oml+QMP8Crbl9lhHF+aR1ic4x7iIA/y/gfzqUZS
    "], {{
      {RGBColor[1, 0.5, 0], AbsoluteThickness[1.6], Opacity[1], EdgeForm[
       None], GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJwNw1dTEgAAAGD3ZCMb2bJlylKG7D1VphNnuRUEQRDEh67rUcuyeXVd156X
mS0rW5alPXief6B/kd/dR45O+qIlRUVF/06WngQUl4BKy4AlpeCycmhFJaS8
AlZZ5TPU4KAwFBCEBoExYAgWAkXUAuBV1RAsJpuBExFIEhJFRqEpaAwVg6Vh
cWYTHgkAMghEJpHEIpHZZAqHQs3naLwGOqEOwaczBAymkMmyWbmSRp6Ux5fx
BdoWpZjNSal0coGwEBOrJFJlk0QtlfkdioRcoxCK6PWEoEtdi6grAQJKQcAy
MKgcAq6EQYsBtWGPtgIKEbHYACQCiEKC0CgwBm2AV1fBYU6/whHQLGaNCCKh
jlCPJBHRFDKKTMJQKSYaLsehEpkMEotJZrMoHHY9vQFLo2ZwsKW8hc7nMQR8
plDAEgnZYpFVzGmUNBEYdL5MKpDLhAp5TCFaLtiblC0SlbKB1yhVqxJqmUMj
T2mVLp1qMuj26Fu5Vk6LtrW7wzoV9na2GfAmXJsR3tVuarfQ+s34GkM1FIcd
sHGjvpqegCPuVKhdqt6QS+tpVTjkHXZx0q2e8DsNbXBTO83aIZ7u8ae92r6I
hyeVzPQGH68NrD8dbBs1dozbbd1WX6fBE9F7w1p/jyOZ0Ji7TLN9YVdIl0iq
55Oa5LxbHBNNzYbjseYru865ePPIKf3kTDAWV6hSyolpvzwhi88528cswV5X
ZEQfGtbBMzBajorL4I1ZeCaLzy7iOTmuvSDOL3Fzea4lTyssN8cKzZd2bNvn
zTsXbemULpXWhvs8gSHN6WjNq03LQlp3/CDyesv+csOYXvAe3Q8d3gusfRlw
B9VP1gfPHUfefxg9exS6tufee+veuG389Hls/533z++pM4eBy8+Gbh34fv2c
/PF94vqLkdWP/Q9Xohe2zbtvnDf2vQc3fY9W+68+H175G926a//6bXzzjuU/
dnTQtQ==
          "]], 
         PolygonBox[{{150, 143, 93, 71}, {201, 187, 75, 109}, {147, 140, 91, 
          97}, {199, 154, 74, 100}, {171, 156, 102, 87}, {195, 151, 71, 99}, {
          189, 203, 110, 77}, {180, 141, 66, 103}, {168, 223, 130, 72}, {213, 
          145, 95, 113}, {214, 146, 96, 114}, {215, 147, 97, 115}, {153, 144, 
          94, 74}, {163, 142, 92, 82}, {178, 137, 58, 101}, {173, 160, 104, 
          91}, {166, 222, 127, 69}, {170, 224, 132, 75}, {179, 139, 63, 
          102}, {182, 174, 98, 106}, {172, 158, 103, 89}, {212, 164, 82, 
          98}, {145, 136, 87, 95}, {191, 209, 112, 81}, {221, 161, 81, 124}, {
          222, 165, 28, 127}, {217, 157, 77, 118}, {219, 159, 79, 121}, {205, 
          214, 114, 111}, {133, 192, 107, 83}, {223, 167, 44, 130}, {135, 200,
           109, 85}, {134, 196, 108, 84}, {137, 217, 118, 58}, {193, 183, 69, 
          107}, {190, 206, 111, 79}, {202, 213, 113, 110}, {181, 162, 105, 
          104}, {142, 133, 83, 92}, {141, 221, 124, 66}, {139, 219, 121, 
          63}, {146, 138, 89, 96}, {208, 215, 115, 112}, {165, 175, 99, 28}, {
          169, 177, 101, 58}, {143, 134, 84, 93}, {167, 176, 100, 44}, {144, 
          135, 85, 94}, {197, 185, 72, 108}, {224, 169, 58, 
          132}}]}]}, {}, {}, {}, {}}, 
     {RGBColor[1, 0.5, 0], AbsoluteThickness[1.6], 
      LineBox[{2, 1, 11, 20, 29, 37, 45, 52, 59, 64, 67, 105, 162, 181, 160, 
       173, 220, 140, 147, 215, 208, 209, 191, 210, 161, 221, 141, 180, 158, 
       172, 218, 138, 146, 214, 205, 206, 190, 207, 159, 219, 139, 179, 156, 
       171, 216, 136, 145, 213, 202, 203, 189, 204, 157, 217, 137, 178, 155, 
       177, 169, 224, 170, 188, 187, 201, 200, 135, 144, 153, 198, 154, 199, 
       152, 176, 167, 223, 168, 186, 185, 197, 196, 134, 143, 150, 194, 151, 
       195, 149, 175, 165, 222, 166, 184, 183, 193, 192, 133, 142, 163, 211, 
       164, 212, 148, 174, 182, 106, 10, 9, 8, 7, 6, 5, 4, 3, 2}]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, Automatic},
  BaseStyle->{FontFamily -> "Arial"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox["uE", TraditionalForm], None}, {
     FormBox["uC", TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 1}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.834514393223419*^9, {3.8422745231164713`*^9, 3.842274546137476*^9}, 
   3.8422746410222263`*^9, {3.842274692872795*^9, 3.842274713104623*^9}, 
   3.842275389397727*^9},
 Background->None,
 CellLabel->"Out[72]=",ExpressionUUID->"a95764da-9008-4c4b-bca9-377f52eb86b6"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"eigen3data", " ", "=", " ", 
   RowBox[{
    RowBox[{"eigendata", "[", 
     RowBox[{"[", "3", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"p", "\[Rule]", "0.1"}], "}"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"r3data", " ", "=", " ", 
  RowBox[{"RegionPlot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen3data", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "0"}], " ", "&&", 
     RowBox[{
      RowBox[{"Re", "[", 
       RowBox[{"eigen3data", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "<", "0"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"uC", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"uE", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"uC", ",", "uE"}], "}"}]}], ",", " ", 
    RowBox[{"PlotStyle", "\[Rule]", "Green"}], ",", 
    RowBox[{"BoundaryStyle", "\[Rule]", "Green"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.834514489270554*^9, 3.834514494944057*^9}, {
  3.842274549048967*^9, 3.842274550913122*^9}, {3.84227471512286*^9, 
  3.84227472092936*^9}, {3.8422753931434317`*^9, 3.8422754166240683`*^9}},
 CellLabel->"In[83]:=",ExpressionUUID->"07ec3a20-4aa2-4a61-9ffb-2d864700b0e4"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJxdmHlUU9cWxqNFLWpoq61iAOtUGy0qKqHqq24kiE/lYfApVdTYpR2oQVCe
Q4P2laUSZ2sxqKWoxIGnqIQqrRUjTqgYHKsQEuYMEAJhCDhRS95lLb8TFv7j
+n3se863z93nnH0zdHnsvK+683i8E914vI7/ixZE+St1lfT2mTMequnR07Z2
m3xYozAynnDOvadFZGZcudAQw7dYGO91O6MTKWsYf5a1MVAqrmVsWxx6WuGw
MT7Uy6efWlXPOOSCPV4naWDcLM01OZ2NjMPnjzuok5RTXsGI8O2ODIJfMPyC
4RcMv2D4BcMvGH7B8AuGXzD8guE37ObA2G2OHOYXDL9g+AXDLxh+wfALhl8w
/ILhFwy/4Oev99xSOAy0dN+DSIsoj/kHwz8Y/sHwD4Z/MPyD4R8M/2D4B8M/
GP7B8N+f3mtXOLTMPxj+wfAPhn8w/IPhHwz/YPgHwz8Y/sHwD/7lZH2LSKkj
y/PTEqXuEcsHjHzAyAeMfMDIB4x8wMgHjHzAyAeMfMDIB4x8wMgnMGzOMIvo
KcsHjHzAyAeMfMDIB4x8wMgHjHzAyAeMfMDIB4x8uvpHXvAPhn8w/IPhHwz/
YPgHwz8Y/sHwD4Z/MPx39Ys84BcMv2D4BcMvGH7B8AuGXzD8guEXDL9d/cE3
/IHhDwx/YPgDwx8Y/sDwB4Y/MPx19QOf8AOGHzD8gOEHDD9g+AHDDxh+us4P
X5gfjPnBmB+M+cGYH4z5wZi/63zwgfnAmA+M+cCYD4z5wJiv6/iYF+ODMT4Y
44MxPhjjdx0P82A8MMYDYzwwxuv6PMbF82A8D8bzXeMxDuLBiO/6dzyHv3dl
xLU7O/5VkLp02gc7HSmM0X+dOCUdq1aVUvfJY3Yoddfon+FXlkjFeoo5cP6h
TvKAccMrr118SyEhHjrGg37uzKMoobqK6Ys/Vz2WSE1Md38rboqcb2H6xcyg
4ypNNdO/iuzfVyuzMr1/T/Nah8DG9Ou/ZpcJtHVMj12aGCKW25nu7R6hlgkb
mY71gA6f6I+QB+5r5In7AecsdMRDx7kGHeNDx7kCHesCHfscOtYJOvYhdKwb
dOwb6FhH6Kh76FhX6Khj6Fhn6KhT6Fh36Kw/eqPjPUBHXULHe4GOOoWO98T6
8Llbr6o0RTTwxIUIOb+Q4v6aL9TKipn+4f8++skhMDD9/rznrwTaEqbHt99e
LpaXMV2YcbBAJqxgum20v0Arq2T6jYTUMKm4iumHity2OARGpq/2XXVR4TAy
PWRzYZ1Aa2L64OKpQ9QqM9Nbx6TPF8stTL+3xWOHTlLN9GP69Vdkwhqmx4+r
aHY6a5guSZw5UqmzMl1Yoo4UqmuZ3u7n+aNGYWN6kSLhpkRax/SzpdYXZlE9
07dOCPeV8+1Mj9x+6Qu+xc708eVDk1WaBqb39N95V6RsZDrq9ND7hQUWUbbr
nn3DWNcW2frDAcqLrO7/qIn1DVDeYfsCjPf09f7AZKn4LttHScfiL2kUf7J9
BsZ7//al0S6RPmHj43sB7znL/HezRnGDcdfzDf0s6oY7DL15PBfj/EP+qEvk
j+9B1I2Hj/FctegEY5y34KlXJrYpdVnMT+/Pjt7kW26x+TSh/RKF6sfsvMK8
mA/+sd5grDfyxXrAP9YbjPVGvjhX4BfnDhj7APlhvtOCS4HFklw23jVb1ddq
1T0Wv7U97vAex07mH+ct1gPPwy/OX6wPxoM/3FfYhxj/synDzqs6/X5QaUyo
lHT63lXsqvDgdfr+Crlz6oSGi49Z0aTJas6gutheT2Qd5/f66+vUXPw+z6+7
eXHxhm+HS+VcfPvvE5J1q7n7LXRl1Fu8JDq27Kebin7cfRn20Tfb4q66+uol
2fvk3PM3bsX24nWut7ebmkXZxbTKzbNNMeE+q+dHpakyKTf+NzEPROpO9Zai
3XpBU11IQaMuLJeaCql16MhE/gw9pSbM91dP15NEfidCeryUXuqe9eClldL8
jtfFq6R5HWXRztXJ+z3+dnLrtVCc+lSSWUVZVwu/Ei4zUc/G6KnypSaSytIf
SDwslJ0yNV3V10J9B274VH6lmpaHeLyjvVxNOTdmpqmirfSuo3yDY6WVomI8
3bVeNso9rK4UcOf7AEHtGgd3jkfPSpglvltHebcuGQTxdhI8k5yXfWenuDU7
xeJRjXQ7baiX8mNXX+3512WhvimDfHze00s67Y+9e8NnyTvdx6iHl5Emg7o5
hyZ8mrVB2Ol+xv619n4x0Yvj4gqPdJ3ZdV/jnh7tP22quhMfWDh7irqqkK1/
sJvXcd0+7vvqF3uJzFrI4hwhucHKgzoWh/ma1HV9eIF69r6PLNasFTYaKM2x
95wmRM/iwo4sG2B5ZmBxyCe0154yyZFSth5t55eGyP9VTq9n+f3AP1nK4k49
41lF88tZHPZPesefX1ew/Tk60dnHUlhJGR3bo5urf9g88vESsb6SxWH/dC/5
JlB+ror1DcV+j9byE420cNPkn3RZVSzO70/VLdl2I4tDX5GlmHRatcTE+gpF
aVqZyN9MPYW925xcfSGu/Pu4scpPzSwOfYd0gns/bR8L6zv8t68JkZq4Onxi
WC58x9Wf7B4lPqipdsWhL+lbro935FSzvsQ4MUitSKqh5T+cKZDkVrO4KYX9
eZbkGhaHviVnR4ZJwNUz+pZ9Ff081dNr6d1PNvnLV1lZnDXBHMWfUev6HehN
XxMl2hgqHmRjfc20XaYEXRO3L4pCD6u8Xf3PAd/fHotaXXHoewZUzflNll/H
+p76gOxa59F6it7s01NbUMfigosTp0iP17M49EV5u70HKzfYWV+UYtw6Tzi3
gQRjG2Ic3P5DnGNLxHHFvAYWh3tr9SS7QjOykVriuFXm6sSvn/Ch+ycl0372
Kv6Hkqs7+Yvxk1u4/Tc1L2ClmKtj4aLBJk3H71vrF3hbknXU9HRLgSJJR8fu
P7eLWg002/ejw7omA70eceiqYl45pW+5HSOcW06+uan5suJKSrBPTNMUVXLn
6Krxym1GGvuzm7uFq7ftH0xL0QSYqTS4cA2fq4+A6x5vWbj3vLPppEFktpA5
ukLG597fpNT1Yun+GkryzHoqCq6l6pkzzyqCaikwL2GqtMVGypaBH6ibbdQQ
G56uOFZPQUet3+vS6inVa9g76vAGapp9qdrJrQ/yxf3xpfnS93zu/kK/0rbo
izHq8WXk0Jfcsaz7nXrv8TyrEbj6k7TuVxdLn+jpy9n/GR2QnU+Zog0GGdd/
oB8JPuOzk7+uiGz/3vSHZsATdr8I/BdZlrVkkKTPcl/19TI60HR0SPCrXyks
88c8xeYSyvzFb5vyxXV239UPHnJCyt0XuI/iMkrIKyCPWi9ci5SuMNDFBYUR
lnW3yRHR4BCNKCbBq9/u61Y/ZP3Dwk0nbSru/sD91bTFT6zm7peUJR/u4Pfg
+qdF58KULx7TktTcpfzkRsoqWCCS7m+kfH2uTjDa1d+g74hcMVSlSGqkUz3G
H9BVllPfV2e7z3ilYve3W9+iekNTiut7c++Cy5L2Rjrv1tJsTnaNg/oa/uPf
J7w61Vf7CLdB4Vx9maKTDzm481t0PW+cljuvwXNKHr/PC3hKg+9dkchXFdPt
0fVOwcpi+m74smHaghIatnnPY1l+CT2N57UKRlfQJv3Y45qRFeRZ+8goCq6i
aM+oHEVQFd1WquZIW4wkyJs0XN1spLWBcdmKY2bKj3XfrUsz07D6IB91eDWt
9zK0OudW06OD/RU6npVG3MlYKuT67f+KzQ1Orq/+M27jbUlmLfk2Zn8uXFZH
CYNDx8mX1pEhJfGaxMNOY7Xeh1R97aQIiRglv9JApevsPO3lBlo9xF4jyS0i
zx+8Djpyiki7Ye8llbeB1hZposSDDDTy4bhdjvgyejBGOkW5oYwC90dlqqIr
SXnDmehYWUkN0yYP1HoZKSjmUYOA+95ItbknOLjviqZBqs/Fd000+4DBKoi3
0JFba67JvrPQy+lnwsWjaih0TdAo5cc1lG7fmCMrtlKbd/8kTZGV5v8cOly5
zUan8k1t5kQbdZ/hs1sTUE8Ra7NX8Lnvgawme6uZ6/vdhiTeE5nt9H+h3dwe

    "], {{
      {RGBColor[0, 1, 0], AbsoluteThickness[1.6], Opacity[1], EdgeForm[None], 
       GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJxNmHucl0UVxs8sIggLXlBA1BTUFJCrgLC6CKLckZsIosg1TUIR3ZQVVjBA
EFbRBV25uobl3exCRkpERkaWmpqZmppaCWUUFabU2nk83/cz+8f5PDPzzpzn
zJl5z/P+fu1nzB13TYmZ1SezRo5NHQ9xu9/bjRz1rIVbS4u25hwOHuJ2BNjY
7UjwULejwCZuxzNPz08A9fwYt8Pcmrn9ya3UrY1bW9riPbYBfztQ/MeB8tsK
nqZuR4Py+wV49Px01iiWk5gn/kN9j/WOW9z+B//Jbqc0iOVUULF8EVQs3RvM
PREe+e0Ij/baCVQsnUHNPQNULHvcOrh1detGW35Pg0dx94BHsfQEFcuZoOY2
8b0ktwe83Ysx7bs3qLj6gIrrLFBxWQr+s936MqYYz2kQ1yB8yUd/9q+4zgUV
10BiEv8ktxFuI93OY0zry/Gn9efjT7FcAIr/r25lbkPchhGXYhkOluO7nFhG
goplFDjA7UJQcQ3Fp3yNZkxxjWVvimUMY+qPY0xxTYRbnN8grvFuE/CrGC8G
Ne8i5gyFvyd5uYRYFeNKt6+5LXG7lFgV12WgYpkCKsbLQcX1TTimuk2jrViu
Yo3mPsiz2W4z2YPOZBaoWL4ETna7AlQsV4KK5cug/E6HR3udAcrvV+DS84fo
V7hdzTzxXwOKfy4o/mtB8c8DxX8dKP4FPNPYHHjk96twaewGUJw3guKcD4qz
EhTnTaD8Xg+P8rcQHo09DEeV2820xbOUmDR3EWPiXAyK8xawknOuhPMRfC2z
OPub8HUr4/J3O/3lbivwK3+3gfK3ivVLOTfdKd2hR1lbzfpFrL+jgb/VoPzd
Ceo+PgZ3jdsa2lq3FtS6u0Gtuwe8y+1x1tTSX0mM9zIuH+tZr3XrGFP/CeZt
cHvS7T63OreNjGvut+hv5rnam7B1+C7xOrbfogbq/dO7OthCu6Qd0pJCw9Qv
NEz9QsPULzRMfWlLK/qFxhR91fMT6J9IX9olvWhjWcPaWtawYy1rWDv4jqdf
aKT60rBj4PuzhUa2pn9SAz7V6k7wSSNOtqxbp1jWrVPhO42++E6nL76O9Atd
Uc2SNkoXO1jWLfHttdCtLjxTPe8Kd3f64u5Bv9CtHpZ1sSPcveArtKqXZa3q
zXrV9nLLutWHOPrSl+5Jw8qIQxpxtmV96k9M5TxTLAN4pjWq00PgVv0fCPcg
+oU+SQ+KGq36+JGFVg1mDwPxK+5h+C10a5hl3RpuWatGWtaqUfgYTV/cY+gr
lrH0Fcc4+oX+ye8Lblst7r24VasLLRKOt6xVE+DQPiYT00SeyZ90YRIxTaY/
gjHN+7fbH92+bVmrlBfVGNUK1Y6X3J6yrFfTLDRL3NPpF1oyHb8z6RdaVfRV
z+dY1pKrLWvVFZa16kpiuYr+y27bLOuh9GI2fHPoF1o4y7IWFv258BVacqNl
rbrWslbNg+96+q+6PW1ZD6UdFZa1qtAS1eebLWvhXMtaNd+yVlXCsYC+OBbS
f81tu2W9kr8qy1rYUJ8WWdanxZb1qeiv4AzFsYRn4lhK/3W3HZZ1TFqzzLI+
rYBjOc+0fqVljZE2SDPecNtpWa90X6ot65PuUPHdM8WyVq0mjlX4e8vtWcva
Ja2psaxbayxrlbTmbbddlvVK+lLL3HX0NW89/XfdnrPQqEL/1jJPGiTNes9t
t4VGac0mnsmkVdIs/a7Sbyrpwwduz1tonX5zfd0+L1/2mds+i98jjX3ggMVv
MdU7vfOqAd+3eL/1Lt/pcwa4PePtH1i8Z8U34YO0f0hbd1738CHaz9DWndSd
eZj2z8hNLXt8gvaPOfNlnNejtH9Cu5r8P0b7p7SLb4nHaf8IruKb6hHaP4dL
+XyRfWqPvyCvhe4/SfuXtOssfsd9x0If9Tvouxba9Cu371nUZ3FuZ48b8bmB
nD0F16/J41Ryto32K+RxNjl7mvZvyGMFZ7IVrt+yzyr2t4P278ij8vYOua4l
fzsZf5OcVpO/Z2n/npzWkKfnWKu87mL8D+RxA3t6kX29Tx43k7PnaW8lR4MZ
e59x5f4F9qQ87cancvMSOdtGvoo79jI5Uz5eJU/byVEF+XiNc9hBjqrIx+vk
aSc5WkaOXyHPyscb5Ek5eKtBbt4kP7vIUQ17fY/YleO3yZNy8y75202+NpDL
d8inNE13qw78gJw15r8Rvb+fWLyb+iY7LAXq++tDi2+11sxRu7nF3fyQ8U8t
3nF9P+lu/o38/53x9szZS1t3+S8W9/mgxfdt8c0hHIKPjxg7SLufuFPwdMXH
XvzsgasD+/jYos4ohn9Y1J9/wtUPHrXLQD17APwX7Y/x1Rz+g6zbh98tzPkP
XAfIUXPy9An53Mf+2xPvp8S7H5/9mKf1OotmKc6g+E4V6jtUtbSesRL+Gylj
fj1nUc+8LfB+xvzSFLlTUS5lbRf2rzF932mshP9bDnA/FE9LxxYNar2wBeO6
P3KrbybdszqLscb857afWMuITXWtNfdH7S3wNiEerW2KT2lIU/wU97I1+VC7
Ddgs5b03I2+6Y03Yyx7uTRewlLaRk67EUEp+VmqMmtuIvZcQW0vycJvbGdTo
Foxr+WTtKcW7qO9D1RXVq8sVS4qatNjxyBTnfoPjLSnu10zHzvg83PGoFHep
M3waX+64wu2/3h6YQi9Vn7oQt2Iud1ydooYtdVyS4tv2Ase7UtSti6XRKerE
hY51KWrPcMf7UtStwY6bU9S2/o4bU9Swexwnpag9wxxrUtS8UY5rUtS28xw3
paiX/RzvSFFruzmuSlEnlKdLyZVimUg89zuOT1H/FNdoYlMsQ4hHuZxKPts6
Lkjxns9wnJfi7vRwrIZLexrBvo5zXJTiPe/teHsKfVBc3Zm/Qe9ECk1Y73hW
Cv1Z53hmCp3R/gaxR+XmXPLTyfHWFHWhLIUv+RmXYm/aV58UfuWzwrEqRW1Q
vN2IYazj2hT6oLyOJLe6Fx25G7pHR3OXWjm2SVFnZqW4E7oP81PcJ83XvTiH
u6G7cD73QefTlzNSPnqRkyNS+JXP6xwrU9Q2cbaDV+c/lDswxbE2hV7d7XhR
Cr1amMKX/Ex3vDeFFml/Y9jjhBRrPv+GcLwMP5o7jfmXpLh/uns9U5yHzuL/
YYKdTw==
          "]], PolygonBox[CompressedData["
1:eJwlk9dvzmEYht9HnDm2FV9tacQqam9qqxZtrdZerdh7E/+G8Q/Yau+9914l
hGgIwoED15374Equ5Pvye9/3fu4nU15ZUFErpZQPtaF+pLQWvuMZmIv3hCr8
PnTBm8At/DDMxwfCKfwxtMTHw0X8FSzAh8IZ/Cm0xovhKv4OluFL4RP+Gebh
feEE/hDq4TPhAL4XKvCJcAV/C3XxctiD74ZF+Bi4gL+EdvgO+IPXgSX4NLiB
f4CF+Ag4hz+HOXg3OIrfhUq8BK7h76EVXgSX8TfQFi+Dm/hHWIwXwCX8NbTB
p8J1vBpm453D2d2GFvhoOI+/gK7h/+i3I9AM7wPH8QdQGD5D39YdJoczURbK
tAO+Ar7iB6Ejvgpq8EPQP5yxsj0JvcMz1mx1Rr/wmTpLM+gRzkRZHINeeFO4
l9yJQXhzeJTcAXWnAfxI7tAWfDPk4d9gRngGyl6ZbQtnqOx+wxp8dXhW6sS6
cAc0e31zejhTZakZ5oQ7pO6oI1PCM9OslPkGfH24y1+gNNxBdU8znRSeqWap
Tg0Od1pdPg3Dwx1Wd89CI3wj/Ex+k7JdGX6bMm6Izwpnvy95VzqFs9fObMW3
h9+6H7Lx/PC3nyXPbnn4rpqhsh0QnpUyzuBDwnd7At3xLLiT3Nld4TP07X+w
M7wj2o2/MDbcOXVNOzIuvDPaFe3shPAOa3fV4VHhHdFuqKMjw3fWXbUzw8J3
0l2049q19uGztHObwpkoi1/J2TUOuzJUl/LC3VKn9JbccPf1pv/zDKGd
          
          "]]}]}, {}, {}, {}, {}}, 
     {RGBColor[0, 1, 0], AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwt0tlSDQAAgOETZWnGCBVCGyrKvhQp2ihZctBiCVFKG0W0qpRSpEQqpcQL
dOXlXPQ14+Kb+R/gT6hqDjaFBAKBv2wSq73EC1FFFO85yCAH2Mpm3tBOC62k
8YFUHhNNJzvoYDvbiKCLbl7yisOMcIgn7OQtu+ghhki20EsfbbzmKKMc4SPH
2M0JZjhONSf5xCn2kM4sp6khgzHOEEsmPzjLU87xmSziOM8c2dRygXFyiCeP
eXKpI58JCkjgEj+5yDMK+UIRiRSzwGXqucIkV9nLdRa5RgMlfOUG+7jJL4I0
cotv3GY/ZSxRyuoE5UxRQRJ3+c0dmrnHd+6TzAP+UMlzHjLNI1J4xwD9DDPE
P28ts8gC88wyxQSjDNBFK3VUEqSADNKIJ4pwNrKB9awjjFDWsub/1yuenkGm

       "]]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, Automatic},
  BaseStyle->{FontFamily -> "Arial"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox["uE", TraditionalForm], None}, {
     FormBox["uC", TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 1}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.834514495672421*^9, {3.842274525902257*^9, 3.842274551576386*^9}, 
   3.842274722766551*^9, {3.8422754027542152`*^9, 3.842275460007201*^9}},
 Background->None,
 CellLabel->"Out[84]=",ExpressionUUID->"fa5c03e0-9faa-4703-8d71-88426a8e8bb2"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{976, 1067},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{1, Automatic},
Magnification:>1.25 Inherited,
FrontEndVersion->"12.2 for Mac OS X x86 (64-bit) (December 12, 2020)",
StyleDefinitions->FrontEnd`FileName[{"Report"}, "StandardReport.nb", 
  CharacterEncoding -> "UTF-8"],
ExpressionUUID->"1ad00490-54fd-49bc-b44c-d2724c774289"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[422, 15, 227, 4, 121, "Title",ExpressionUUID->"a73c02a0-3198-4cee-a01e-9ddb9794a7ac"],
Cell[652, 21, 309, 7, 73, "Author",ExpressionUUID->"955d9765-5966-49fe-9395-d4bef08eb8c0"],
Cell[964, 30, 283, 7, 38, "Department",ExpressionUUID->"a0f81619-815d-4d6f-89c2-4637a266c928"],
Cell[CellGroupData[{
Cell[1272, 41, 321, 6, 135, "Section",ExpressionUUID->"28badf05-170b-427f-b6c8-cc4a2b50f70a"],
Cell[CellGroupData[{
Cell[1618, 51, 593, 9, 35, "Item",ExpressionUUID->"763b3740-75ea-4318-bd0d-d9abc7667534"],
Cell[CellGroupData[{
Cell[2236, 64, 2277, 59, 156, "Input",ExpressionUUID->"f039ab92-f90a-48ab-a750-f5331c5afff0"],
Cell[4516, 125, 1068, 29, 78, "Output",ExpressionUUID->"7e16a359-969f-4700-82a0-1f622b41ae60"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5621, 159, 297, 5, 51, "Input",ExpressionUUID->"d0cdc25a-626b-4417-b973-3fb5a7bb112a"],
Cell[5921, 166, 243, 3, 51, "Output",ExpressionUUID->"2376c29c-241c-40cd-b9c5-76b123ffc64d"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[6225, 176, 319, 5, 87, "Section",ExpressionUUID->"ca20d7d9-cb7a-458b-a457-4d93f93bd2ac"],
Cell[CellGroupData[{
Cell[6569, 185, 1074, 26, 78, "Input",ExpressionUUID->"55d4307b-2def-480b-909a-8fa64afb3df7"],
Cell[7646, 213, 2116, 41, 615, "Output",ExpressionUUID->"c1841426-9856-4c3e-957d-625ad9a262a7"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[9811, 260, 222, 4, 87, "Section",ExpressionUUID->"792ae558-5cd9-4cb4-9098-f15685b3a8a0"],
Cell[CellGroupData[{
Cell[10058, 268, 558, 13, 51, "Input",ExpressionUUID->"ec9f6ea3-d58f-43b9-bfb0-6e78c6ee5dc3"],
Cell[10619, 283, 593, 17, 51, "Output",ExpressionUUID->"191d46de-c43f-4e33-a5c6-5f23a4a3cb58"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11261, 306, 275, 4, 87, "Section",ExpressionUUID->"0dac2df4-84b9-41e3-93d3-7bfc871b7768"],
Cell[CellGroupData[{
Cell[11561, 314, 489, 10, 51, "Input",ExpressionUUID->"f0335e49-428b-48c5-824c-82cbc6e4009e"],
Cell[12053, 326, 1766, 52, 104, "Output",ExpressionUUID->"7c7bc343-4570-46c0-9caa-fa4da1db0314"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13856, 383, 241, 4, 51, "Input",ExpressionUUID->"22eaa91f-3917-46f7-9d83-baf60c961f48"],
Cell[14100, 389, 1927, 61, 86, "Output",ExpressionUUID->"2be04c0a-2bf7-4bf1-9ccb-1285bc3809c6"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[16076, 456, 330, 5, 87, "Section",ExpressionUUID->"6e8de002-011d-4fd8-8fff-c3e3d9524dcc"],
Cell[CellGroupData[{
Cell[16431, 465, 480, 9, 51, "Input",ExpressionUUID->"f3d8facb-adf3-4ead-bfaf-960ba1c0f366"],
Cell[16914, 476, 1266, 41, 78, "Output",ExpressionUUID->"cbfdf903-3001-4670-8722-a3c13c157685"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18217, 522, 218, 3, 51, "Input",ExpressionUUID->"592d3cf7-6d47-4849-8720-ea76a052f7e4"],
Cell[18438, 527, 2267, 67, 70, "Output",ExpressionUUID->"d67a3d1b-bfa8-4985-948d-677405659664"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[20754, 600, 366, 6, 87, "Section",ExpressionUUID->"5c727f8b-5d37-4f37-a49b-e9bdeb9dc426"],
Cell[CellGroupData[{
Cell[21145, 610, 388, 8, 51, "Input",ExpressionUUID->"ca284c07-f2f6-41ad-8388-0ee0d89debcf"],
Cell[21536, 620, 779, 24, 51, "Output",ExpressionUUID->"83c90568-ee4d-4a92-9ab4-91c05cf4f791"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[22364, 650, 572, 8, 87, "Section",ExpressionUUID->"24cd3001-ea45-47ab-bb96-5c6765007add"],
Cell[CellGroupData[{
Cell[22961, 662, 418, 7, 51, "Input",ExpressionUUID->"e6436050-4ae4-4e7b-ba6d-22cd79f0763c"],
Cell[23382, 671, 432, 12, 51, "Output",ExpressionUUID->"09110649-9296-4262-b952-8afd8ec53130"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23851, 688, 473, 8, 51, "Input",ExpressionUUID->"969a80a6-e6b9-4a50-bc52-e27f08dc7073"],
Cell[24327, 698, 910, 29, 51, "Output",ExpressionUUID->"648b9271-c75f-42f8-ad44-88fffa7aa3b8"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[25286, 733, 430, 6, 87, "Section",ExpressionUUID->"72b4bc04-7676-4749-91f2-3fd8b4907447"],
Cell[CellGroupData[{
Cell[25741, 743, 578, 8, 35, "Item",ExpressionUUID->"8aa5c24f-f858-4959-8683-e9356406a9f1"],
Cell[CellGroupData[{
Cell[26344, 755, 1669, 36, 104, "Input",ExpressionUUID->"f56a7601-52fa-4249-8237-485543004148"],
Cell[28016, 793, 10009, 181, 468, "Output",ExpressionUUID->"313a5eee-569e-40b1-a654-e775860f404d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38062, 979, 1191, 29, 104, "Input",ExpressionUUID->"c7e86c4e-3b40-43e0-b4ff-f51cde07196a"],
Cell[39256, 1010, 6324, 114, 468, "Output",ExpressionUUID->"b168c029-3e40-4119-bb00-4ee6276a91ae"]
}, Open  ]],
Cell[CellGroupData[{
Cell[45617, 1129, 1142, 29, 104, "Input",ExpressionUUID->"068e9cfd-920e-4a4b-8c13-9ff540e8dd83"],
Cell[46762, 1160, 10253, 188, 468, "Output",ExpressionUUID->"853e50f5-f1a6-4178-9cc5-1de55872be0c"]
}, Open  ]],
Cell[CellGroupData[{
Cell[57052, 1353, 238, 4, 51, "Input",ExpressionUUID->"667ba917-4317-4dec-b2df-542439b636e8"],
Cell[57293, 1359, 23823, 409, 468, "Output",ExpressionUUID->"0cc1c95b-36c0-4b1b-9751-bbb2e00791a6"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[81177, 1775, 421, 8, 135, "Section",ExpressionUUID->"08f23630-6205-47fe-bef3-851508e50197"],
Cell[CellGroupData[{
Cell[81623, 1787, 593, 9, 35, "Item",ExpressionUUID->"dfe2d88b-f43b-4b63-b949-281ef52c42cf"],
Cell[CellGroupData[{
Cell[82241, 1800, 2460, 60, 182, "Input",ExpressionUUID->"60605db5-98e4-425b-998f-80caa384fb53"],
Cell[84704, 1862, 1377, 35, 104, "Output",ExpressionUUID->"8eac1b0e-f123-4191-95fd-d15f9f7f0272"]
}, Open  ]],
Cell[CellGroupData[{
Cell[86118, 1902, 298, 5, 51, "Input",ExpressionUUID->"c6a2e8ce-9921-4772-9e4c-9aaf537add02"],
Cell[86419, 1909, 344, 5, 51, "Output",ExpressionUUID->"f1794859-77c8-48d4-88e3-e4d85584e52f"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[86824, 1921, 319, 5, 87, "Section",ExpressionUUID->"854c782d-2b8b-444b-a73a-b526a1e58fec"],
Cell[CellGroupData[{
Cell[87168, 1930, 1316, 30, 104, "Input",ExpressionUUID->"5e9f4592-8b37-4f66-b7f2-a07dcded9e65"],
Cell[88487, 1962, 2123, 41, 645, "Output",ExpressionUUID->"c243464a-18b0-45b2-a6c0-f99eff34106b"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[90659, 2009, 222, 4, 87, "Section",ExpressionUUID->"6e04284b-303c-4a0c-842f-7bbae89eb16b"],
Cell[CellGroupData[{
Cell[90906, 2017, 619, 13, 51, "Input",ExpressionUUID->"1ae1de4c-6ed7-4458-8a86-d2722aeea5de"],
Cell[91528, 2032, 653, 18, 72, "Output",ExpressionUUID->"d30075e1-90ce-48f7-a9c5-d056d0302794"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[92230, 2056, 275, 4, 87, "Section",ExpressionUUID->"d1a22f61-c8d3-47a1-9758-bea908fb52dd"],
Cell[CellGroupData[{
Cell[92530, 2064, 548, 10, 51, "Input",ExpressionUUID->"b5ec988d-29df-4e49-b816-498bf993d405"],
Cell[93081, 2076, 2136, 60, 138, "Output",ExpressionUUID->"0c17ad1a-0f92-4b11-973e-7a3e8a837c04"]
}, Open  ]],
Cell[CellGroupData[{
Cell[95254, 2141, 297, 5, 51, "Input",ExpressionUUID->"68341285-01db-4a82-ab17-14188cd97bd0"],
Cell[95554, 2148, 2091, 63, 86, "Output",ExpressionUUID->"97402d29-f18b-49b4-a545-606ba8e61535"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[97694, 2217, 330, 5, 87, "Section",ExpressionUUID->"f80e7fa8-76db-458f-84cb-46a26f53cc60"],
Cell[CellGroupData[{
Cell[98049, 2226, 542, 10, 51, "Input",ExpressionUUID->"c8cf73a1-4ff4-48c6-8851-4b7969d6d44d"],
Cell[98594, 2238, 1637, 53, 104, "Output",ExpressionUUID->"c6ff30d0-d33a-478e-9f5a-3888f2dca2f7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[100268, 2296, 271, 4, 51, "Input",ExpressionUUID->"abefd052-73e1-4d38-83c5-06a05ccbdaf1"],
Cell[100542, 2302, 2649, 79, 121, "Output",ExpressionUUID->"30839153-fd82-4a63-a016-e74593002774"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[103240, 2387, 366, 6, 87, "Section",ExpressionUUID->"e39c5931-fac6-4efb-beb0-441a6e3d99ea"],
Cell[CellGroupData[{
Cell[103631, 2397, 446, 8, 51, "Input",ExpressionUUID->"ce0bb349-5cb2-4f0b-83a4-be5de002e619"],
Cell[104080, 2407, 1022, 33, 78, "Output",ExpressionUUID->"56b0a5e8-cffb-4f25-8eea-17001a3fd636"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[105151, 2446, 572, 8, 87, "Section",ExpressionUUID->"5baf18f4-4931-481f-ad99-f5aaddf7e5b1"],
Cell[CellGroupData[{
Cell[105748, 2458, 479, 9, 51, "Input",ExpressionUUID->"c28a2f38-4007-43cf-b4c9-f755f30d66d4"],
Cell[106230, 2469, 878, 28, 51, "Output",ExpressionUUID->"4f37fb08-8ea6-403d-b301-89ecf38f4878"]
}, Open  ]],
Cell[CellGroupData[{
Cell[107145, 2502, 534, 10, 51, "Input",ExpressionUUID->"5517b3df-8a28-4367-999e-17e479c8c77a"],
Cell[107682, 2514, 927, 29, 51, "Output",ExpressionUUID->"fb871d7b-af88-4454-9324-07fd9b1edabf"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[108658, 2549, 430, 6, 87, "Section",ExpressionUUID->"1badad71-e6c0-4744-b8c1-3bf58078cf37"],
Cell[CellGroupData[{
Cell[109113, 2559, 578, 8, 35, "Item",ExpressionUUID->"75a5b281-a86d-4e96-91f2-c41dbc6d5f85"],
Cell[CellGroupData[{
Cell[109716, 2571, 1907, 41, 130, "Input",ExpressionUUID->"033cc442-147b-4f00-aba2-3b966159717d"],
Cell[111626, 2614, 10083, 182, 468, "Output",ExpressionUUID->"e58b5636-c999-4422-9be0-5f425259cbbd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[121746, 2801, 1326, 34, 130, "Input",ExpressionUUID->"ce31771c-ea40-40a8-a381-a1f411569b3c"],
Cell[123075, 2837, 6350, 115, 468, "Output",ExpressionUUID->"a95764da-9008-4c4b-bca9-377f52eb86b6"]
}, Open  ]],
Cell[CellGroupData[{
Cell[129462, 2957, 1298, 32, 130, "Input",ExpressionUUID->"07ec3a20-4aa2-4a61-9ffb-2d864700b0e4"],
Cell[130763, 2991, 10301, 188, 468, "Output",ExpressionUUID->"fa5c03e0-9faa-4703-8d71-88426a8e8bb2"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

